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Chapter F2 J . . 

, , ' I ^ 

INPUT-9UTPOT- AND ASSIGNMEOT^STATEMENTS 

• . f . J . ^ . « • 

F2-1. Intro.giuction '^ 

In Chapter ^ V? developed an appreciation of input, output and assignment 
stej)s as components ^of algorithms expressed in the form of flow charts, ^o 
^ far, we hav^* viewed^^flow irharts as a means for conveying a sequence^pf corapu- 
Wion rules primarily Trom one person to another. We have tacitly assumed'. 
. . that only man c«j read, underst^d and carry out the intent of such flow 
^^charts. Naturally ve want t^^ include computers in the set of^all thi^igs- 
which p^n fead, understan^^and-'carry out procedures. . , ^ 

K)RTRAN II"- I^nguage and ^ processor * ; ' ' * , 

'•Fjogrammi'ng languages like ALGOL and FOROiRAN accomplish .this objective. 
The steps of a \p|t>gramniing language are Called statements . They correspond 
roughly the boxes of a flow chart, ^ 

, ^ ..FQRTRAJf, meaiiifig' FORmula TTRANslation^ _is an English-like prograHmiing 
language developed 'by IM about ten years, ago. It 'was hf storioally the 



FORTRAN ll is an improvement^over 



fipst 'o*f a se3:ies of slpiilar languages 

the original vfersiop and is* probably the most widely used of available pro 
graimaing languages. More advanced versions are aifso in common use. They 
' ai:e known by varipus names, the most widely used naitfe' being FORTRAN IV. 
. _ - FORTRAN II was designed with th'ese- objectives. 

1, A'^T/ide-variety of algorithms can be described with this language, 

** . f"^^ chief area of application is for expressing algorithms which 

' 'j ^ ' '"^'deal with scientific and engineering computation. Algorithms ;for 

many bther types of problems can also be expressed satisfactorily 

/, . '^in FORTRA^I II..* 

^, *- ' ^ 

I ' 2. ;The rules ^or "grammar", of the language' are defined precisely 

-\ (jjnlike English). ^The h^l result' is that an 'algoritll^ written^in 
\ ' \ " FORTRAN II c&n mean precisely the same thing 'to each person who 
_ A % -reads it^ i.e., It means the same thing ^ eve ryvhepre to people* who 

\ ^^"^S; ifamed 'this lan^ge,^^ Hence, FORTRAN II provides\a means for 
communicating" algorithms' via Ncorrespondence and publications, from 

... 1 



on.e person ^to- anther. 
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With minor additions or modifications,. VoRTRM II can be ""implemented 
on many type^ of digital computers^ By Impleinentation we mean tHat 
a computer can "be progra™e^^to accept a^orithms as "source" pro- 



grams that are writfer^-'i^h a^xrersion of FORTRAN II and automatically 
convert them to sequenoes^'bf con^ut^er instructions, often called 
"target" programs, wh|fch can then bfe executed by the coppu^er. 

These features of PORTRAC/ 11 -have led to it's wide acceptance Wong 
imthematicians and seientistsT In 'th^s chapter, we shall' "begin to describe 
FORTRAN II as it is typical w implemented for use on a computer, calling it 
simply jtiRTRM. 



1 



Each state^tent in 'FOmJRM is yritten s.o that when typed 6-r'. punched on 
a -^ard it can be trans fer/ed to tme computer's memory. Here it can be 




anned character by cha/acyer '/and analyzedfoFT^^ftrl l Ir i Leu tv^ y - 

Programs which anajiyz^ t)4ese stateme^its are called compiler or processor 

I J. J " ' ^ ' " 

pifograms. A typical. M)RTRM processor reads statements originating on punch 

^ * / / . * * 

cards and analyzes thfem b/. converting each statement Into an equivalent 

sequence of c omputerp/st rue t ions • W4r^ these in*structi'ons to be executed,* 

it would amount to having the computer carry out the intent of the statement • 

which is itself distinct from the machine's particular injtruction repertoire, 

' • „ ? . • 

Target programs and source programs 

I" * , * *■ ' * 

^ The processor program will read, and analyse ' all ' the statements of a 
FORTRAN Program to generate a complete set of instructions, or "target'* 
' program, sufficient, if execyte'd,' to, carry' out the intent' of the^„ entire 
process. The "target" progra?! gets its name becausoi' it is the target or. 
objective of the processor program. In turn, the processoi* has received as - 
input a "source" program written in th^ K)RTRAN languag^'^i-^w*::::.ij|^^ 
prograH| fin^ishes geaetating the' target psogram, the computer* may execute ''^-^'^ 
these instfuctions fight away. This' is .feasible because the target program 
is developed and kept in* the computer's- memory. When the target program is 

/ ' * . ^ : 9 . . « 

too big to fit in high speed memory along with the large processor program, 

- -> ♦ • • 

^ -the jgqnerated target is stored ^temporarily in Some form "pf* auxiliary iliembry 
. ' ,1 ^ * . ^ ' 

media sucTi a^on magnetic drums, discs, or tape^, 6r in the form of decks 'cyf 

$)unched cards. When stored in auxiliary memory med^a, eapecialTy jJinphed 



cards, target programs ^can be recalled for execution at any subsequent time, * 
as suggested in Figure F2-1. We would rarely wish t9 read or ^itudy the 
'target' pro-am burselves/but in. principle this can always be done by causing 
the processor to print the target progran^. ' • • ' 



■TARGET PROGRAM • * 
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/processor \ ' 
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SOURCE 
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ON CARDS 

STORAGE SUCH AS 
TAPE. DRUM OR 
DISC 



Figure F2-1. The "compiling" prpcess y 



OK ON PRIMT^O 
PAPER, 



It may be intriguing to you to i^earn how the processor prograin^ dpes its 
Job. After all^, it is also a flow-chartable process ^and hence could easily 
be wit^n our abilities to understand it. Chapter 8 wlU^^shed some Jight on 
this interesting process. For the present^ however^ we will avoid any head- 



on discussion of thife'^-tS^pic because our first ^inlyerest must be to.learh^to 
irrite simple algorithjns for* solving mathematical problems in FORTRAN* We- 



will^, however^ be making occasional comments -which* bear .iudir'bctly' on 'the 
nature and structure of the processor. * 



\ 

V 



General at^pearance of ^a' FORTRAN proRrapi . . n 

' Recall JbHe process for computing * ' ^ « 

D VA + B .+ C . , . 

^wh^se flow chart ve displayed in Figure 2-2. Each l^ox can be writte^n as a 
FORTRAN statement as shown in Figure F2-2. ' ' ' } 



'Flow Chart. 



' , FORTRAi? statement 



input^ 



agsignment 



A,B.C 



I 




output 



I 



100,. A, B; C 



D = SQRTF (A^2 +' ]^'^2 + C^^2) 
PRINT 101; A, B, C/ D 



Figure F2-2. FORTRAN Statements compared with Flow Chart Boxes 

Notice the 'similarity between the^boxes aA the FORTRAN statements. The 

differences^ are largely superficial /that is, the ISrmbols used in eacfi caee 

may be different /'but the ideas appear to be the same, 

♦ ' • ^ 

We-4oii*t have to coniject the statements with arrows, because, yhen we 
write FORTRAN statements one below the other, we imply ihat they are *to be 
carried out one after the. other from top to bottom. In order to suggests,^ 
repeating the process »for many -sets of data, we drew a line J!rom Box 3 back 
1^ Box 1 in tXiS flow chart. There is an* analogy., in FORTRAN to accomplish 
t^ie^, ^me'' 6b j^t ivV. We» can simply give the READ statement; a numei^c . label , ^ , 
'corresponding to 4ts box nujnb^r in the flow chart, which is 1 in this 
case. Then we aid after the PRIiJT a statement which sends control to' the 
designated Statement. .This is shown in Figure' F2-3, • 




LABEL 



STATa^Nl" 



READ 100, A, B, C ^ , 

• *" • 

D = SQRTF (A^-^2 +' B^^2 + 0*^2) 
PRIl^IT 101, A, 'b,)C, D 

/ * 

GO TO 1 



Figure F2-3. The GO TO statement directs control to any desired -s tat ement. 




.<N-;\.-r:r>^^ Other vo¥^a, we give the READ statement the label 1 and then 
"!ace a F^)RTRAf^ -statement^ " » ' . 

GO TO 1 

J'or^the purpose of indicatin^^ a .lump or transfer to that stateme;it. The jump 
statement c?onsis1?s of the symbol GO TO followed by. a numeric label. Because 
it reads like English, the jump statement is easy to understand and we shall 
^ V say no^more about it at -this time. ^ - » . 

As we focus on an gxamp^)^ of a new , language, it' always tahes a while for 
its features, its special symbols and-punctuatio^n patterns to stand out - 
clearly. V Y©u have probably^observed, and co^rrectly, that "fn PORmAN 

1. Statements are separated by virtue of being written on •* . 

separate lines. ^ ^ - ^\ 

'2. Alphabetical characters in each statement appear to be' 
, exclusively upper ease. / * 

* * ■ ' \ 

3. A statement may have a label. If so,^the lab^L precedes the 

statement, and. is separated from it 'in some way,, apparently . . 
'by some blank space. 

h. The as^signme?!^ symbol is the eqoial sign ( = ), an unfortunate 
« state of af fairs^^^ince a left point Jng' arrow wou,ld ba,V:e ^ 
V ' be^n our choice? ' 

5* ThF^symbol for Exponentiation is the double asterisk (-^*), and 
we appear to have ^lost the ability to use a symbol like '7 , 
having now to use SQRTF which is a curious abbreviation'. ' 

^ 6. Certain nijimbers, specifically . 100 aod 101, which appear 

following the words READ and PRINT, have no apparent 
^ . counterpart in the .flow chart*^ * f • 

Now, before taking a more methodical look at input, output and assign*- 
menf statements and rules for forming these correctly, we show in Figure I2-k . 
a complete FORTRAN program as it might be written on a coding sheet prelim- 
inary to-being punched on*"cards, one card^per line of FORTRAN coding. 

'We can now notice a number of features. Some have to do with layout of 
the^cards. Others relate to the language of FORTRAN itself. 

Card layout ^ . ,^ 

A statement label will ^ be punched in card Co).uinns. 1 through 5. The 



ERIC 9 



1 



statement; proper will te placed in card Columns 7 .through 72, We can'add 
comments to identify our program or to help make it self-explanatoiy; 
comment will be recognized if we place the^letter C in card Column 1. So 
the title, */ " 
> . ' . * EVALUATION OF 

is treated -as a commdnt instead of a procedural statement, because tl^e letter 
C will be ^punched in Column 1 of this card. We don^t expect the computer 
t(^^ay any attention to comments, They* ^re there mainly to help u^ document 
what we are doing.' . ' ' . 

By crommon convention. Columns 72 - 8o are not considered part of the ' 
program. You might say the computer. would ignore any g^rt of a program 
which happens to qe puached in these columns. 



o 



•o 

FORTRAN .CODING/ 





' ' / 


Progrom 


Grophic P 








Progjromnr^^ec ^ 


Dote ' 


Punch 









isTATEMEHiT 
NUMBER , 

> '5 


c 
Q 
VJ 

6 


K . I • ^ \ 

FORTRAN STATEMENT 

7 10 >5 ' 20 25 . ^30 '/ 35 1 » 40 > 


c. : . . 




i EiVALUAiT 1 0N ,0F D , , , > / 


' ■ 1 




READ, 1 00,, A B,, C , ,,.•,/ 






D = ,S ORTF,( A^(^t2, + B#,*.2 + ,C*^t2), '.X 






P.R.I N,T. .1 0.1,. A^.'JB . CD..,.,-. . , \ 






G0 T,0 .1. , . . V ', • , , . . ^ , \ 


. 100 




f0.r'm,at.(3p,i b'-s.)," .y 


i:o.i 




F.0.R.M,A.T.(.4.F,I 5'.,8,), ', .'. .,..,;;./ 






tM.D. ,.<.,.,...,,,.■..,..,.,.,..(' 






. . , 1 . . , . I . . . . 1 . . . , 1 . . 1 . . . . ^ . . . . 1 \ 



'Figure F^^h, Compl&te program written on a FORTRAN gDding sheet 
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Column 6, although not utilized in the example .given 'in .Figure F2-k; 
will have ^ special significance to be explained later. • ' 



Language 



Notice that a statement 



END 



is placed at the end of the program to m'ark the end; Thi^ is charact^eristic 
of al^^ FORTRAN programs. They. are not complete without this terminal state - 
.ment. • - ' . . , 

There are two n^w statements shown labeled 100' and lOl, Statement 
100 provides a format required by the, READ statement. By format Ue mea^i 
certain coded information pertaining to the Vayout ot da^a as it might .be' 
punched on one or more data cards or as it mikht; appear printed or typed out 
on one or'^more Unes of paper. ' Like^S^ise statement 101 provides a 'format * 
required by the ^ PRINT statement'. We wtir defer further discussiorf of 
fonriat codes until Section F2-j.- (if ybu are wondering what special signifi- 
cance the formal numbers 10t3 and 101 have— don't bother'. You "will see 
later that these are^qhosen quite Arbitrarily.) 

So far, then, we see that a FORTRAN program appears to consist of a 
^group^of st^ements and'possibly' comments eading with END, £;a,ch ^t^iteraent ' 
is punched on a card of standard layout so Coding' fo^s like that shown in 
Figure F2-i| may prove helpful, A picture of the cards produced from this^ 
coding form is shown in Figure F2-5. 
' / 



OQUcea rr< 



/ 



GO 70 i 



R£AP too. A, 3, C 



£\<4LaATiON OF P 



COL,S C0L.7 \ 



C0LZ2 



Fifewre F2-5. The Program as"'a "deck" of cards 




fi2-P FORTRAN language elements 

In this' section we shall deal with the character' set or "alphabet", of 
FORTRAIJ and their use in -tion^tructing some of the simple componehts of state-^ 
merits, specifically the construction of numerals (constants and labels), 
variables, names for functions/ and*-operators . ' . , * 

]e jhara-^ter s^t \ ^ " ' ♦ 

t • . • 

TJie char^ter|^^ich are used in FORTRAN statements are limiteci to 

those shown Table P2-.. . Notice the absen'^e of lower case letters. » 

■ • V, ' '' ^ . 

, . ^. Table F?-l 

^ ''^■/^"TlTenPOtlTRAN Character Set ^ , 

(a) Letters- ABCDEl^FGHI 

J K L- M N 0 Q R 
. ' • • S' T U V W .X Y Z ^ * ^ 

^ (b) Digitsr- . 0'1?3.^^56789 

^ ( 0 "Special Cb^araeters + - ^ / V ^ 



Constructing num^tal s (numerical constants ) and labels 

To understanci how numerals are written in a FORTRAN statemenlT it helps 
to re -all the two scheme? ve discuss in Chapter One for Representing numbers 
interna 1 1.V , i.e.*, it\ a word of rrte^mory. TT^ere^was th^ xntegei' scheme and the 
floating point . Scheme. Externally, i.e.,, in a FORTRAN staterrtent, there gre 
two types of numerals, called integer and real. v 

For example, the numerals 

h, +i-f*, -15. 1001 
' ' \ ■ , \ 

Q.re0^1 examples of integers in FORTRAN. ' • 

The connection between the external and tj^e internal » representation is 
this; If a numeral, U, for instan^?«, is written as 'a FORTRAN integer 
(extel^nai^y ), then it will be represented internally in the integei* scheme . 



52-2 



Expressions involving FORTRAN integers are then evaluated in integer arith- 
metic which is, of cpurse, characteristically different from real^ arithmetic 
as -explained in Chapter One. • " < ' * • ' 

, . >' *' ' ' * ■ 

If a numeral; kj_, for instance, is written as a FORTRAN real, -then it 

will be represented internally in the floating point scheme. Expression's 

which involve FORTRAN reals are then evaluated using real arithmetic .'"I" « ' . 

- We' see then, that the way we write numerals as numerical constants in 
FORTRAN statements (i.e., with or without the decimal point) determines by 
inspecti6n, and without ambiguity, whichrtyt^ of •internal representation the 
number is .to >hav$ duying .commtation. • . ^ 

There is one more, way we can write a FORTRAN real. This is to u^ the 
"E" notation mentioned in -Chapter One.. You will rejiall that "E", as it 
xs used in ^ numeral, is merely a convenient symbol whose meaning is "times 
. 1:eh to the power of". Some examples are: ^ ' ' >^ 



Number 

52 X 10"^ 
5.2 X 10' 

-2 



6'. 023 X 10^^ 



Method 1 

(Without "E"T 

.00900652 - 
' \ 000000052*^ 



Method Z 

^- ^(With "^E" notation) ' 

' 52.E-8 
5.2E-8 

-2.E0 
or -.2El^ 
or -.2E+1 ^< 

or 



602300000000000000000000 . 



1. 587^011 X 10^!^ 



1587^10110000. 



29.E-I 

6.023E23 
fAvogadro^s number^, 
the number of mole- * 
cules in a gram mo]^e 
of gas at standard 
temperature and 
pressure. ) 

1. -587110^10, 
(1 light year ingles.)' 



The "E" notation provides Ois Vith some flexibility. It is especially 
useful in expressing number^ which are v^ry small or very large in magnitude. 



'Strictly speaking, the arithmistic used is floating poj;nt whiah, you. 
recall from Chapter One; waB 'said to be a way of 'approximating true 
real arithmetic. We take the liberty of . calling this type of arith- 
metic "real". / • ^ ' , c ^ ' ' 

> / • 
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Moreover, any one real number can be expressed in a variety of ways. The 
decimal point in the first part of the numeral c^ be shifted at will without 
changing the value of thfe number, provided a corresponding change (in the 
dpposite sense), iia^^de to the exponent. Far example, moving the decimal 
point to the left is equivalent to multiplying the number*by a powei* of ten, 
so we have tp subtract the same power from the exponent. 



Exercises YZ-Z Set A ' ' • . 

1. Show three ways (other than those showft* above) to. write -2." in FORTRAN. 

2. Use "E" notation to show 1 light year in centimetjers . Use a precision 
~^of \ decimal places. _ / 

3. Write in customary notation the values represented by the following 
FORTRAN constants: . * . ' ' 

5.$6eo, ^002E+8, 88.E-1, * -522.4, 5.E5 ' 

V 

V/rite real constants usin^ "EV notation for the numerals ^ive.n: . ' 

3.91 XJLO'^ • 

9.09 . ' . 

^ , .6.67 X 10-^ ^ ' ^. . * . 

176. U X 10"^' . ' < ^ 

5. 1)0 all the FORTRAN, constants given below on, the left correctly .express 
. the^uggested vialues sltOwn on the right?, ^Correct those FORTRAif constants 
that do not. ^ 

/ • ' 

179E-i^ 1.79 X 10/^ . ' ^ , 

. • 6179. ErT 6179 X>iO"? - 

». • I 

, l6.'79E 16.79 ' 
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Statement labels ^ * - • , * 

In /FORTRAN, statement labels are special ^ forms 'of PORTIiftlN in^tjegers/ 
namelj^, unsigned ir^eger constants other than zero. Thus ih,^ 1000, 393 
would be perfectly good labels,;; but -h, 0, Ht/ pr SI are not , acceptable . 

Constructing variables ana name^^for functions • 

A variable is formed according to- this rule^: 

^ A letter, or a fetter followed by a limited sequence of letters or 
digits^ Tfie total number of characters in Ihe name is lirilited.* Most V 
implementations (FORTRAN processors) perniit^ variab^^s to have up to six , , 
characters. We will assume this limit of six characters for our text. 



Examples 



HARRY TEMP 
. IKE COLUMN 



X Tii6 
A15AA 



of number represented *• ^ 

What tyi4 of nuri^jer is represented by a given variable?. In the FORTRAN 
expressions *ljke T ' " 

* A + B or A / 6 

do you think values for A and B are represented internally in the integer, 
or in the float'ing i)oint scheme? It makes quite a difference, because the 
kind of ^arithmetic which will'^be performed in evalfuating A + B or A / B 
depends on the internal re.presen-t?ati^n o'f A and B. 

Consequently, a rule for determining the mode -or type of internal repre- 
sentation -.is required. ^'It is this: ' ' ^ 

> All variables which- begin with the letters^ I, J, K, L, M or N - . 
correspond to numerical values which ^ are stored in Fhe integer scheme.. We 
*sh9^ refer to such variables as ^integer variables. ^ ' • * 

AU other variable names then correspond to variables sto?'ed in t^e 
floating point scheme and can be referred to as real variables . 




Examples of integer variable names are 

ISTEP3, 11^ ITEM, KTl, .NUM, -KCOUNT 
Examples 'of real -variable names are • 

IMP, X2, FNUMBR; SMELE, ABEL, BAKER 

11 

15 



r 
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Exercises F2-2 Set B 



1. WKich of the following inscriptions (or character, strings) ^are' invalid 
for variables'? * - ^ 

ZZFfZ, 19763; '^•51976^ ' 19S76; NN, 

ITIF, VAR, T15TT; - MAR.Y, MAR*Y - * 



( 



2. Which pf the followi^ character strings are valid for integer variable*? 

' ..'I15Z., N-G, J97r G2tT; ABI,- ' ^ 



^37; N37^ 



Names of predefined functions * . n 

'Several kinds of functions ma^ Be named ia FOHTRAK. In" this \;hapter we 
shall use a group of predefined functions. The rule for naming these func- 
tions is similar to the rule for naming variables except that all names for 
predefined functions end' in the letter F'J For example;, 

. ABSi' for absolute value function '\ 



■ V 



and SQRTF for, square root func/tion.. 

Other functions of this kind are given in Table F2-S. Like variable* ^ 
nameS; these^ function names must begin with i& letter ^nd must be followed by 
a limited sequence of letters (ending in ^e total number of characters 

in the function name is not entirely standard among the various 4)rocessors» 
For this book we sjiall asfeume a predefined "function^^Bie consists of from 
^ to 7 chai'acte^i^s including the terminal F. The letter F, by itself^ is « 
not a legal name for a function. ^ . * 

In. order to avoid confusion "between these function names and variable 
nameS; we shall adopt this arbitrary, but very |)ractical rule: 

i ' - 

Never end a variable name with the let^tet* F, if the name _is to 

consist of four or more characters . Thus, ' ^FF, TMFF, and JEFF would 

be avoided, but TIF, TUF, JEF and F ai^e perfectly '^'OK" para^^s for _ 

variables • ^ ^ ^0}, 



.0. -.- • . . J6 
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.Table F3-2 
Pre.defined Mathematical Functions 



Name 



ABSF 
SQRTF 

logf"* 

EXPF 
SINF 

COSF 

atAnf 



Meaning 



absolute value » 
square root 

logarithm to the base e 

povers of. e or exponential 

si^e of an angle whose measure is 
- given in radians 

cosine of an angle whose measure 
is given in. radians 

arctangent or principal angle in 
radian? of a given tangent 
'vaKIe. That is, ATAN(X) 
gives a value in radians 

. 'corresponding to'the principal 
angle whose tangent is X. 



^^pst of the functions listed in. Table P2-2 should already be familiar 
to y6u. When we use one of these special pre(lfefined functions in a' FORTRAN 
statement, the resulting machine code .automatically carries out the evalua- 
Son. of the specified function. The ^rgunjents we use with the functions 
listed^^ in Table F2-2 are all real and th^ function values ti^t are developed 
are also all real values. ' - * 



Operators ' ' ^ s , 

To write arithmetic expressions and assignment statements need 
operator symbols. Considering the limited character se\we have with FORTRAN 
it should not be surprising to find some, compromises with conventional matl^e- 
matical notation. Table F2-3 sh9ws a' list 'of the symbols we shall be using 
for the varioxis arithmetic operators. Fbr convenience in later discussion 
they are given in hierarchical order, that is, in descending order of, pre- 
cedence which is, the same in FORTRAN as ^t is in our flow chart language* 



We have also included the assignment symhol in this group. It is a "binary 
operator^ hut^ of course^ not really arithmetic' in nature. 



' ^Tahle F2-3 ^ ' 
FORTRAN Operator ^hols 



itoRTRAN 



• \ 

\ 

same level ( ^ 

of p^-ecedence \ / 



same level 
of precedence 



A- 



Meaning 



exponent iat ion 
raising to a power 

multiplication 
division 

stddition 
suhtractiqn 



assignment 



Example 




A**B means A' 



A*B 
A/B 

A*+ B 
A - B 



B 



A = B means A <~B 



r . 
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F2-3 Input - output statements 

^Now thatr'^ou hav% become somewhat a'ccustomed to the appearance of FORTRAN 
characters/and 'to their use in the elementary components of the language, 
like nimierals,^ variables, function natnes, and operators, you are ready to 
study/the three important statement types: input, output a^id assignment. 

% 

'You have* already seen examples of input and output statements, namely: 

^ f, READ 100, A, B, C . ; 

and '"/'*" 

, . ' PRINT 101, A, B, C, D 



■The Statements are simp^ in form. They always begin with the key word , 

* : ' " READ , # 

' ^ ' . 6r PRINT 

I ^ Following the key werd we provide aif unsigned integer constant which* is 

the label of another jPOgTRAN statement. This label tells where to get some , 

needed fontiat "information. We call this label a format' number. A comma 

' , rr; 

follows next. «Then ve write the list of variables whose values are either 

to be read or printeii, This list is simply' a set of variables! separated by ♦ 

commas, if there is mo're than one. 

s • ^ J 

We will not limlit the number of list elements of an ingut or output list, 
but we will prohibit pertain things from beirlg on the list, k constant for 
example. . > 

Thus-, ' ' ' J ' . \ 

^ . READ 100^ 2.5, ' A ^ 

•* or PRINT 150^ A, B, 2.5, C,.152, D 

are invalid , because constants are prohibited as list elements. In the sam^ 
vein any expressions which require evaluation, such as A + B,^ ar^^ such ^, 
f orms 'as (A), (-A), +A, or -A, are prohibited. ' o. ^ . ' 



Thus in the two statements 



If 

that follow 



READ ' 15c, A + B, -A, +B 
or PRINT i51, ^*A + B^B ' ' , 

each list element is invalid.* ^ ^ 

In this ^ay vp arrive at the "general form of an it^t %t output state- 
ment, i.e., 

^ . READ format number , input list 

or PRINT format number , output iist \ 
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For a REAP statement, the list elements are the variables^ j^hose values 
are tQ be assigned from th^ input data, while for a PRINT -statement the list" 
elements are the variables whose currently assigned values a'^e to.be printed. 

< * ' ' ' ^ \ * 

Executing a READ' statement " . * 

~~~~~ ^ . ' , 

* ' * ^ * ' • , * 

In this disciiission we shall assume, as in the flow chart text that inpi^t 

data originate on punched cards. 3?he effect of executing' a READ statement 

is as follows. • . ' ; . , t 

1. Fir^t we assume a card is in position to be'read by tHe computer's 
input device. If not, the execution of the program ceases 
immediately. In some but not all systems, the Qomputer might then 
print some message like ' , ^ f 

i "YOU HAVE RUN OUT OF DATA" - ^ \ 

or "ALL INPUT DATA HAVE BEEN PROCESSED" 

/ «^ 

2. The format number of the READ statentent is then used to jconsul^^_^ 
(look up) the format code stored in memory which has this formal^"" 
number. The jformat number is simply .a way of identifying the 

* particular Fbrmat code you ^spt among several different formats 

that may be employed in the same program. 

FoAiat g/Sie^ yhen properly ititerjoret/ed, will, in effectf, tell the 
computer where on the card to find thes desired -values. The, code 
itself is written in a special language which will be*- described 
momentarily. • • . * ^ 

3« Th^ content of the card that is ready to Toe read is then transported 
to the computer memory, where it is examined one character at a time 
^. , beginning with Column 1. The format code whiqh was f9Vind in Step 2 
is then used *to decide which numerals (from>the card) go with-^hich 
r variable (of the input list). The purpose, is to^aqijieve ^a orie-to- ^ 
_ .one matctjing between .the numerals on.lhe^a^rd and the variables pf . 
^he input list with which they are to be&ome associated. 

At this point some examples should help us see how this matching process 
is cari'ied out, , ^. 

Example 1 

Study JS^gure F2-6 where you see a READ statement, an associated FORMAT 
statement, dnd a picture of a card that might possibly be read as a result 
of executing the given READ statement.. 



The format number is 22. THe associated format code is (115; 115; 115). 
The CQde is interpreted to meaf^-in thisjcase--that: the data are^ arranged • 
on the card in tkree group^s/ each fifteen columns .-in vidth. .The* letter ^ I, 
for integer; tells us tho/t each group of 15 columns is expected to contain 
an integer which can stored in memory at a place associated with an integer 
variable. 



READ 22, WUM,KPAY, KAMT 
ZZ FORMAT (ri5;II5,II5) 




Figure F2-6. A data card with three intege/s punched on it; 

together with related READ and FORMAT statemfents. 



, We shall call each Q\},§h group of columns a "field. Codes like II5 are ' 
then called field codes because they tell the type of number; I; (for ' i 
integer) in this cas^e; and the width of th^ field in wl^ich the numeral can \ 
be found, • - , . ^ ^ 

You may wonder if there 'is a,ny special reason for choosing a f teld width- 
of 15 columns for -these integers. We want a field which iS; of course; wide 
^ enough to contain the largest expected integer value. But ^it doesn'-t hurt to 
use a wider field if room is available on the data card. The main, thing to 
remember; howe^^er; is that the unitarportion of the data value should b^ 
punched irl the rightmost position o^the card field. This keeps you out of 
trouble because; in many computer systems; any blanks pside the field^arid to 
the right of the numeral will be interpreted 'as trailing zeros. v'ThuS; if the - 
integ'er .1502 were .punched in colunms 9-12 instead of 12-15; on the card 
shown in Figure F2-6; it would be erroneously interpreted as ^^1502000. ' 

j ^ As the scan of the data from the card proceeds; column by column; the 
information is picked off ohe ^ at a time to be matched with the next 
^item in the input list. When a match -is achieved; then that numeral from the \ 
card is converted to its internal form as.' a sequeace of bits representing 
^ ^that number in tlie appropriate ^cheme. This is done by a special program 



• F2-3 ' , ' , • , - . 

provided by th^ prog^ssor. For I- fields the conversion is made to the 
* integef -scheme . 

Several general rean^rks are always true concerning /orma^. st'^atements . 

y 1. 'Every format statement in a program jmist have a •(uniquej. statement 
' * number. * * ^ 

2, Every format ^tateAent has .t he form FORMAf (format code). 



Example ''2 



^ Nov study Figure^ F2-7 where you see another READ statement al6ng with an 
associated FORMAT statement. Here we show a data, card with numerals for three 
real numbers punched on it. The field^cod^s, which show where these numerals 
8^'re to be found and how they are to be converted to their internal form,' are , 
each FI5.8: , ^ ' 




100 F0RMAT^i5.8,Fl5.8/l5.a} 

if 

2.5 



4.3 



At 



4- 



-6.1 



V 
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Figure^-?. A data card with three values punched .on* it, 
, V together with related FfEAD and FORMAT statements. 

_ • ^' / ^ 

♦ .The letter F (for floating point)' t^ells u$ ,that\ the 'group of 1^. columns 
.is expected to^ contain "the numeral for a real number, ^_nd v):\^n converted to 
internal form, it should be represented in the floating point scheme. . • 

In Figure F2-T each numeral on the data card contains a decimal point. 
It-'s easy and quite natural to punch the decimfid. point' with each numeral, and 
yhe^ it is punched, the positioh"^^ the decimal point on the data car4 will, 
- deci4e its ''representation in the floating point scheme. There %s provision, 
how^y^er, to omit the punching of the decimal point and have the^ computer fig- 
ure out where it ought t^ be in each field; that is, locate the place in the ' 
numeral where the decimal point is intended to be. Use of this feature requires 
a^^littl'e ifiore care than is probably warrantedr'-^.^th^(|^eginri^ig. However, it is 
worth knojfing about for several reasons . It will help you, for one thing, to 
understand the meaning of the" 8 in the field code Fl^.8 . 

% 

18 



22 



Properly interpreted this code means: The numeral fount?4in the next I5 
^, columns of the card is to be converted for storage in the floating-point scheme 
If the numeral does not have a decimal point, actually punched on the cardj^ 
then a decimal point is located by counting 8^ columns to the left from the 
15th, ^r rightmost; column of the field.. The decimal point is intended "S^ go 
iimnediately to the left of this eighth column. You can see thx^ .illustrated * 
in the card picture ^found-in Figure F2-8. 



F0RMAT(Plp.8,..0 



I 2 5 4 5 6 7 6 9 K) II' \Z 13 14 15 
I 1^1 I I I 



I I I I 



"413 ^ 



0EC1MA.L ;JblMT 
BELONGft ^ERE. 



At 



\ * Figure F2-3. F-type dat& ^without "punching the decimal point* 

The numeral ^3, which is punched in Columns^ 
' . and 8; will be interpreted as ^4-^3. 

The numeral ^4-3 ' is punched in Columns 7**and 8. *We assume that no other 
punches are present in the first fifteen columns Because no decimal point is 
.shovfii on the data^^card, the numeral will be interpreted as h.^ whe;i it ^s 
finally stored in memory under "control" of a field code F15.8, ' ' 

' j ^ ' / 

. We. sei^that F- field codes, which arWsed to identify numerals represent-, 

ing rear numbers punched on data cards, are somewhat more complicated than 
I-field codes* An easy way tg^ view the fom of these codes is. suggested in 
Figure F2-9* 




and 



field width 



a r^al 
number 



/ 

field 
.width' 



decimal point 
locator (d 
digits from 
right end of 
the fields) 



Figure F2-9. Genera^ forms for I- and F-fj^eld codes 
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Strictly speaking a field baa be of any width, from 1 to 8p, 

, the width of the card, not just 1^ columns as we have been illustrating, 
gimilaply, the deeinaal point l,ocator d can vary i a value, and, not just 
equal 8. as illustrated thus far. p'or the time B^ng, however, in order 
not to have to keep so many new ideas in our head at* once it , will be simpler 
to imagine the, field' widths are fixed at ' 15 and a value of d fixed at^ &.* 

Another point about writing field codes--and tMs one 'is worth remember- 
ing--is that you don»t have to write a series of identical codes. You can 
g^oup them and use a repetition indicator to tell how many are intended in ^ 
sequence. Thus in Figure F2-6 we could write 

■\ ■ ■ . . 

22 FORMAT (3115) 
inste«Ld of . " " 

22 FORMAT (II5, 115, 115) ' . 

Similarly, in Figure F2-7 we could write- 

100 FOgMAT (3FI5.8) 

in place of the tore tedious ' " . . ^ 

100 "format Cfi5 . 8, F15 . 8, F15 8) 



Example i ' ^ 

( 

A card cAn contain any combination of I- . and F~ fields desired and in 
^ ' \ 

any.sequence as may be seen in Figure, F2-10. Notice that the one-to-one 

match is properly achieved employing format** number I5. *lhe first and fourth 
'l^list elements in the READ statement are real variables and the second and 
ithirS items are integer variables (according to the leading letters in their 

names). Correspondingly, the first and fourth field codes of format number 
jl5 are fields while the second and third are I- fields* This ^arantees 
|that assig^ied values for the rea^^ariables will be stored in the appropriate 
.floating-point scheme, while assigned values for the integer variables fire 

stoi:ed in the desired integer scheme. 



I 



I # 
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READ \%{ 
15 F0Rb^X(F\5.^'^^^ 




Figure F2-10., Two. field types on the same data card 



A 



Exercises F2-3 Set A ' v - • ^ 

- We Imagine a class of veiy simple problems tp be solved in the computer. 
Ut the flow chfert for each of these have a structure identical with t'Ke one 
in Figure F2-3, In the following exercises you are given Box 2/ Your^job 



in each case Is to:^' 

A. Decide what should be' in Box l7 
* B. write an approp:^iate READ stat^Wnt. 



\ 



C. Write a companion FORMAT statemenW^ 

jD. Draw a pidture of a typical data bard which could be read 
as a result^ of executing the READ statement which you have 
written. 



, -I ' " ' c 

Z<-((axx+t)x X +\c) X X + d 
' ■ f -I 



Q <- /(u - |) •: 



'X V 



5.,.; X 2/(Y +|) 
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AREA Xjir^ - (V3 X /r^ - + -r^ x PHl) 



Executing a PRINT statement ^ . ^ , 

Printing urfder contror of a given format is analogous to ^reading, but 
in a reversed sense. The value that is currently assigned .to each element of 
the output list will appear typed or printed across the page by the output 
devic^-. When each nunlber' is (Copied from its place in memory iti is conve^^ed 
for printi^ng to the external form that is dictated by the matching field ^ .-^-^ 
code. The spacing of the numerals on the printed line, is set automatically - ' 
by the* widths of the corresponding field codes. As you 'can see, the fie^' 
codes, in a sens^ control or govern the appearance Of the tirinted line. 
Several examplesWill help you tb see this process. 

^Example 1 * ^ . \ 

Examine Figure F2-11 where you will see an exampl^ of a PRINT statement-- 
and what it might accomplish when executed. ITotice the format is essentially 
the same as that used in an earlier example for a R^AD statement. This 
illustrates the 'fact that the same, fgrmat can te used in connection with any 
READ, -PRINT, or group oi* such statements, so long as efich refers to th? 
'•format by the same number. - , ^ * 



For most simple problems, it is helpful to obtain a computer-produced 
printed copy of •iVery piece of data that is read in to the- computer. This 
pripcess is sometimes referred to as -^gho checking the data . Of ^ourse it 
isn't essential to print the data values in precisely the same order as they 
may have appeared on the data card. Thus, a program may begin with state- 
ments like * ' ' ^ ^ . ' , 

1 ' READ 22, NUM, KPAY, KAMT 

• . * ■ V 

2 PRINT.22, NUM, KAMT, JCPAY 

^ \ 22 FORMAT (3I15> 

Data read as a -result of executing statement 1 will be printed as a 
result of executing statement 2. The echoed data can then be checked 
visually to see if the computer got the right information.' Both statements 
refer to the same format. 

Now since statement 22 has a passive or reference role, ratjier than 
an active or "executable" role, its position in t^e sequence of statements 
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which forms a progr&m isn't really critical. . To emphasize this distinction, 
we shall henceforth call non-executable statements like the.fprmat, ' 
declarations, because, in a sense, they provide "declarative", or descriptlv 
' information only, in this we can reserve the word statement for the - 

class of action or. executable steps of a computer program. FORTRAN has 
,^i6.veral o|fer declarative statement tjypes. You will be introduced to these 
^in later^ chapters. 'You might" now check Figure ^-h and notice how the 
FOR>IAT: declarations in that program are placed at- the very end (just ahead 
of END) t" They could -^just as well have been placed ahead of statement number 
1. Many people prefer to keep FORMAT declarations near the input or output 
statement which refers to it (or near the first of these). 



PRXOT 22, NyM, KAMT, KPAY 
22 FORftAT (3113^ 
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* assumes : 

current value of NUM* is I502 

current value of KANPT is 17I+ 

current value of KPAY is II+572 
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• Figure F2-11. Example 'of a PRINT statement, its referenced 
-.^ ' ' FO^^'r^ and ar possible line of printed Results 

^ Claused by the execution of the PRINT statement 

• r. ^ ' ' 

Ma^Cimum line^vidths ^ 
♦ • * • 

You may have wondered-^frow many numerals can be printed on .one line, ^and, 
if more than on« line is needed^ how the format can written for use in 
controlling such a process, Unlike ijie punched card which has a fixed number 
of columns, the maximum width of a printed line is not so well standardized. 
For a given class of printing equipment there is some degree of uniformity, 
however. Thus characters or 13^; characters are common "sfkndards" . 

^For our text^ we will assume a maximum line width of 120 columns, or 
characters. Hence, if'we use fixed field widths of 15 columns, we' will be 
able to prin^; up^ to eight numerals. p^r liri'e. . * 
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E»ample 2 _ 



Now suppose we want' to print, say ten values, as a result of executing 
a sipgle-i^HiT-statement. Since many printers' use type spacing of 10 
c ha;ractg TS,.tp^he \oj:i26nta^l inch, the 120 character width line is 12 
inches--across. For many purpo^s this is wider than we actually need. To^ 
Ijeep our page width narrow enough to fit in our notebook or In an 8^ X 11 
report, let's suppose we agree to print* not more than four numerals on any 
one line. Figure F2-12 illustrates one way this anight be done. We use as 
our format (UFI5.8) v^ich is interpneted to mean: 'U^ to four items of the 
form Fl^,8 ( per line ) 4 The numerals on each printed line are to acbount 
for to four lis-C elements of the PRIOT statement; Wh^n ^the PRINSJ, state- 
ment is executed, values are copied out of memory one at a time from positions 
associated with A, B, C,^ etc. Each copied value is then converted for 
printing in the desired^ extepial real form, as illustrated. 

\ 



PRIOT 31, A, B, C, D, E; F, G, H, P,. Q 
31 FORMAT (to5.8) 




Shows Appearance of 
printed value for and G 



s * 



Figure F2-12. Showing printing o^ up to four items per liile^ 
each item under oontrx^l of an FI5.8 'code ^ 
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If a line printer device Is used, then, when four such numbers have been 
converted, a complete l«lne is printed -at once. However, if a typewriter 
device Is used^ each number is typed as soon as It has T^een converted to , 
output form. This difference In behavior ^s no vital consequence as the net 
effect Is the same. In any casj^, when four items have "been dispensed with, 
tbe :5fermat, {hFl^^Q), is reused and Jhe next four num*bers are matched with the 
four F15.?*Sj^ converted and printed. The format is again reused. This time 
the last two list elements are matched with the first two of the F15.8*s, 
converted and printed. Since the last of .the ten values has •now been copied 
out, and the number properly converted under control of^its matching field 
code, the list is discovered to be exhausted , i.e., no more ^tems rertiain to 
"^pe copied. This discovery then signals the end of th^* process. Execution 
of the PRIIJT statement is terminated and the codputer is free to execute 
whatever statement happens to be next in the program. 



Exercises Set B 

For the following output lists, write PRIIIT statements, each with a 

* * * • 

suitable FOBMAT declaration to go along with it. Use only 115 and tl^.Q 
fie^d codes. Assume a maximum line width of 90 columns, i.e., six items. 

1. /f, B, J, K, and L * 

2. A, B, C, J, K, L, Al,, Bl, CI * 

3. IKE, JAK^ BAKER,' CHARLY, D, B 

Exercises Fg-B Set b . - 1 / 

1-6. .In these six exercises, you will continue with\tlte d(Svelopment begun ' 



earlier, in Section F2-3>, Set A.^^Jour job now in ea^h 'exercise' is to: 

a. Decide what sfiould go in Box 3 of the flow chart. 

tatemlht iwt 
intent of Box 3. 



... _ ^ . ,^x-^' 

b. Write a PRINT statem^t jwhich, if executed, would carry out the 



c . Sh'ow an ac c om;^ny ing FOF MAT^ de c la rat ion*' 

d. For the data you used inlthe first set of exercises, compute the 

inted and show its expected formatted 

on the printed pa'gfe. 



result which would be prj(: 
value as ♦it would appear 
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2-9 



J2-U 



P2-lj Agsifinment statements ' ' , 

- We shall first examine the- parallel existing between our previously 
developed flow chart concepts of assignment steps'and those of RJRmil 
statements. 



4 



FloV Chart 



specific^ 
example 



general 
^ form 



variable 



arithmetic 
expression 



\ 



variable 



arithmetic 
expression 



The TORTRM variable is a character string built up of letterarand digits 
as described in Section F2-2, * ' ' 

-The^ arithmetic expression in FORTRAN corresponds to any meaningflil 
computational rule which uniquely defines a single numerical value. There 
will be a few important restrictions to observe in writing such expressions 
correctly. Before considering these, let us see how several flow chart 
examples are rewritten in FORTRAN, 



Examples 



Flow Chart 



X <-2.5 



.2: 
3. 
\. 

to 

5. 
6. 



aj X b 



Q <- Am - |) X n 



PORtRAN 
X = 2.5 

Z = 2.5 + T 

.T = ABSF(A) ^ b/c 

Q = ^QRTF((m-4.5)^EN) 



AREA 



AREA = 3.14159/2. * R**2 - (S*SQRTF(l^2 - S**2) + R**2 * HJl) 
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We have already become acquainted with three of the basic components of 
arithmetic ^expressions; namely, cbns1:^nts, variables, and ope^ors* 

^ Two otKer important components are parentheses .and function referencefc ». 

In FORTRAN we shall use parenthesfes in several different ways. If you 
look at Example 5, you will see the 'parentheses* used to form subexpressions 
that force a desired 'ordering to the computation. This is the use which 'is 
' perhaps most familiar to us. , 

- / ^ ' . . 

In Example 3 you see the parentheses used in a new way--to enclose the 

arguin|ht of a function, A reference to a function value, qt function refer- 
ence, for short, consists of the name of the i\inction followed .by a pair of 
parentheses enclosing tlie argument. This convention, . in FORTRAN, which is 
used for all flections, is often used in mathematical notation, 'For example, 

•„ ■ " • ^ - 

may mean "the function f evaluated with argumei>t x?'; or simply, "f of- x^ 

Since some functions are given special marks or symbols, like | |, 
^ , which are simply not available in the FORTRAN characte^ set, it makes 
good sense to' use the parentheses. We then get a uniform way ^f or writing . 
function references. 

< 

How are the parentheses used in Example There are two pairs used 
'here, one nested within another. 'You caa probably see easily that the outer- 
most parentheses serve to enclose the argument of the SQ3RTF fu|iction. This 
argument is^ the value which' will .be obtained upon evaluation of the expressior 

^ ^ _ '(EM - lf.5) * EN , - , 

Here the parentheses are employed for orderJLn^ coprputation. ' 

As you can see, a function argument may, itself be a FORTRAN expression 
of decided complexity. A similar example is given in Example 6. 

Number types in arithmetic expressions ♦ " 

In ge^al, <^ expression in FORTRAN must be homopseneous . this we 
mean the variables, constants and function references which are joined by 
arithmetic operator's in one expression, all must be of the same type— either 
all real or all integer. (There is one epcception which is treated a bit 
later,) , ' ' - 
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5^or example: ^ A + 5* 

6, * T 



*.sqrtf(x) 

67absf<t - P) . ' t . ^ 

are examples of homogeneous real "expressions. * 

i ' K + 4 

. J/k ^ 3 . _ ' 

^^^^^ • ^ . • ' 

are examples of homogeneous integer expressions. ^ :< 

In the case/Of real expressions, /each arithmetic operation is carried 
^ ' ^ out in real arithmetic. The value- resulting in any of these expression 
*^ evaluations is developed internally in floating point representation^ 

][n the 'case of integer expressions, each arithmetic operation is 
conducted using integer arithmetic and the resulting value is\ developed' 
internally in integer representation. / 

Target code produced from an arithmetic expression is generally more 
efficient when the expression is homogeneous. Moreover, the results of each 
operation are easier to' define. For these reasons, and other purely pragmatic 
ones, the^ original FOI^TRAN languages and compilers were desigped and con^ 
structed to reject, as illegal, any expression*«that vas not homogeneous . "/ ^ 

Many of the later versions of the FORTRAN language retained. this 
restriction. As a result, expressions like , , • -* 

'I + 1. ' ^ 

A + k 



7/t + 5. * U 



\ 



■ (M - * N '■ 

are illegal FORTRAKI You can see why, in Example k, names like EM, 

' 4 - - . * . 

in place of M,' and EN, in place of N, were used in constructing a 
hpn(ogeneous real expression as the desired argument for SQRTF. 
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■ _ F2-l| 

Because in customary mathematj.cal notation di not normally impose this 
restriction against '.'mixing modes" of the operands, you may ^ind yourself - 
very frequently writing illegal FORTRAN^ without- meaning to. 'jf you fail to 
make the necessaiy changes to correct the errg^, fear not! • "when your .progr^^n 
is read by the computer, under control of the processor program, it will 
certainly l»e rejected fir this reason. In many cases, In faqtlthe processor 
^ will find most or all such illegal expressions the first tim4 tL program is--' ' 
examined. Ifot only will you receive- a 'printed rejection slip, bVit al Auch ■ • ' 
violations which h^ve been found will ^ identified rather cleart So, it 
^ will be vejcy easy to correct this type of errorJ \ - * ' 

; ' - ■ ■ A - ' 

Special case of exTPonentiation vith integral ppvera - \. 

You may have pondered if indeed all homogeneous expressions are\ illegal 
. ^ in FORTRAN. Are there no exceptions? Some of you may recall th^ di\cussion 
in the main text on' exponentiation. In order to distinguish between 

3 

^^-^ ' , a with the meaning a x a x a 

and 

a"^ with the meaning e-^^"^ 

^ we agreed that two different ways to write a^ were needed. In 
^ notation i 

A 3 

* I 

is permitted so that ve can imply the computation a X a x a and distinguish 
it from - < 

> 

^hich.is an order to carry out the computation e^^^. Note^at a more 
cumbersome. expression equivalent to' A 3,, can be^^^^men in'^FSKTRAH. /Jt 
is: -EXPFCs. ^ LOFGCa)) which also means ^^^^^ . 

\n short, raising a real expression (like to an integral power (like 3) 
* is. the one and only form of nonhomogeneous expression which 'is permitted.! 



'In more advanced versions of FORTRAN like F0R1?RAK IV, there are other forms 
of nonhomogeneity which are permitted. . ■*' 1 
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F2-; 



Which of the two expressions, 

' ' *A**3 or AJH*^. 
I " 

^do you think is easier to carry' out manually? By computer? The obvious 
answer is A ^ 3 in each case, . • ^ 

To compute e^^'^ requires the determination of a logarithm, and the 
raising of e to a power. These operatipns are carried out in a computjer 
-with Ibhe use, of the LOGF and the EXPF functions, which are separate programs. 
Though automatically supplied when expressions like A,** 3« are used, each 
ordinarily involves from 10 ' to 100"' times as much computer time as a simple 
multiplication . We can conclude therefpre^that, \^en a choice is available, 
expressions like A 3 are to be preferred to A 3*, because a more 
efficient target program will result. 



Integer division and its relationship within the greatest integer functi6n 

Integer division in FORTRAN plays an important role in algorithmic pro- 
cesses because it is related to' the TRUNK function which we defined in Section 
2-5 and which is in turn related to the logically powerful greatest integer or 
"bracket" ([ ']) function. Spe<?tfically, the FORTRAN expression l/j is 
equivalent to the^thematical expression .TRUNK(|) for;all integers land- J.^ 



1. 
2. 

3. 

k: 

6.* 
J. 



FORTRAN 

Expression Computed Valtie 



§710 

'-X0/{-10) 
^ 11/10^ 
10/1 
-5/10 
-15/10 
10/(-l) 

; l/(-10) 



/ 



0 . 

1 
1 

10 
' 0 

-1 

-10 
0 



MathemSi^cal Equivalent 

[9/10] 
[10/10] 

[11/ro] 

[10/1] 

-[5/10] 

-[15/10]' 

-[lo/i] 

-[i/io] 



Notice that when I 
the FORTRAN expression 



and J ^uj^e both, positive or are both .negative, 



I/J, 



is equivalent to the mathematical Expression, ^ [^] . On the other hand, if 
the sign of either I or J, but^not of both, ds negative, th6 FORTRAN 
expression 
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Spec 



is eiuivalent to th$ ta^thematical expression, -[|-|]. 



si "Q"^e on the exponentiation operation ) 



You Inay have Pondered how^we would write expressions like 



F2-if 



r 



in FORTRAU. Ei-ther we mean 



E 

(1) A (B E)/i.e., A^^ ) 



(BXE) 



or we mean 

' (2) (A B) E, ks/l (A^) which is really A 

As you caa/^ee, they'll only%e the same when B^ is the same as' B X E. 

We must therefore conclude that the^ operator is not associative, 
and in FORTRAN it is invalid to write 

^ A B E 

becaiise it is considered aijjjDiguous. One of the two unambi^ous forms (l)-9r 
(2), whiqjiever is, desired, must be explicitly written. 



Functions which have function values as their -argumentg** 

In some of the expressions you viH see in the next exercis?, the argu- 
ment of a function is expressed in terms of the value of 'another function. 
This is perfectly pennissible in FORTRAN, as there" is no restriction on the 
complexity of an arithmetic expression when used as an argument of a functions*' 
Ano,ther ex^ple of this^^ich we saw earlier was the ejcpression 



6XPF(3. * LOGF(A)) 



^-^'gument of EXPF ^ f *• 
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Exercises F2-^ t?Set A 

3/2 

1. We vould like^o express 



in FORTRAN, where A > 0. 'Keep in mind 
that, A-^'- is, in this case, the same as vA-^, 
the jTollowing corre^ly express K-^^ in FORTR^. 



.3/2 



or (a3)^'/^. All of 
Some are awkward. 

Some have supeiffluous operations. Coraiaent on each and choccse the one 
which appears to' be the most efficient computationally. * 

ABSF(A 1.5) 
(A **'3) 0.5 



..a . 

b'. 

c. 
d. 
e. 
f . 



I 



•SQi^TF(A *^ 3) 
' ABSF(SQRTF(A 3)) 

(A 1.5) 
^ABSF(A) ^* 1.5 

sqrtf'(absf(a 3)) 



2. 



If A can have negative values, which of the ^even FORTRAU' expr 

/p 

given in the preceding exercise correctly expresses |a| ? 
than one, which is simpler computationally? Explain. 





ssioris 
' more 



- Unary minus ! 

In Secti6n of the flow chart text you learned to distifguish ^between 
the*unary a;id number- naming minus, on the^one hand, an^ the binary minus, on 
rthe'other hand. If a m^nus sign_ appears at the very beginning of an expressi 
or immediately following a left -parenthesis, it is either a^ unaiy minus or h 
. number--naming minus. It cannot ,be a binary minus. ^ ^ ... 

Here are- ^ome examples in FORTRAfJ statements. 

^ i. Q = -5. ' ' ' ' ^ 

' • 2. ^ Q = -(A + B) , 

3. / Q*^ -A ("C) * . ' 

U. * Q = i-h.) 2 • * ^ • 

Q = -(if. *^ 2) 



5. 

6. 
7.' 



Q ^^k. (2) . 

Q = f.k. ^* .2 » 

Q = SINF(-A + B * (-COSF(C))) 



Jlote tY]^t {k) assigns a value of +I6 to Remembering from Table 

^ of the flow chart text that exponentiation takes precedence over unary minus, 



.ve see thaf ''(5')/t^)~and Xl) each .assign -i6 to Qi. Similarly, -A (- C) 
^ . is the sama as '{£^^^(^9 ' * 

Of 'course,, nq tvD ' operators "may^e^ ;^^^ Thus, A X -B 

is invalid. We^ aa5±ea_4 A (-£) " ot perhap's -A ^ B . Similarly, - 
/V*5][4_C or- rA^^» - C ^ are boUi nivalid. * ^^^^^_w 



Exercise's F2-^ Set B ^ ' ^ * - 

Correct the following three- invalid FORTRAN statements.. What pro"6lems 
arise, if any, it^- correcting tj^e second an^ third Statements? 

1. T = T * -A ^ \ A 
\ 2. F = C/.t3 + 

3. G = A + B ^iC -F/D) " . . c 
^ . . * » 




F2-5 2^ order of computation in a FORTRAN expression ' • ' 

^ ^ "a ' . ' ' 

We have not deferred thi^/questfion to this. late point because tfiere is 
something special that must be saijd here about FORTRAN that is different from 
what we havjs already said in Chapter 2 of the flow chart text. On the 
contrary, the fules for interpreting the ordjsr of computation are precisely 
the s§me. If you don't recall these, clo^e the .book now and try to recon- 
struct tTiem, Then compare them with the following. 

1. When- parentheses a^e used to nest one subexpression inside another, 
evaluate these nested subexpressions in the order from the innermost 
to the outermost. , » . ' 

Within any one subexpression evaluate: in, descending order of precedence: 

Highest function references 

* 

(exponentiation J 

/ 

^west +^ - 

3. In case of a tie in precedence level (within any subexpression), 
perform those operat,ions in the tie from left to right. 



F2-^ Meaning of assignmerit vhen the variable on the left Is of^ differenf 
type from the expression on. the right 

Is it possible tp convert a number from integer to real - representation 
or vice versa? Based on our discussion so far^ it would «seem not. To answer 
this question properly, you should notice that .until now all the assignment 
statement's we have illustrated were homogeneous in the sense that the variable 
on the left of the ( = ) sign and the expression on the right of it were both 
of the sam^type (real or integer). There are two other, obvious possibilities 
in FORTRAN. They are not only bo,th legal, but highly u^efiilrk 

*( i 
In short; we have four cases: ; ^ 

(a) real variable = real expression i 
'or . J 

• (b) integer variab].e . = integer expression ^ 

( c ) real variable = integer expression ^ / 

(d) integer variat)le = real expression . ^ - 

We have nothing rfurther to say about cases (a) and (b). It is (c) 
and (d)^ we are interested in because, such statements can be used to convert ^ 
integers to reals and vice versa. ^ ^ 

Let's first^consider case (c). The number assigned to the real 
variable simply has. no fractional' part, < 

Example _ ^ 



T / 1'+ J " ^ 

■ , / , / 

Observe that this sequence of FORTRAN statement^ will lead to a real ^lue 
for T equal to 22.0 . In other words^ the integer sum of I results 
in the real vilue^for T having a zero fractional partv-"* 

Now consider case '(d). Here the integer vjlue 
is tnihcated 



led to'^the iva^iable 
le'^SBTiSE-dBsCTibed in Chapt^jv^^^ection 2-5/ In other 
words, when the expression is real and t!{e variable is integer , then! i' 



variable = expression 



really means 
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variable -= sign of expression X 



\ 



expression 

absolute 
"value 1, 
_ marks 

greatest 1 
integer ^ 



] 



What this t?oils down' to is that a FORTRAN statement liHe 

' . ^ I = A 

is the equivalent of the »f low chart assignment: 



I TRUNK(A) 



t which is by no means the same as 




A. 



Why? Bec^iT^ ^the-ilow; _chart variabiles I and A do not have specific - 

them/ hen^t^"" 1io-roundlng ( lopping-of f 



digital representations aaaociatec 

-y . * — 

can be iiHplied in the simple flow chart assignment. 

* ^ * 



Example 1 



r 



DIAM 5.9 

ICIRC = 3 •1^159 * DIAM 



This sequences of FORTRAN statements leads to a value of I8 for ICIRC, 
and^no't l8;535 which would be the ^actual circumference of a circle with a 
diameter of 5*9» ' ' 1 ^ * - 



Example 2 



BALNCE =52*51 . 
WITHDR = 92*i+9 
IOWE^= BAWCE - WITHDR 



Assuming BALNCE and WITHDR refer to bank balance and withdrawal , .the 
overdraw is' i39»98» The fractional part, •98, is loat when storing ,-3$, 
in lOWE, ' ' . . 



n-6 

Exercises F2-6 

■ ' > 

1. Which of the following statements -are invalid in FORTRM? Explain. 

(a) Tl = Tl * VA|?3A * ' ' ' % 

(bX^A = EXPF(A4 * Z + A3 * Z 2.) 

' .(c) Y = LOGF(SINF(f)) + 22 ■ . • 

* ** » * 

(d^ J = J + 1. , . 

r^fe) I = I/K ■ • V V . " • 

V 

2. 'Assuming 1(b) above is valid, what changes would you propose in ^he 
/ interest oi* efficiency? 

In each of the following wxite a sequence of one or*more FORTRAK assign- 
ment .statements to. accomplish the indicated task* 7 

3. Assign to IPASS the integral part of the real value now assigned to V. 



4. Assign to FPART the fractional part of the r^al value , currently assigned 

• , to V. ' ' * ■ 

* , 1 , • ^ 

5. Provide ah alternative real representation for the integer vaiue currently 
assigned to the variable called IKTV. .Call this new representation RINTV. 

6. Compare ^ the 4*ollowing FORTRAN statement wiyi this accbmpanying flow chart 
box ' • 

- . ^ J = F/K 



Are they the same? If not, xihange tlfe flow chart box to conform with 
the FORTRAN statement. ' ■ . " 



F2-7 Writing complete programs 

Remember the ^0 bill problem? You will now be able to write FORTRAN 

« 

statements, indeed a complete FORTRAN progi^am, whicVi will show hqw to compute 
the most I^O-bilis obtainable from the given ( real ) PRICE of Dad's Jersesr 
cow. Before looking at the version given below in~ Figure F2-13, write your 
own first and then compare the two. 



Label 



FORTRAN PROGRAM 
Statement or declaration 



C WHAT THE JERSEY COW WILL BRING 

READ 101, PRICE 
C PRICE IS GIVEN IN 

C DOLLARS AND CENTS ^ , 

PRINT 101, PRICE , 
101 'FORMAT (F15.8) 
C NOW CONVERT "fO PENNIES 

IPRICE = PRICE ^100. ^ * " 
C COMPUTE NUMBER OF 20»S /. 

NUM20 =: IPRICE/2000 



102 



15 



NOW PRINT ANSWER, I.E., 
VALUE OF 'NUM20, AND STOP 
PRINT^IC^, ipi20 • 
FORMAT (115) 
STOP 

END ' - ' . 




0 


T's.Tsnonnooj 




0 


0 




it- 


0 


0 






0 


0 






0 




1 ^^^^^^^^^^^^ 




0 








0 



Picture of a Data Card 



Picture of Results 



Figure F2-13. The Jersey Cow - programmed in FORTRAN 
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Ex^rcices F2-7 

1-6. In the exercises of Section E2-3; Set A and C, you worked out the 
input and output details needed fo^^i^c FORTBM programs each having 
the simple loop structure shovn in Figure F2~3. You're to finish the 
job^now by writing out on a coding sheet each of these six simple 

FORTRAN' programs . ' ^ ^ - - - 

' X . / S 

7. What single assignment statement can replace the two that are used in 
the pro'gram given in Figure F2-13? ^ 

8. Recall the problem {in Section 2-5 of the main text) to simulate winning 
points on a carnival roulette wheel. It's presumed you have already, 
drawn a fl9w chart for this situation. Now write a complete J'ORTRAN ' 

t program which is equivalent to your flow chart. 



/ 
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F2-8 Some clerical 'details 



This section introduces you to some' additional details concerning the 
5)reparati9n of punched cards for FORTRAN programs* 



Length of a statement ^ • 

A FORTRAN statement (or declaration) can be, for all practical j)urposes, 
as long as necessary. If it is necessaiy to continue any one statement on 
a series of lines, provision is provided in the coding form to indicate 
continuation . This is the purpose -of Column 6 on the coding' form shown in 
Figure F2-4. Thus, suppose we choose to write the statement 

T = SQRTF(A ^ 2 + B 2 +'C 2) 

on two line^s instead of one.^ Its appearance on a coding form might then be 
modified as shown in ^'igure F2-lt, 



Label 



Statement 



V 



T = SQRTF(A ^ 2 + B 
2 + C 2) , 

Column 6 



Figure F2-ll|. -^Showing use of the continuation code 
on a coding form 

i 

The digit 1 is used here on %he continuation line . Any digit is 
satisfactory except zero. Notice we leave Column 6 blank on the first line. 
The cards punched from these instructions are pictured in Figure F2-15. 



COL, 6 



/ i/[2 -t-c^^iy 



, Figure F2-15. Showing use of the contin\iation code on 
the punched' card 



F2- 



Blank spaces (columns ) in a statement 

A blank 'column in any statement or declaration punched on a card is, 
with one exception, alvays ignored in FORTRAN . To illustrate this point we 
can let a 6 represent a blank card columft. Then the statement 

* A(3DEt3IUD=QBCJAoKDEdRD^llX^ 
will be treated in the same way as 

ABEL=BAKER+10.1^ 

So'blanks can appear anywhere ^nd will be ignored. Moreover, the A in 
ABEL need not be punched in Column 7. A statement 'may be punched beginning 
in any column which one chooses, between 7 and 72 inclusive. The same 
principles 8pply to READ, PRINT and 'FORMAT. - The exception, where blank 
columns do have significance occurs in H-field codes used in FORMAT 
declarations. This subject will be treated in .the next chapter. 



i 
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The printer carriage *^ ^ 

f The carriage of a line printer/ which holds and rolls the paper, must 
receive instructions as to how im^ch to roll, or feed the papervforward prior 
to the printing of the next line- Note this is analogous to an^lectric 
typewriter receiving a carriage return impulse which causes the paper to move 
forward one un:v"t or line in the vertical direction. Also note tha't vertical 
paper movement on either a typewriter or on a lin^ printer is forward , only ? 
If we are going to space the printed material ^o as to get various vertical 
arrangements, we will want to have the computer issue spacing commands to 'the 
>f)rinter carriage telling how many lines are to be moved, prefe^abl^ before 
each printing action. 

How is this control achieved in the FORTRAN language? To see how this is 

9 

done, we^should first visualize a"string of characters which has been developed 
for printing as one line* Call this a line image . These chara'cters are, -in 
general, the numbers converted to their output appearance governed. ^?y one or 
more field codes. The printer which receives this line image is so*wired that 
it can (and does) "shunt off" tlie leading character (i.e.^ the left-most in 
the string) and interpret it in a sp*ecial way, Instead of being treated as 
the first of a series of characters to be printed, this character is received 
as a coded signal whicfi then activates carriage moyement. Depending on the 
character, the carriage moves various distances (i/e., feeds the paper forward 
various amounts). For this reason we ^call the leading character of a line* 
image the parriage control code. It is never printed . 

The three codes which all printers are 'wired tp ui^derstand^ and th^ir 
.resulting effect are: 



Carriage 
control code 



Effect 



space . 



□ (blank) Rolls paper forward 1^ 

0 (zero) , Rolls paper forward 2 spaces.' 

1 (one) ^ Rolls paper forward so as to* position it 

* at the top of the next page* ^eay^.ng a 
, one-inch margin at the top.' '5B'"page" oil 
a continuous printer form roughly corre- 
sponds to a single ticket 'on a roll of 
movie theatre tickets or a single "square" 
# * on'*a roll of paper towels. In oth^ wor.ds^ 

a page consists of the. paper between two 
horizontal perforations. In the cQse of 
the printer we shall assume these perfor- ' 
atlons occur every ele-^en inches. . 



'This section can be omitted if you communicate with your computer via a type 
writer rather than by punch card for input and by line pointer for output. 



The next question is, what technique can we use to force the character, 
hlank, zero, or one, as desired, into the leaditig position of each line 
imQge which is to be printed? ♦ 

For the present suppose we limit ourselves to single space printing . 
For this we will only need to force 6 blank into the leading position, 

You have probably notice^ that upon printing a number under control of 
an I- or F~field, if the niameral is* small and does not fill the field , J 
positions to the left of the leading digit print as blanks. This is our 
clue. A blank ;>ill therefore automatically be pr^ent at the left^ end of 
the first field of a line image if • , 

'1. this field is an I- for F-field, and if 

2.. the field width is irider than necessa^Nfor printing the 
desired number and its 'sign. 'r'ti ^ ' 

^ Now if you check the examples used so far in this text, you will find 
that with a field width of 15' columns, most of the simple numbers we have 
dealt with are in this iiategory. 

For a field coded as M 

F15.8 . . , . 

what is the largest number which can be printed and stil^ have a blank at 



Jbhe left-Tnost position? 



To answer this, interpret the X's in. Figure F2-l6 as digits. 
^ 15 

T 

.OTOOOQCX , 




Figure F2-l6. Printed form of a leading nymeral under 
FI5.8 control * 



Clearly, a number as large as ^ 

± 9999.99999999 

can be printed. 

For printing integers ther« is even less problem with a field coded II5. 
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P2-10 Input' and output of alphanumeric data , . - ^ 

Computers can store and manipulate alphanumeric data (letters, digits 
and special characters). This idea was introduced to you earlier in the 
flow chart te^ (Figure 2-1? ). You my have wondered how, i'n FORTRM, we 
can describe the input, output, and manipulation of such data. 

Since*FORTRAN was originally designed to describe algorithms for 
computation <Vith numerical data, it should not surprise you that its scope 
is limited when it comes to describing processes for handling alphanumeric ^ ^ 
information. 

Nonetheless, we can perform a few elementary alphanumeric processes 
via^some FORTRAN statements. These actions, though simple, will permit us 
to carry out some surprisingly complicated information processing. We» can^ 
easily describe in FORTRM how to input alphanumeric data' (to memory), as* 
we fehall see shortly. Assignment statements can be used to move alphanumeric 
data to new locations in memory. Thus, if a variable T, for instance, has 
acquired an alphanumeric value via input, then subsequent execution of the 
^statement,- 

• • S = T 

would assign to S • the alphanumeric value of T. Of course, we cahnot 
perform any meaningful -^arithmetic operations on alphanumeric data; but, as 
we*^ shall see in Chapter 5, it \?ill be possible, to compare two alphanumeric 
quantities for equality. Finally, ve c^n express in FORTRM how to print 
out* alphanumeric data from memory. » • 

Sto rage of alphanumeric data 

Up to now we imgined FOR^HRAN permittitlg^^^^ two classes" of variables; ^ 
rpal and integer numbers. But alphanumeric data is n^thej? real nor integer. 
H^w,then, can alphanumeric data be stored in memory? FORTRAN II processors 
permit you tp use/storage locations that are'assigned to numeric* variables 
for other purposes . 

We^shall in theofolloving'discussions imagfn^ a computer whose word length 
in bits, is such that, one word stores six whole characters, (if ypu have f9r- 
gotten how the Characters, may be coded:, each as a group of "binary bits, you 
should review Section l-k,) 
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^ample * • . 

Suppose the instructor who posed the original problem to compute 

' . ^ D = Va2 + b2 + c2 

has since been introduced to computer programming and. has some inkling as 
to th6 power of computers. He now poses the problem this way: 

"Imagine that several different sophomore geometry students have given 
you values of A, ' and C corresponding to the edges ;of a rectangular ' 
prism. You are to compute "for them the distance D, which represents the 
length of a diagonal, according to the formula suggested in Figure 2-1. 
Write a program which prints values for A, B, C, and the computed value 
for D, ■ and then also prints for identification purposes the full name of 
the student, hi^ room number, and seat location; like 

BOB JOHNSON, ROOM 3^2, ROW C, SEAT h." 

.1 ' 

If identification of this sort were, punched on a card, might look 
like that shown -in Figure F2-17-. ' ' . - 



/ 



aBOBajOHNSON,D R00MD342pH0WQC,o SEATcAoonQ 



1 



Figure F2-17^ Alphanumeric identification (l.D.) card 



With respect to the structWe of the algorithm 'first flow charted in 
Figure 2-1, little -has changed when we add, as we now must, the steps for 
reading and printing the alphanumeric identification, as shown in Figure P2-l8. 
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. 0 



ALPHANUMERIC 
IDENTIFI CATION 



A, B, C 



A, B, C, B 



ALPHANUMERIC - 
IDENTIFICATION 




! ^ 



Figure F2-l8.' Flow chart with provision ^ior "reading and 
^ ... ]^rinting alphanumeric identification . 



The question is, how do we write Boxes 0 and h in FORTRAN?* 

♦ 

To write Box 0 in FORTRAN we must first imagine the informat^Lon on 
the I»T>» card^^divided into six-column fields. Now, recall that ye are 
supposing each group of six characters can be st^jdfed in a single memoiy 
word, say, one that is normally used for storing real numbers. (Any other 
character packing capacity would not change what we are saying here in 
principle.)/ Then, any r«al variables^ we choose--you name it— such 



as B, S, T, .U, V, W and X or Rl, 
as elements of th^ input list. 



IC 



be suitable 



/ 



R2, R3, etc., 
So, in effect, 



will 



k6 



we want to think of alphanumeric information that has been punched in the 
^^card as being partitio^d into groups of six characters each, fot storing 
into locations associated .with some group of variables. This idea 
suggested in Figure F2-19. ^ 



15 



READ 15, Ri, R2, R3;.Ri+, ^R5; R6, 
FORMAT '(7A6) 



U )8 



P « t8 24 36 



42 



Rl 



RZ I 113 

I 



R4 



RS R6 I R7 



At 



Figure F2-19j. Illustrating use of AljJhabetic, or 
* " A- field codes 



The ffifernal representation of alphanumeric^ data follows a scheme, which jLs 
^ different :erom either the. integer or the floating point scheme. So, a 
^pecial conversion procedure is needed and hence a special field code is also 

'reeded to call for this conversion. 'The field code that is used has the form 



A .for 
alphanumeric 



width in columns 



During execution of a READ statement, if ^A^A-field code like a6, is 
encountered wheif consulting the referenced foraat, cthe charSbters matched. to 
this code, are converted, for s^orSgTiH the ^ .alphanumeric scheme. ^So, executing 
a READ statement when linked to a format, as one shown in Figure P2-19, will 
prove suitable for inputting alphanumeric data. r ^ 

^ Many other possibilities, some even more subtle, woiild also work. We 
shall consider one of these alternatives presently. Right now, we want to 
get^ an idea how we might write Bo5c the output box, in FORTRAN. We can 
take advanl^age of the natural symmetiy of input with out-pubT The statement* 



•Will be suitable. 



- PRINT 15, Rl, |2, R3, T<k, RJ, r6, R7, 
FORMAT (7A6) 



Rir 



hi 



'0 i. 



F2.10 



In executing this PRINT statement the . a6 fields encountered in the 
consulted format will force the reverse conversion of the alphanumeric data 
from the internal scheme, to the external scheme^'and printing will then be 

accomplished in groups of six 'characters per field. 

* ■-' ' ^ ^ 

The complete program requested by the instructor might appe&r as shown 
in Figure F2-20. 



Label 



Statement* or declaration 



C ' EVALUATION OF D 

C EACH SET OF DATA^ A, ^B, AND C IS PUNCHED 

C ON ONE CARD,^ BUT IS PRECEDED BY A CARD 

C ^ CONTAINING ALPHANUMERIC IDENTIFICATION 
1 ■ READ 15, Rl, R2, R3, Kk, R5, R6', R? 
15 ' FORMAT (7A6) 

READ 101, A, B, C 
101 FORMAT ikFl^.Q) 

.J) = SQRTF(A**2 + F^*2 + C**2) 
PRINT 101, A, B, C, D 
PRINT 15, Rl, R2, R3, R^^, R5, R6, R? 
GO to' 1 
END 

Figure F2-20. Program 'showing- facility for alphanumeric 
input and^utput . 



In examining Figure F2-19, you may havfe wondered about the form of the 
information on the I.D. card. Except for the fkct that the information 
will be grouped -into "words** of six columns each for storage, the name, 
room number and seat loca^on is in free f o5rm « >That is, name, room 
nmliber, etc.^ don't have to be punched in certain fields on the card. More 
over^, with the excej)tion of a blank in Column 1^ whic^ will be explained in 
a .moment, no special attention ia given to the number of blanif colupis be- 
tween words. 

^ ' / ' 
How then would' our program, (Figure F2-2o) handle the I..D. carcJ for 
Algernon Thistlewhaite whose I.D. card is shown In* Figure F2-21? ^ 
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30 



36 



42 



46 



HI 



R2 fa> j R4 

I * I 
I I 



R5 
■Ar- 



R6 R.7 



Figure F2-21. I.'d. card Vith a longer name"^ 



The last six^ characters on this card will not be read by the compute r^^ecauge 
provision was made for storing information from only the first forty-two 
columns of each I.D. card (7A6). If we t/arit*to be safe, we had better make 
provision to store all 8o columns of the card. This means 13 full words 
of six characters each. plus a parti^ly filled word containing the last two 
characters from the^ card (Columns 79 a;id 8o). In other words, we need an 
I input list of Ik items with a governing folmat code like* 



To revise our program so it \;?ill both read ^nd print the full I.D. ca^d^ 
we can* replace the affected READ, PRINT, and FORMAT statements with those ^ 
shown here: * / 



4 



READ 15, HI, R2, R3, Kh, R5, r6, R7, R8, R9, 
RIO, Rll, R12, R13, Rli^^ 
15 FORMAT fl3A6,A2)^ 



FRlm 15, Rl, R2, R3, R^f R5, R6, R7,,-.R8, . 
R9, RIO,' Rll, R12.,, R13, Rl\ 



Now to clear up the question of Column 1 on the I.D.. card, and why we 
suggested i? should be blank. The xdystery is easily solVed if we recall our 
disicussion of printer carriage control from Section F2-9-. If ,our output comes 
from a line printer we need to ensure that the first character of an output 
line is shunted-off and used as a carriage control code instead of being 
printed. By insisting on a blank in 'Coluran 1 of each I.D. card, we guarantee 
single spacing of each printed line of I.D. information. If you us^ a type- 
writer for output--no harm— the first position of the typed line will be blank.^ 
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taiapter F3 
BRANCHING AND SUBSCRIPTED VARIABLES^ 



F3-1 Oonditional statements 



/ 



-Ixi,:gection 3-1 of. your flov chart text you studied techniques for 
branching by means of a two-exi't condition box^_ ijn Section 3-3 you will 
study multi-exit con*ttion,|^s. Bra^iching instructions may be written in. 
"FORTRAN ^y means of an, IF statement. The FORTRAN IF statement can be used 
as either a two-exit or three-exit conditional statement. In practice, how- 
ever, the IF statement is most frequently used for two-way branching and less 



frequently' for three-way 'branching, 
statement. 



Figure F3-1 



>ws both uses of the IF* 



FORTRAN, 
Statement 

English 



IF(I - 5) 2, 3, h 



\ If I-5<Q go to Statement 2 
?* . IfUJ"5=0 go to Statement '3 
i<Xf I-5X) go to Statement k 



One .way of 
drawing the 
condition 
box 



Another way 
of drawing 
the condi- 
tion box 




as a three-^way 
conditional 



^*IF(I - 5) 2, 2,^k . ^ 

If I-5<0 go to Statement 2 
If I-5=:0 go to Staitement 2 
If 'l-5X) go to Statement 




' as a two-way ■ 
conditional 



EMC 



Figure F3-I. ^ Examples of FORTRAN IF-statement$ 



Da 



An IF statement consists of the vord 

. ' • IF ■ . • 

followed"by a FOBTRM arithmetilc expression enclosed in parentheses,^ for* 
example, ' ' • 

. • . ' (I - 5) ' \ 



followed by thiree statement numbers, for example. 

The CQniplete Statement is 

' IF(I - 5) 2, 3; ^ - 

It tells the computer to evaluate the expression 

I - 5 

and choose among the following alternative^: . , 

if I - 5 < 0 go to Statement Nimiber 2; 
^ ^ if I - 5' = 0 go to Statement Number 3; 

if I - 5 > 0 go to Statement Number 



Thus the w statement is a test of whether the value of a given expression 
is negati-ve, zero, or positive. .The three statement numbers specify where 
'to go .In each casev . * 

The secoh4^ example in Figure F3-1 shows how the three-exit IF statement 
can be used as a tyo-yay branch. There is really no requirement that, the 
thTee specified statement numbers in the IF statement be all different. Two 
' of them can be the same ,-^yand when two are the same, we get a t^Ho-way^ branch. 

Thus, we can test whether the value of I - 5 is non- positive (l - 5 < Q) 
' by wri"feiflg. 

. IF(I - 5) 2, 2, U 

Jt is important to observe that this gives the- ^ame effect as if we ask 
whether I < 5 is true or false. 

" - . ' ' . \y 
Exalfnples ^ ' Consider the statement IF(X - ^O) (>)\1) 7 ' 

■ l*i What will be the ne:ct ^statement executed if ji = i|.7? (Answer: 
Statement IBumber' 6, since h,7 - 5»0 = -0.3 is negative,) 

^2. . What vill'be* the ne:ct statement executed if X'^'^ 5.^? (Answerl 
^ Statement Number since' 5.0 = +0.4' isi positive* ) ) 
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3. 



What will be the next statement executed if X = 5.0? (Answer: 
Statement Number 1, since 5.0 - 5.0 equals zero.) 

Based on the above IF statement, complete the following: If X 5.0, 

the next, statement eScecuted will be Statement Number 6; but if 

X ^5^0, the next statement executed will be Statement Number 7. 

(Answer: If X < 5.0, the next statement executed will be Statement 
Number 6; -but if X > 5.0, ' the next stai;ejflent executed will be 
Statement Number 7.) 

Draw a two-exit condition box with exit 'arrows labelled "t" and "F" 
corresponding to the IF statement given aljov^. Give two alternate 
solutions. Answer: 

, t- 

(x-F.o > 0^ * 
6- 7k - 





Which is easier to understand at a glance? Answer: The author thinks 
it will be the one on the ri^ght in most instances. 



6. 



Refer to ;^he lower right hand condition box ^in Figure F3-1. Redraw the. 
box without changing* its intent that the T exit now leads to Box k 
while the F exit leads to Box 2. 'Answer: 




How has the corresponSlng IF statement changed? Answer: It hasn't. 




/ 
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F3-1 ^ ^ 

Exercises F3-*l Set A '"'^ - 

In Exercises 1-5 use the follcfwing IF Statement: 

IF(K - 8) 20, 30; 20 

1. What will be the next . statement -executed if K = ^i? 

2l What will 'be the next statement executed if K = ^? 

3. What Vill be the next statement executed if .K = 9? " 

,h. Complete the following: The next statement executed will be 

.'statement Number 30 if and only if K 8. v> 

5. Complete the following: The next statement executed will be 
Statement Number 20 if and .only ^ if K 8. 

6. Write an IF statement that wj.ll go to Statement Number 9 if X is 
posi-Plve and to statement Ni^mlJer 10 if X is non-positive. 



In Exercises 7-12, transform each of the given equations or inequal- 
ities int*o an equiyalent form having **0" as one ^member.: ^ — ^ 



7. K >7 (Answer: t - 7 > Q) ^ ^ ' . ^ 

8. X > 8.1; "ii ^ ' ^ 

% Y<ii.2 \ . . \' ^ " 

10. A < B . 

^ 11. X =-Y ^ ^- ^, ^ . p 

♦ 12. X + 5 >5 Y . ' T . • • 

• ' • . » ^ * ' , <^ 

Converting -^o a three-way algeb-ijaic test ^om any one of six 2-way: • 

relaiions. ' W ^ * * • ' 

Figure F3-2 shows how one might go about writing IF statements corre- 
spom^-ng to t^ra-exit condition boxes involving one of the six relationallW k 
symbdls - 

W ' .' J '. 

In every case illustrated we convert a relation whose form is .A relational 
symbol B, to a difference whose form is ,A \ B, , wher^^A ^ and B . ai*e the 

•e3^essix)n8 in the conation 'box tWt lie op either side of the relational 

••symbol. '\ . • ' * 



Flow Chart 




'What is True? 



Reasoning 



'iC-7 < 0 



Is K-7 > 0 T or F n -^^F 



Go to Statement # I lO" 
T ' I 



K-7 = 0 



10 



Fortran 



K-7 > 0 

T IF(K-7)10,10,20 



20 




What '"is True? 


Xw8.if<o 


X-8.i^= 0 


X-8A> 0 




Is X-b.i^^ T or F 


F 


T 


^ T 


IF(X-8.U)11,11,21 


Go to Statement # 


11 


11 * 


21 


What is True? 1 




Y->.2 = 0, 


ft 

y-i)-.2>o 


• 

J 


Is Y-k.2<0 T or F 


T 


: F 




^IF(Y-U,2)22,12,'l2 


Go to Statement # 1 


22 


12 


12 



What is True? 


lA-B < 6 . 


A-B = 0 


'A-B >J9 


is A-B<0 T t>r F 


T 




' F 


Go to Statement #1 


1 23 


23 ? 


13 ' 



IF(A-B)23, 23,13 



^What is True? | 


X-Y < 0 " 


*X-Y = 0 


'^Y > 0 


Is X-y=0 T or, F 


F 


T 


' F 


Go to' Statement # | 


ih ^ 


2ii • 






IF(X-*f)lU,2lt,U' 



What Is True? 



* ^ f y^ Go-tQ Statement #1 \^ 



X+5-Y=0 



15- 



'J ■IF(X45-Y>25,15,25,. 



35 




25 
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Figure 13-2. Six' varl^tlong-on^he-IF-^hemfe^rf gr -Uo br anches 
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Table F3^1 shows the six joattprns possible in us,ing the IE statement 
as «a two-way branch. / 



' • Table F3-1 

• t 

The six patterns for statement numbers in two-way IF^s 
going vto Statement 1 if the/ relation is true, 
^ going to Statement 2 if the relation is false 



Form of Relation 


IB 


Statement 


A > B - 


. IF(A - 


B) 2, 2, 1 


A > B, ^ 


IP{A - 


B) 2, 1, 1 


• A < B" 


ifCa - 


B) 1, 2, 2' 


^ ' A </b 


IF(A - 


B) 1, 1, 2 


/ A y4 B ^ - 


IF(A - 


B) 2, 1, 2 


/ aU b 


if(a - 


B) 1, 2, 1 



Among the importajit things to remembe:i^about' the form of the IF 
statement are the following:'' 



(a) 

(b). 
(c) 

(d)- 



The. arithmetic expression to be tested is enclosed in 
parentheses folloWjirig "IF". 

The three exits are rtegat^^ zero , and positive Jn that 
order. . ^ ^ » . j - , 

A statement number must be vritten for each of the three 
exits, in thp. orde.r indicated in (b),,.even though on^ of 
the* conditions may never occur. * ^ 

Th^re must be commas between the statement numbers but 
none between the right parenthesis and the first state- 
ment- number. , ^ ' 
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P3-1' 



* ■ You may recall the Ruritanian Postal Regulati6ns problem. Compare* the 
executable steps of the following^ FORTRAN program with the flow chart in 
Figure' 3-1 and verify that the two are^-^qtilvalent. 





• 


Label 


Statement or declaration 






1 




READ 100, N, A, B, C 


* 








D = SQRTF(A**2 + B^^*^ + 0**2) 


/ 


■ o 






IF(d - 29?0) h, h, 1 ; - 






h 




PRINT 101, N 










GO TO' 1 






IDO 




FORMAT (16/ 3^6.2) 






101 




FORMAT (16) 










END 



Exercises . F3rl 'Set B 

Write a FORTRAN IF statement for each of the following condition boxes: 



1. 



o 5 



W > 0 




20 




(x < 9*1) 



2. 



i 




15 




■ 4; (k - I. £ -q) , 

30 y 



6'. * (a B > CJ 



Identifying - remarks in FORTRAN output ^ 

In Section 3-1 you worked with output boxes which included identifying 
remarks as elements- of ^ the output list. For example:^ 



.ERIC 



"A is the larger. The 
value of A is", A 



1)0 



Note that an identifying remark is a constant, more precisely an alpha- 
numberic constant, as defined in Section 2-6 of the Main Text^ There is nq 
way that the remark itself can be changed in value as a result of carrying 
out the steps of the flow chart. 

In FORTOAN an identifying \remark such as 

"A is the larger. The value of A is" 

does nqt appear directly 'as an item in the output list. Constants are not 
legal elejnents of an output list. Instead, ve insert the remark in the FORMAT 
declaration. Thus far, you have encountered the following types of field 
specif icati'ons in FORMAT declarations: I, F, and A, To these we now add 
the H-f ield, 

^ life H-field (l!d3J.eritli fl,eld) may contain the text of any remark, message 
or h^^ding,"^ i.e,, any dlplian^umeric constant. Instead of enclosing the text in 
quotation mJbks, we precede it ^y the letter H (identifying the type of field) 
and precede that letter by 'uStsigned integer indicating the length (number* 
of charactei^) ofj^he teaj that immediately follows the letter H, Figure F3-3' 
gives seveiral illustrations, ' ' ' , 
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Flow Chart 



FORTRAN and the 
resxiltlDg printed line 



J, "IS 
THE 

LARGER" 



"THE 

LARGER IS 
J=", J 



X; "IS 

GREATER- 
THAN", Y 



"THE . 
LARfGER IS 
J=", J 



PRIKT 100; J 

f 

100 ^ FORMAT ( l6j l4Hc3ISaTHEoLARGER') 



XmXX IS THE LARGER 



—4 

PRINT 102, J 
102 FORMAT (l9Ht7rHEoLARGERQlSOjat=0, l6) 



I THE LARGER IS J = XXXX5CX 



PRINT 101, X, Y 
101 FORMAT ("F10.3, ITHQlSoGREATERatniANo, F10.3) . 



XXmX.XXX IS GREATER THAN XXXXXX.XXX 



PRINT 103^ J 
103 FORMAT (l9m:THEdLARGERcCCSoJo=o, 3^ 



THE LARGER IS J = XXXXXX 



Figure F3-3. Use of ^H- fields for Identifying remarks 

- ^ Note i;he 'dlf ferefice hetveen , FORMAT statements-^Og^ap^^ 103^ , In, 102 the, 
first character of output for the line Is a blank, while In 103 the first 
character Is the numeral "1". Tf you recall th^ dlsciisslon of printer 
carriage control in Chapter 2, you wiH. remember that the blank initial char- 
acter calls for a single space before printing, while the "1" 'calls' for a 
skip to the top 6f tlfe next page before printing. If you begin a line with 
an H- field, and if you do your printing on a line printer Instead of a type- 
vriter,' you must remember to put in this first carriage control -character* 

for each line. If the first fleld^would not nomally be an H-field and you 

/• 

want to be certain that the leading character will signal the carriage control 

' ' i 

of your chgice, you can use a one- character H-field at the beginning of a 
FORMAT statement. , . ' * • 
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Carriage Control Desired i 

singly space, 
double space 

skip to top of next page 



Initi-al H-ffeld , 

—7 ■ ' 



■■1 



,0 

IHl 



For example, we might rewrite the first format in Figure F3-3 as 



100 



FORMAT (iHo, l6, llfHoISoTHEolARGER) 



in order to gueirantee a.blalik' as. the leading character for use in carriage 
control. ' , 



Exercises F3'-l Set C 



/ 

V 



Write FORTRAN PRINT and FORMAT statements for each 9f the following* 
exercises. Assume that flow chSrt variables beginning any letter I 

through N are integer variables (and use l6 fields). Assiame that all 
other variables are floating point variables (and use- P 10. 3 fields). ^ 

Insert a one -character H-field at the beginning of each output record to 
provide carriage control on the assumption y.ou are using a line printer for 
output, 'In EXeacises 1 and ^2, aingle space before printing. In 3 ""double 
space. In \ skip to the 'next page before printing. The choice of Format 
numbers is ^p^to you. - \, ^ ' < 



1. 




2.* 



- - < 



X,,'^'IS^ 
VALUE OF X"' 



\ 



3. 



\ 
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Exagiple 1 



The second form of the flow char^ in* Figure 3-5(b) could be"progratnmed 
in PORIRAN as follows: •. 



Label 


statement or Declaration 




READ 100, A, B, C 




ALRGST = A 




IF(ALRGST - B) ^1, 5, 5 


k 


ALRGST = B ' • * 




IF( ALRGST - C) 6, 7, 7 ^ 


6 


ALRGST = C 


7 


PRINT 101, ALRGST 




STOP' ' " 


100 


FORMAT (3F10.3) 


101 


FORMAT (22km{raLARGE^TOVALUEDISa, F10.3) 




END 



Example 2 



The flow chart ±r\ Figure 3-6 could be programmed in FORTRAN as follows: 



-L bel Statement or Dec^laration 



k 
7 

100 
101 
' 102 
103 



READ 100, a; B, C 
IF(A - B) 4, 3r3 
IF(A 1 5, 5 

PRINT 101 A- - 
STOP . \ 
PRINT 102, C 

GO TO 8 ^ . . ' ' ; 

IF(B - C) 6,v7. 7 
PRINT 1D3, B^ 

GO TO 8 ^ * \ ' 

FORMAf (3F10.3) . " ■ ' -^^-^ 

FORMAT (32HlAOlSdTHEdLARGEST,DlTaaVALUEOISa, FIO.3) 
TORMAT (32}milSDTHEnLAR6EST.Dl!I^ FIO.3) 

FORMAT (ssffiBoisaGram^ FIO.S") 

END . i ' 



Example 3 \ ^ ^ ' . * 

The tallying problem in Figure 3-7 could be programmed in FORTRAN as foiUowB: 

Label Statement or Declaration » 
"^IREAD 100, N 
KCOUUT = 0 
LOW = 0 
MID = 0 

KHIGH - 0 ^ 
READ 101, T 
IF(T - 50.^)) 9, 9, 5 
LOW = LOW + 1 
GO TO 7 
5[' IF(T - 8o!o) 10, 10, 5 
'10 MID = MID 1 
GO TO 7 

6 KHIGH = KHIGH \ 1 

7 KCOUlPr = KCOUNT + 1 
IF(KCOUOT - N) 3,11, -11 

11 PRINT 102, iflCOUNT, LOW, MID, KHIGH 
STOP 

100 FORMAT (l6) 

101 FORMAT (f8.3) 

102 K)RMAT (43HpVALUE£dOFCfiO:^OUNT,OI/^^^ 
AREq, hl6) 

Exercises; F^ Set D - ^ ; * 

1-5. Ref^r to' Exercises 7- 11, Section 3-1 Set A. in your flow chart J:,ext,. , 
For each of these five exercises you prepared a flow chart of a simple *^ 
algorithm- Now write a FORTRAN program corresponding to each chart. Choose 
statement numbers , where needed, to correspond tb the box numbers used in 
the flow charts* Specifically, let 1 be ^ the label for a statement that 
corresponds to Box 1, 2 for Box 2, etc. * 

- • . \ 

Exercises F3-1 Set E * ^ ' ' 

1-6. Refer to Exercises 1- 6, Section 3-1 Set B. For each of these^six 
exercises you prepared a flow chart of a simple summing algorithm. Now 
write a FORTRAN program corresponding to each flow^chart. Choose statement 
nxaribers, Vhere needed, to correspond t6 the box numbers used in the flow * 
charts. Assume a suitable input format code for €he data Is F10«$ . 





^ • - 1 ■ 1 ■ ' ■ 

' . \ F3-2 

1 ' 






F3-2 Auxiliary variables * ' ^ ! 

^ i ^ '! 


The use auxiliary' variables in FORTRAN programs mirrors what you liave 
already learned in the flow chart text. • ' i 


1 

Exercises F3-2 Set A ' \ 


1-5. Refer to Problems 2 through 6, Section 3-2 Set A in vour flnv nh«v+ 
text. Write FORTRAN programs for each of the flow charts you have con- 
structed for these exercises. Choose statement numbers which correspond 
to the box mujibers used in the flow charts. ' 


f 

You might like to see hov the flow chart for the Euclidean Algorithm 
given in Figure 3-l4 m(ght work out in FORTRAN: 


* ^ • JuaDei 


Statem^t or Declarat'ion 


• 


'READ 100, KB 
IF(KA-KB) 5, 5, ^ * . 


k 


KB*= KA ^ * ' 




M = KB ' . 




KB = KB- ^ ^ \ ' 




IF(KA) 6, 7, 6 . ^ . ^ \ . 




KE = *KB - KB/KA-^KA ' ' * ' 
' ' ' 




KB = KA - " ' * ;r - 
KA,= KR ^ ; j 




QO TO 5 ' ' f ^ ^ 


' f 


PRINT 102, Kb' ^ < . ' ^ ^ ' ^ - : 




STOP - . ' 


. 100 


FQRMAT(2I6) 


101 

"102 


F0mT(l2HmfflbcK:ro0H^^ I6, 5H3AI©D,l6,lfHDIS.) 

format(ihd,i6)' . , : 














a' • \ 
ERJC . , \^ 
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Exercise F3-2 Set B , ^ - , 

Write B FORTRAN progrM corresponding to the*aow chart that you constructed 
in Problem 2 of Exercise 3-2 Set* B of your flow chart text. 



acercise F3-2 Set C , ^ 

1-8. Refer to Problems 1 through 8, Section 3-2 Set C ir5>$kir flow, chart 
text. For each of these eight exercises you have prepared a flow chart 
of a simple algorithm related to coordinates of points on a straight <> 
line- Now write a FORTRAN program corresponding to each flow chart," 

Choose statement numbers which correspoiid to the box numbers u-ed in the 

I 

floH charts. ' ^ 



V 



ERJC 
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F3-3 Compoimd conditions and mu5^j.ple branching 
* Writing compound conditions using IF statements 

^ In Section 3-3 you encountered condition boxes involving more than one 
decision, e^g., ' \ 



2 < X m> X < 




,and you saw that this single decTsion box was equivalent to a pair-of boxes: 



X < 5 





This pair of boxes can be coded in FORJRAN as 



IF(2.0 - X) 3 

IF(X - 5.0) 2, -3, 3 . - ■ ^1 

Similarly, the following pair of equivalent flow ch£^ts ' 

1 - ' 1 ^ _ I y\ 

(x < o) -T-^ » (x > ly 



f 

(X < 0 or />^^ 



-2 



- T 




i 3 



can be toded' as, 



IF(X) 2, k, h.' - 
lF(X - 1) 3, V, 2' 



"V, 



■ ) 




ERIC 



08 » 



Exercise*^ FS-J- Set 'A • - . ; - 

1 Bef^er to' Exercises 1 - Section 3-3 in your flow char^ 'text. Eor 

eAch of these exercises write j^fORTRAN sta"tonents equivalent to the /low 
^ ' charts you pre]paredX. For example, -the flow ch(art example would probably 



be. coded in FOR^^ as:* 
* • ' , Labiel 



IF(X1 - X2) 2, 30, 30. 
IF(? -Uj) 3> 3/ 20 ^ 
IF(V- S)" 30/' 20, 30 



Writing mult iple branching ins'truct-ions using IF sxat^ents 



In Section F3-1 you ledrne* that the FORTRAN IF statement was a kind of 

V • **> * 

three-exit (i.e.., th:;ree-branbh) conditional statement. That is 
IF(a£itbde€ic expression ) 10, 26, 30 



is equivalent to -the flow chart 




In^ certain situations you will l^d it useful to use the IF staitement 
this way — as 'a three-way branch. , Frequently, however, . when you need to do' a 
three-way branch, you will find that your problem does not^ neatly fit the 
pattern of a single IF statement. Usually you will fiiid it necessary ta ^ 
write two or more IF, statements c^orresponding to a i^ltiple branch conditi'on 
box. * , ' ^ • ■ , 
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Example; Flow Chart 



X+Y<5 




X+Y>100 



r 



FORTRAN: 



kl 



S(X + Y - 5.0) 21, ko, ko 
IF(X + Y - 70.0) 22, 22, kl' 
IP(X + Y - 100.0) 23, 2k, 2k 



I 



> Exercises F3-3 Set B 



1^ V7rite FORTRAN^ IF statements tov 



I>20 




r 



Write IF stater^nts>^foit 



3. Write a FORTRAN program corresponding to the flow <ihart you pr'epared for 
Problem 10, Exercises 3-3« ■ ^ ■ 

* Write a FORTRAll^rogram cori^espondin^o the new flow chart which you 



^ 'prepared* for the caxnival wheel proWLem as the answer to Problem 11, 
^Exercises ^^^^ 



^F3-^ Logical ejq)ressions 

' It would *be nice if it^were possible in FORTRAN to code the two-way 

condition boxes of qui* flow charts in a^^Ilore straightforward fashion, !ro*do 

this one would need to add relational operators (and symbols for thenl) to the 

set of allowed operations. For exanrple^ it w(^^d be nicQ if a condition box 
like^ / • 



'3 



could be coded as the single statement: 

) 

IF(A >■£>){, -f 

< 

or possibly in some oth^r^way^ like 




• meanl|;ig : go to 2 if true ^ 
-meaning: go tu 3 if false 



IF(A > B) GO TO 2 - 



- 1 



where here the W statement shows what to do if the condition inside the 
parentlj^es is true^ implying that ve should otherwise proc*eed to the next 
statement if false. ^ ■ - 

You may have surmised or have learned from other soiirces that FORtTRAN is 
more a set of similar dialects than a single^ language. The dialect FOBTRAN II 
which we are studying, has no relational operators and hence it has no "logical 
IF statement like tlje ones suggested above. 'It has onl^y the "arithmetic" IF 
statement. As a matter of fact, there ar^e FORTRAN dialects^ notably a group 
referred to as.JE;ORTRAN IV which do exhibit logical IF statements like those 
shown above. You can gue§s t^iat coding .'a flow chart into one of these 
, dialects of FORTRAN is an- easier job. " -r- . ♦ 



This, section can be skipped without loss o^" continuity. 



r 
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^J3^5 ' Subscripted variables 
. ' Representation of subscripted variables in FORTRAN 

Figure -F3-i^ shows how subscripted variables are represented in FORTRAN." 
Flow Chart Forjn. ^ FORTRAN Form ^ 

h . X(N) 

' . ^1+2 , . B(l+2) 

^°^0E-6 ^ • J0E(M0E-6)- 

' ' CHAR(5^I+i^) 

^^2^+15 . ^ / ALPHA(22*J+i5) 

^^lTxJAY-9 • , BETA(l7*JAY-9) 

Figure F3-4. Representation of subscripted variables in 
FORTRAN 

As you can see from the figure, ^bscripted variable^ aie represented in 
FORTRAN by enc'losing the subscript in parentheses and 'writing it following the 
. variable to .which the subscript is affixed. This is another exajnple of a 
^ notatioit'that enables FORTij^N code to be written "on the line. » Other examples 
you have seen include ';a**2" for "A^" and "SQR3F(X)" for >^,. Since 
iJ»lied ^tiplication (omission of the multiplication sign) is strictly 
forbidden, ^there is 'no danger than Y^^ will be confused with Y X N—the 
lat'ter is written using aii asterisk between the variables. The only sburce. 
of confusion in the notation is that a similar notation is used for function 
• references. You will see later in this section how this potential ambiguity 
is resolved. ^ ^ , 

• In your flow chart text, just atout any integer-valued arithmetic 
expression was pennitted as a subscript. In the ;yersion of K)RTRAN you are 

studying, the rules are much more strict: 
•i , 

(a) A subscrip-t may t^ike on onOy positive integer values (zero, 

' . negative numbers, ^and numbers other than integers i^ay .not occur 

as subscriplj values). ^ . 
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(b) The most general forms permitted are , 

* constantl * variable + constant2 

and 

constantl * variable - constant2 

vhere constantl and CQnstant2 are unsigned integer constants and 
variable is a non-subscripted integer variable. Other joermitted 
forms are , . - 

constant 
variable 

variable + constant *^ 
\ variable - constant 

constant * variable 

vhere constant is again an unsigned integer constant and variable 
is also of non- sub scripted integer type. 



Exercises Set A ^ / ^ 

Examine each of the f ollowing ^flov chart subscripts. ' If the subscript 
is legal in FORTRM, writ^ the subscripted variable ^n FORTRAN form. Other- 
wise point out the defect. 

3. CHAR. - . 8. AIPHA_j^ . 



5-. , ' 10. ABC 



>5xJ 



Allocation of memory storage for arrays 

As yoiiklalow, the jDomputer must have a storage location available corre- 
sponding to each varia'ble ±h a given FORTRAN program. If all variables were 
of .the form « . . * ^ . 



* ' : * . ' . * ,^ ' a;X; char, x^, etc. - 

this would be a simple problem indeed. The processor could merely assign a 



* 

ERJC ■. ■ _ 
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storage location to each variable occurring in yjour program. 

But what about X^? 'flow can the processor know for what values of I 
to aspgrt storage locat.idns? It^Ciinnot tell merely by looking at' the 
oCcurrenc"es of in your progr^.- But if it weie to wait until *the 

program was being ex^uted, it might find that it needed locations for 
X^^ X^, X^, X^^ and had assigned locations only for ^> X-j^q- 

Since storage of arrays s-uch as X-j^, X^, X^^ in consecutive locations 

is of overwhelming^ importance in* efficiency of program execution in most 
, computers, we would like to have advance knowledge of the range of values 
possible for a subscript before we start executing the program in which 
that subscript occurs. More precisely, we -require at least knowledge of 
the maximma range so that we vill allow enough locations. 

In FOI^TRM this problem is solved by means of a DIMENSION declaration. 
It is used to specify an upper limit for the subscript valine for a particular 
• subscripted variable. Figure P3-5 shows the form and use of the DIMENSION > 
declaration. .These declarations must precede all executable statements of 
^ program. 
f 

Now we cscfr explain how the computer can distinguish between such things 
as GHAR(J) meaning the J-th element^ of the array CHAR and CHAR(j) 
meaning a function reference consisting of the function name CHAR followed 
by the argument expression "j". if a variable CHAR appears in ' 
DIMENSIOll declaration, then a subsequent occurrence of CHAR(J) will be 
interpreted as a subscripted variable. If jCHAR does not appear in a 
DIMENSION declaration, ^then CHAR(J) will be interpreted as a reference to 
a f\inction named CHAR. (A common source of error in writing, FORTRAN 
programs is to omit a necessary DIMENSION declaration. With this omission 
the processor will^dutifully interpret all occurrences of CHAR( j) ,as 
function references. The result is usually a nasty bu-t misleading error 
message telling you about your misuse of the fi^nction CHAR.) 



V 



^ Declaration > 
DLMENSION X(5) 

DIMENSION CHAR(IOO) 

DIMENSION- X( 3), Y(k) 

DIMEi^SION 1(2) 



Resulting Storage Allocation 

five locations^: one each for X^, X^, X^, 
X^, and 

one hundi'ed locg^taons: one each* for 
CHAR^, CHARg, ^^loo ' 

one location each for 'X^v^, X^, 
Y^, Y^, and Y^^ 



two locations; one for I. --^and one for L 
t 1 . " < 



Figure F3-5^ Storage allocations for FORTRAN DIMENSION; 
declarations 



Input and output of arrays ^ 

The form *of the input and output lists for arrays in FORTRAN is quite 
similar to that in your flow chart language, ^vith a few small, but important 
differences in detail. Figure F3"6 shows how* arrays can be input and 
output in'FORTRAN* ' * 



Flow Chart Box 



FORTRAN Program - 
Label Statement or Declaration 



{KAY^, i= 3(2)N) 



100 



{X^j=l(l)l|} 



101 



READ 100^ N, (KAY(I), I = 3,N,2) 



F0BMAT(l2l6) .<^ 

PRINT 101, (X(J), J<= 1, -h) 
FORMAT (JHO, PtP12.6) 



•Figure F3-6, Input and output of arrays, in FORTRAN 



You should note these points concerning the form of oior notation in 
FORTRAN: '"^ 

{a) 'The increment for the subscript is placed after the upper botuid 

instead of between the lower and upper bounds, ^s practiced in the 
flow chart notation^ 



. ' Cb) When the increment is equal to 1 , it may be omitted. 
^ (c) Parentheses rather than braces^ are used; 

(d) The upper, -bound for the subicrii>t is set off by commas rather than 
by parentheses. ^ ' 

Example Figure F3-7 

^ Draw a flowchart knd write a FORTRArj program to find. the largest com- 
ponent of an n-component vector (n < lOCT) . 




2 




/ 


l(l)n) 



ALf^GST^A 
K <- 2 » 



( K > n , 



(alrgst > 



Ay^GST ^ A, 



K ^ K + l' 



8 



"LARGEST 

coMPoirar 

jIS^ ALRGSTj 



. FORTRAN 

r 

DIMENSION A(IOO)^ 
READ 102, N 

READ 100, (A(I), I =1,N) 

ALRGST = A(l) 
K = 2 

r * IF(K - N) 5, 5, 8 

5 I1f(ALRG6T - A(K)) 6, 7, 7^ 



6 ALRGST = A(K) 



7 K = K + 1 

GO TO 4 ^ 

8 PRINT 101, ALRGST - ^ " 

" STOP 

100 format(Ufi6.8) 

101> F0RMAT(^2H LARGEST COMPONENT fs, FI6.8) 
102 F0RMAT(I3) 
►END 



Figure F3-7. Find the largest number problem 
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Exercises Set B 

In Exercises 1-3 vrite '.the necessary DIMENSION statements' and READ 
statements to input the arrays indicated. 

1. {A^, i = l(l)k} Assume k < 50 ' ^ 

2. -CB , j = 5(2)n) Assume n < 125 

J 

3. {A., i = l(lO)n), (B., i = 10(2)n} Assume n < 50 



Exercises F3-5 , Set C ^ 

1. Write a FORTRAN program corresponding to the flow chart in Figure 3-2U(b) 
(The carnival vheel pi4blem using subscripts). 

2. Write a FORTRAN program that corresponds to Figure 3-25. Choose your 
format codes using 15 and F10.5 fields as needed. _ Assume values 
for B are punched on data cards vith up to seven values per ca rd. 

3. Write a FORTBEaN program that corresponds to the flow chart you drew for 
. Problem B of Section 3-5 Set A. Assume the value of n will never be 

greater than 50. You will have to write the program so that subscripts 
for the coefficients range- from 1 instead of from zero, because zero 
subscripts are not allowed in FOKiR^. . For format codes use 15 and 
F10.5. fields as needed for input. As§me polynomial coefficient values 
are punched in order on data cards having up to five values per card. 



Exercises F3-5 - Set D ' 

Write a FORTRAN program for the flow chart you dre^f<for Part C of the 
Exercise 3-5 Set. C in the Main Text. Assume by "anv^ize orchestra'^ we mean 
one that will never exQ^eed^. 500 players. Assume the ag§s x>f the players are 
punched as integer data (I5 fields, up to 10 per card). 
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F3-6 Double subscripts . • ' I 

' — 

^Representation of double subscripts in FORTRAN *' 

In Section :i?3-5 you learhed that a subscripted variable like 

r - ■ • 

could be written in FORTRAN as 

y' " ' • X(I) ,i ,, ■ 

Nov, iti Section 3-6, you have been introduced in your flow chart language * 
to doubly-subscripted variables such as 



The FORTRAN 'representation for such a doubly-subscripted variable is just 
about what you would expect it to be; 



X(I,J) 

The two subscripts are' separate.d by a comma. 
Thus the assignment. _box. 



WIN-^WIN + A 



could be written in FORTRAN as 

WIN = Wlif + A(K,J) . 

Allocation: of storage for doubly - subscripted arVays 

A DIMENSION declaration -is required to allocate space for doubly sub- 
scripted variables. For example ' 

DIMENSION A(3,h) / . 

This'-particular declaration allocates 12 spaces for a matrix having 
.3 rows and h columns. The subscripted variables of this matrix are : 



♦ 1)1 - 


h,2 


\3 


\^ 




^2,1 


^2,2 




\^ 




.'3,1 - 


- \2 


^,3 


A3 1^ 





; and fof the second subscript they are 1, 2, -3, and h. No other subscript 
values We significance with the above DIMENSION declaration. 

4^ 



ERIC 



P7S 
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To input or output an entire matrix, we can employ , either of two 
approaches 'suggested by the flow charts in Figure F3-9» 



{tAj^^=l(l)NiI=l{l)M) 




Method •(a) 



1^1 



/ 




= l(l)N) 




h \ 








I 1+1 






Method (b) 



Figure F3-9. Two ways to input values of a mat;rix having 

♦ 

M rows and N columns. 



The FORTRAN coding for each method is given below for comparison. In 
each case it is assumed that data are punched cm "cards four values per card 
using FIO.5 field codes. ^ ' e 



Method -(a) 



.Methord (b) 



ERLC 



READ 50, ((A(i;^)^J = 1, N), I = 1,M) 
50 FORMAT (too. 5) 



' . 1 = 1. 

. 3'?' READ 50, (A(I,J), J = 1,N) 
pO FORMAT (1;F10.5)' 
I.= I + i 

IF(I - M) 3, 3, -6 ; * 

: • 76 * 



In Method (b); Box 3 is an order to input one entire row of the matrix 
(the I rov). The other boxes define a range of values for I from, 1 
to M inclusive, under which Box 3 is repeated* In Method (a) Box 1 is an 
^•order to input the entire matrix rov-by-row* ' 

It would certainly seem that Method (a) mi^ht always be preferred 
because it is easier to write* In the next chapter we will show situations 
where the approach used in Jjlethod (a) h^ cer^in practical advantages* 

The FORTRAN program equiVjalent to the game flow-charted in Figure Z-Z^ 
is given below. The one-row-at-a-time approach for input of the 46x6 
matrix (Method (b)) is chosen* 



I 













1 — 1 






iDi?A'n inn ( h(T T^ t — i 
rsjtil) lUU, \i\\LfO)j 0 =5 1 






T — T J. T 
1=1+1 


• 




liJ^l - DV 1, 1; d 


0 




KituAi; lUl, jX; 1) 






WIJM — U*U « 






ATAQTi* — n n 
i\iJJoct — U * U 






± — ± 






u — ± 


\ 




mj -6) 5, 5, 6 


5 




WIN WIN + A(K,J) 






J = J + 1 






GO TO 


6 




IF(I . 6) 7, 7, 8 


' 7 




ALOSE = ALOSE + A(I,L) 






1=1+1 






GO TO 6 


8. 




AHET = win - ^ALOSE 






PRINT 102, ANEJT 






STOP 


100 




FORMAT (6f8*3) 


101 




FORMAT (212) 


102; 




FORMAT (150, F8*3) ' 


•0 




END 
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Exercises F3-6' ^ ^ « ^, 

1-5. Refer to Exercises 1-5, Section 3^^* For each of these exercises 

you h^ve dravn flow charts describing certairk."rov" or "column operations 
on a matrix/ For each of these flow ch&ts your job is to write the 
equivalent FORTRAN statements preceded by' the necessary DIM3SNSI0N 
declarations. No'FORMAT declarations will be required in your answers. 
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' ^ Chapter Fh - / 

LOOPING *. 

, * • » ft 

The DO statement ' ' 

It should come as no unpleasant surprise jbhat the lit^e wonder box we 

discovered and called the ^'iteration box*' has an almost perfect parallel in 

'FORTRAN. This shorthand is called the DO statemefnt. An example is shown in 
Figure Flf-1. 





Kr-l' 


I < 1000 


F 


I <- I + 1 







DO 15 I = 1,1000 



Figure F4-i, The iteration box and ^ equivalent; DO 
statement 



tlttedly, the p^aliel is not yet obvious, but we promise it will be 
show eac^in the larger context of a loop. 

To display a parallel between a flow chart and a FORTRAN "DO loop" we 
present Figure F4-2. The algorithm^ usedr^ you will recognize, is the Fibonacci 
sequence generator * (Figure 4-6). 

Our next discussion will focus on the loop counter tfiat.is used in a DO 
statement. It's "personality" will now be described. ^.Yoti may find this 
description somewhat long and detailed. Don't be discoiuraged. Read it once 
anjd then we will look at seme examples . 

1, It's- always a non-subscripted integer variable, like I, IKE, 
IBALL, J, M, .etc.^ and never lfli^'^(5), INT(K), etc. 

2. As in the iteration box, the loop counter must always be given 
■gf an initial value (> O). ' 
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Label Statement or declaration 




Figure" Fk'-^. Parallel between a flow chart and a FORTRAN " 
loop 



■4 ■ 



\ The po statement is ar^ order to repeat a certain task*,, once for each 
value of the counter I in the range 1 to 1000 inclusive. .The task 
which is to 'be repeated T^egins with the first statement that fpllo^ jbhe DO 
stQ;t^m^nt and terminates after executing the statement labeled 15, To mark 
the teminatipn^ of a "DOJLoopr we can use a Special "dimnpy", -statement or- 
marker for th*is purpose. ^It fLsveimply the word ' * - ' ' , 



HI 



ERLC 



• ■■ r 



^ONTINUE 



8o QO 



(Althougn we frequently use the COiJTINUE statement to mark the terminus, an 
assignment, RKAD or PRINT statem^t would also be suitable.)' We now see that 
the sl^itement 

DO. 15 I = 1,1000 , ' 

should be read this way: 

"Repeatedly execute (or do) all statements which follow this one, down 
to and including the (COmNUE) statement that is labeled I5, Do this once 
for each value ot I in the range 1 to iOOO," 

3.^ A highest value for the counter»s range mtisf always be supplied. 

h. An increment must be supplied if it is other than one. We shall 
.see examples presently. ^ 

In supplying Items 2, 3, and h in the DO statement,, we are privileged to 
give each either as an unsigned integer constant or as a non-subscripted 
integer variable. The rules are strict on tbhis point. If you forget and use 
some other Yorm, don't worry, the processor ^ill^ politely remind you by 
printing an error niessage for each poorly formed DO statement. ^ 

- 0 . ■ • I 

Example s ' * 

1. IX) 51 I = 1; 2 . , ^ 

means: repeat the execution of the immediately following statemejnts 
through the one labeled 51,. once for'esch value of the loop 
<^ountfer I,- in steps ' of 2 until I > If . Thus if Had ' 
^ a value of 7, the loop would be repeated \h times) for values 
•of- I of 1, 3, 5 and T. When the counter Is incremented 
again its value exceeds 7 and .we ; exit from the lo^op by 
. proceeding to execute the first statement that follows tne 
one labeled 51. f " 

2. DO 51 K-= 3, 8, 2 

means: repeat through 51 f once foi* each value of K until K > 8. K 
starts with a yalue of 3 and is ihcremented each time by ' 2. ' 
The loop is repeated, for K'-^^^ 3; 5 .aad 7= Ithree times). ' ' 



3. DO 51 J = I, L, X 



■J "AW 

loop y for / 'for \ N^for counter 

doimter initial' loot)' inc^ement^ - 

^alue limit r 

means: repeat Jihrough Statement 51, once for each value of the counter 
J until J > * J is given a starting value equal to that 
current]^ assigned to I. ' Prior to each proposed repetition 
- J is' incremented "by K. 

DO 51 ' J = 1^ L ^. 

} '\{ 

loop / for / / Vfor loop ^ . 

counter initial' limit . . • . 

value * • 

means: same as Example 3, except that no increment Is supplied. In' 
such an everrt we (and the computer that pust^ read this state- 
ment) always assume the increment S one. ^ • • 



Exercises F^-1 ' V . - ^ • 

Pro'blems 1-5: For eacii of the flow charts you drew in answering the exer- 
cises 1 through 5 Section k^-l, now Wijlte an Equivalent ^FORTRAN program. In 
order to write a complete FORTEIAN progrto, use the informatioji supplied "below 
fpr each of the e^xercises. " - • , 

■ L \' ' • • • 

4. Assume ID is 'a; 6-digit integer,^ and that values of A, and C 
can^"b'^ punche^d on cards- using F10.5 fields. ' * 



Use an 3;iO field-code for ID and F10.5 codes for A, B, and I)., 
f Assume N a number less than 10000. * • 

3. Choose a wide I-field like 115 for. printing values of I, LTERM - and 
' ^ S. (Is 5 ,a good'^name for our putposes in the FORTRAN program?) 

'Use 115 fields for all integers. Let the four value fi| of P^l^e^punched 

on a single card. I , . * , 

1 ^ • *^ ..." . ^ 

5. The Timekeeper records his data in hour^ ^d hundredths of an hour. 

Personnel keep data on rates in dollars and cents j)er hour. 



Illustrative examples 

There are a niimber of simple examples of loops in Section h'2 of your 
flow chart text which can be easily transliterated into FORTRM code with 
the aid of the DO loop. We shall use theses^'6* fu^^r illustrate some of 
, the details in the proper use of the DO statement.* Figure shows FORTRAN^ 

coding equivalent to the flow charts 'in Figure h.-Q. * 



Label 




Statement 






DO If I = 1, N 






Y = X(l)^^ 






PRINT 51^ X(l), Y 


k 




CONTINUE . • 



Label Statement 







DO 5 J = 1, N 






REAP 51, X 
■Y = X**3 










PRINT 52, X, Y 


5 




CONTINUE 



Figure F4-3. FORTRAN for flow charts in Figure k-Q. 
♦ • 

7 

Similarly we see in Figur^ m-k a FORTRAN equivalent of the flow chart 
in Figure U-9(a). • . • 



Label 


Statement 






SUM = 0. 






DO 10 J'= l^^^N 






SUM = A(J) +'^SUM . 


10 




COUlTOTJE 


k 




PRINT 50, SUM 



Figure 'Fk-k. FORTRAN for Figure k-9(a). 

/ 

Before looking tat our next transliteration^ please try your hand at 
writing the FORTRAN for Figure' lf-9(b) . Now compare your code with that 
found in Figure Fk--^. 
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Label Statement 



2 
'3 

h 
IQ 

5 



FMAX = ABSF(A(1)) . ^ ' 
DO 10-J = 2. N 

if(absf(a(j)) - mfx) 10, lo, k 

FMAX = ABSF(A(J)) 

CONTINUE 

PRINT 5"^/ IKAX ' 




Figure FU-5, FORTRAN for Figure U-9(b)' 

There are several lessons to "be learned here. If you look at the flov 
you see a line returning from the false side of Box 3 to the increment- 



ation gortion of Box 2, 

How^do we express this return line in FORTRAN? 



) 



You may have been tempted to!^rite code something like that In Figure 
Woe is you if you did!* , , 



Label Statement 



FMAX = ABSF(A(1)) 
DG 10 J =r 2. N \ \ 

tF(ABSF(A(j) - -FMAX) 2/2; V ' 
pIAX = ABSF(a(j)) \ « 

CONTINra 
PRINT 50, FMAX 




Figure lU-S. ^A mi^talSel 



A 



In that case you have falleriv^into. a great tig .FORTRAN trapj .Unfortunately; 
^when we order a return €0 St|atemen"ti 2 it is equivalent to returning to foy 2 
bf the flow chart at the initialization eompartmeni jnd not at thtf l increment - 
ation po!rtion ? This means we start th^.loop ^11 oy< 



instead of ' continuing with the c 



et again with J = 2 
nting ^nd testing process. In fact the 



.only ; way tb say in FORTRAN *go b^ck to the" inc^ementatio*);*^ portion of the DO 
stat^ement" is to send contrc^l forward ^to the statement* which ""is named as 
the terminus in the DO statement, ^ou ca^ see "how* why the 60NTINUE statement 
Is so handy to use as a terminus for a l50 \oop*. ""ESy giviilg Vt^e CONTINUE the 
label 10 in "this case and "funneling" all- flow to it, we guarantee a ret] 
;he inoreihentation portion of the IX) looj^. \ I 

With this idea in 'mind you should have no trouble writing ,phe%ORTfiAN 
fori the flow charts in Kigure i^^-lO and U-11. ^ . / 



ERIC 



81^ 



87- 



FU2 



Exercises Set A 

In each of the following four exercises, we present FORTRAN code for the 
flow charts in Figures k-l^ia), k-lQ(h), kAl{a) l^-Xl(b), ^respectively. 

Your job is to indicate what errors, if any, have been made in Ihe coding 
process. The necessaiy declarations are being disregarded here. ^ 



1. ^or Figure ^io(a) 



2. For Figure llrlO(b) 



3. For Figure i|-li(a) 



k. -For Figure ii-ll(b) 



Label Statement 







FMAX = ABSF(A(l)) 






INDEX = 1 


2 




DO J = 2, N 

IF(ABSF(A(J)) - IMAX) k, 3 






FMAX = ABSF('A(J)) 
INDEX = J 


k 




CONTINUE 


5 




PRINT 51, , INDEX, MAX 



Label Statement 



k 
10 



MAX = ABSF(A(2)) 
DO 10, J = ^ N, 2 
IF(ABSF(A(J)) - MAX)' 10, 10, k 
MAX = ABSF(A(J1) * 
CONTINUE 
PRINT 5, MAX 



Label Statement 







FA 


7 ' , 

:T =r 1 ' 






DO 


10 K = 1^ 1, N 






FACT = K * FACT 






PRINT 5, K, FACT * | 


la 




CONTINUE * * 1 



h4 •' 

ERIC 



' Label Statement 



10 



LTERM = 1 



NLT' = 
DO iC 



L, -N 



COPY ^ LTEBM 
LTERM = LTERM + NLT ♦ 
NLT =1 COPY ' . 
PRINT,' 5, K^' LJERM 
CONTINUE f . 




There is one more point that i? worthy of note about the DO statement. 
Many FORT^RAN II processors interpret the DO statement slightly differently 
than the iteration box would imply. * The iteration box implies that the 'test 
for termination is made immediately after the counter is initialized. If 
.the test happens €o ^ail this first time, then we never do execute the 
computation portion of the loop-^ot even once. Not so with DO statement in 
many FORTRAN II processors. After' initializing the counter, the test is 
bypassed and. we do enter the computation portion. It's only after completing 
an incrementation of the loop counter that the test i-s made. -Hence the 
commutation portion of a DO loop (for these processors) is always executed ^ 
at least once . Ask your teacher if the processor you are using in the • 
laboratory behaves this way. 

Supposing for the moment your FORTRAN processor is of the type Jus*t 
described, under what circumstances would the FORTRAN code in Figxire Fk-^ 
fall to agre^ with our flow chart in Figure il— 9(b)? 

The^ answer is-- if by some chance the value of N were 1, then, the 
printed value for MAX at Box 5 (Boxes 3 and k not executed) would be |a^| 
(unequivocably ). On the other hand, in the FORTRAN case, if [A^l ; were 
larger than |a^|, then the printed value for FMAX would be U^U since 
Statemfimts 3 and 'V are executed once. 

, JWell, then, what must we do to bring about strict equivalence between 
the flow chart and the FORTRAN^ under these circumstances? One way this can 
be done is to precede the DO statement with the necessary test to permit 
skipping over the DO loop. Thus we might write in this case 



10 
5 



FMAX = ABSF(A( 
IF(2 - 'n) 2 
DO 10 J = 2, N 



(A(l)) ( . 



,t^s is the extra- step 



CONTINUE" 
PRINT 50, FMAX 



Tq add this step is cumbersome and totally unnecessary in most instances. 
YoMT judgment is relied on here to decide, in these mattery. Obviously in 
^hls case, ,you are not goin^ to be looking for the l^ges*t pf 1 Qumbeys - 



with, a DO statement j 



'Fk-2 



Exercises Fk^2 Set B 



1 - 17. For each of the flow charts you drew in answering the exercises of 
Secti6n 4-2, Set A now write the equivalent FORTRM statements as ^ 
partial programs only. Do not bother to , write declarations unless it 
helps you to see what is going on. iint— the tricky where you 

should be .wide awake are: h, 5, 6, 7, 12, I3 and 17, 



You may be interested in seeing how the flow chart of the factors -of -N ' ' 
, algomthm ^Figure h-lh) is coded in FORTRAN. Here* it is in Figure Fk^-T. 

You Vill notice that all the data of this problem are integers, but 
the SQOTF function requires a real argument and produces a real result. So 
* the value of N is converted to a real value by 

' • FN = N -i-^ ' 

and FN i's then used as the argument of SQRTF in the next statement. The 
rest of ,the program should be easy to follow. You may recall that [n/k] 
can be written in FOI^^RAN as simply N/K, because "for positive values integer 
division produces ''.the same result as the bracket function. 



Latiel Statement or declaration 



1 
'50 




FACTORS OF 'N . 
READ 50, N 
FOFMAT (110)' . 

FN = N ' 

IBOUND = SQRTF(FN) 

PRINIT 51, N - ' ' 

F0RMAT(l6H3im3FACT0RSC30Fn,I15(^l^H3AR^ i 

DO 10 1, ifOIJND ' 

IF(N - K * (N/K)) 10, 6, 10 . ^ 

L = n/K ^ H 

PRINT- 52, K, ll 

formatCiho, 2115) 

continue . " , ^ 

STOP . ! 

EMD * . ^ 4 



Fig\are Fk-T . Algorithm for factors of N ' equivalent"' 
to Figure k-^lk. 



, . ^ ThaiFORTRAI/ co4ing Of the polynomial evaluation algorithm Figure 14-17 
presents an interesting problem. ' m the flow chart w^ ahow the polynomial 
coeffixxients^as' components of a ve/tor A beginning at f A^. In FORTRAN we j 
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r 



cannot have a zero subs-cript. The smallest value is 1. Consequently' in 
the Equivalent FORTRAlf program you will see how we have had to position the 
coefficients beginning .at ' A^. To do this a slight change in the algorithm 
is required, as you -will see vhen'you study Figure fh-Q. -r 



Label- Statement or .declaration 



^50 

10 
51 



POLYNOMIAL 

DIMENSION B{k) ^ 

READ 50, (B(I). I - 1, i^) 

FORMAT (i;F15.5) 

READ 50, X ^ / 

VALUE = B(l) 

DO 10 K = 2, 

VALUE = VALUE » X.+ B(kL 



CONTINUE 
PRINT 51, VALUE 

FORMAT( li^HbTHEnVALUEDlSa,F15 .5 ) 

STOP 

END • 



Figure Yk-Q. FORTRAN equivalent of the polynomiaL evaluation. 
■ (Figure h-ll) ^. 



Exercises Fi;-2 Set C' 



1-3. For each of the flow chart^ you drew in answering'the exercises of 
. Section h'2 Set B, now write the equivalent FORTRAN— ( full progranS^. 
Note special information supplied below: 

(l) Assume that N will never exceed 5O. Assume the values for 

X and A may be input \jsing Fl|),5 field coqles and that up_t2__ 
5 values- of X may b^ p.unched oji a single d^taf card. 
|2)* ^Make same ^sumptions Sbout the .data you raa^de iri^th^e preceding 

■ exercises. ' 1-^ ' , . . \ 

(3) If yoi/'continue to use integer variables, es^)ecially for your 

winnings, you may need to know how to^c^ll for t\e absolute value 
an int€iger expression. 'In FORTRAN II we use the ABgF function 
Tjhen the argument is^-real but we use t>ie XABSF function when the ' 
ai'guraent is integer. The resulting value bf*say XA:^^(k) is an 
' . , integer, ' . " 1 * . 



Table-look- U2 

In this section ve shall do two things: ^ 

1. Verify our ability to write the FORmN equivalent to the fairly " 
involved table-look-up (Figure 1^-24) using the- bisection methoi. 
Figure FU-9 shows the program. 

2. Learn, a little more about input -out put, statements and format code. 

The assumptions we have made in order to write the program given- in 
Fi^re F4-9'are: 

1. The table to be stored will never contain more than 200 X«s 
and 200 * Y»s. . r 

2. All input data can l5e formatted usinjfe I5 fields 'for the. integers 
^nd FIO.5 fields for the re^s.-* 

3. A data card read at Box 2 contains only one pfeir of valued for the 
X-Y table. 
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Label Statement or declaration 



2 
51 

52 

ifi 
11 

'53 



6 
<L2 
5^^ 



10 



1, N)- 



DIMENSION X(200),Y(200) 
REA9 50, N 
F0RMAT(I5) . 
READ 51, (X(K). Y(k), K 
FORMAT (2FIO.5) 
RKAD $2,A - . ' ^ 

FORMAT (FIO.5) ' • 
IF (X(1)-A)1H, hi, 11 
IF(A-X(N))5, 5, 11 
PRINT. 53, A <f 
F0RMAT(1H(,F10.5/ 

26HDiSalJ0ToINbRANGEd0FaTABLE; ) 
STOP , 
XOW =' 1 , , 
IHIGH = N I 
IF (IHIGH-L0W-1)7/12, 7 
PRINT $1^, XLOW, YLbW, A, XHIGH, YHIGH 
FORMATdHa 5FIO.5) 
STOP 

MID = (LOW + IHIGH)/2 ' 

IF (ArX(fflD))9, 9, 10 

IHIGH = MID ^ 

GO TO 6 

LOW = MID 

GO TO 6 I ' 

END 



Figure Fl^-9.^ FORTRAN eqiilvaleni of table-look^-up using 

bisection" (Figvire h-2h) ' / 
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i Implied IX) loop as list elements 

^In^ studying the program in Figure\F^-9; and in comparing it with its 
* corresponding flow chart you may have. not iced that in Statement 2 we used 
the input list element - ; ^ ^ 

(xCk), y(k), k i 1, n) . . . 

This 1^ the FORTRAN equivalent of our new flow*cha;:t notati/on 



Executing the statement: 



{Xj^, Yj^, K = l(l)N) 



READ 51(X(k), Y(K); K = 1, N) 
51 F0RMAT(2Fia.5) 



is a shorthand 'way of calling for execution of the statements; 

DO 10 K = 1, N ' , 0 ' * 

• READ 51, X(K). Y(K) . ' 

. ' . 51 F0RMAT(2Fi0.5} ' 

* 10 CONTINUE 

You can now see why we frequently refer to a notation ^llke 

- , - (X(K), Y(-KhrK ^-i-^-fO , - ' . 



as '"implied DO^ loo'^p" potation. 

In fact; -if we look at such notation In this^way^ wfe can perhaps begin \ . 
to understand it for the firat time. Other examples which shqUld help are '/^ 
given in Table J^-1. Study these Garefully, * 

5/ ' . 

4 Table Fl^-1 [ 

Implied DO loop notation 



Exaniple "Shorthand" , ' , Equivalent "longhand" *. |\ • 



^ READ 51,(A(I), I .N) I ■ DO 10 I - 1, N, 2 * ^ ^ 

51^ F0RMAT(2F10.5.) j READ 51, A(i). A(I+i) '^x;/. 

51 F0mT(2F10.5j 



^ I 10 ' CONTINUE* 



Input values for all A^ fi'om A^ through Inclusive. Mta 

card contains up to two values per\card. If the value of N is "Jj 
for example, then h cards will be read, the laa-t ciard xiontalnliig the- 



seventh value. 



JTable Fk-l (continued)" 



Example . "^Shorthand" ' • Equivalent "longhand" 

2. READ 51, (A{I) I = 1, N, 2) ' ll>b'^0 1 = 1, .N, k 

51 rORMAT (2F10,5) J <READ 51, A(i), A(I+2) 

I 51 FORMAT (2F10,5) 

' 10 CONTINUE 



Input values f©r odd-subscripted values of A^, beginning with A^ 
and up to or including A^ (depending on whether N 'is even or odd, 
respectively). The first card contains values for A-^ and A^, the 
^ ^7^ *' ' 

3. ^ ' READ 51, (A(I). I = 5, K, 3) ' to 10 1 = 5, K, 9 « 

51^ FORMAT (3F10.5) | READ 51, A(i). A(I+3), A(i+6) 

I 51 FORMAT (3F10.5) , 
■ ' ^ * 10 CONTINUE 

• I 

If K has the value 15, for instance, then values will be input 
for A^, Ag, A^^, and A^^. The first three values will be punched on 
the first card that is read and the fouirth value will be on the second 
card. ' . • 



h. READ 51,' (A(I),B(i),C(I),I=1,N) ' " DO 10 I = 1, N 

51 .FORMAT (3F10.5)- . I - - " READ 51i'^A(l), ^(l), C(l) 

I . FORMAT (3FIO.5) 

^ 10 CONTINUE 

' I' 

Inputs values for A^,, B^, C^, A^, B^, Cg, A^', B^, C^, etc., up to and 

^ including A^^, B^^, >C^. The values ^ for A^, 13^,. for any one value 

of,^i ire punched <pn a single card. - , j 



In the rtext section we will^loob^at impli^ed DO loops once again. 



The X-field ' jr ' 

\ 

We would now like to give you a tidbit aboyt format code which might 
occasionally prove usem to you. It is the so-c^led /'X- field". There 'ar j 
several interesting uses for it. Nbrmally,' however', you can get by without" 
using them. 



f * 
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For an input format J the X- field is useful for skipping^ over (or 
ignoring) <the data which may be punched in certa^,n jolumns of a data card. 

For an output forpat the X- field allows you to insert any number of . 
blank spaces between siaumbers or comments that -ere to be ]irintecl on a line. 
When an X-fi*eld appears at the very beginning of a format code^ it guarantees 
^ cblank which can be used for single-line spacing in carriage control. 

We write an X-field code using the f;orm: , ^ 

w X 



;igned^ ^1 
iger , 

For example^ 6X or. 12X or IX (Never do we write X6 or X12 or 



an xinsigned ^letter X 
integer 



XI !) 

We now return our attention tp the table-look-up problem and again 
consider that the deck data cards containing the ,X„ and, Y^^ (Figure 
Fi^--10) shows four possible "card designs." 



Example Format Code 
1. ( (2F10.5) 



10 



20 



3. 


















. f 



10 Ik 



2h 



(I2X,F10.5,8X,F1^) 



WHO 

KNOWS 1 



WHO 
[KNOWS 2| *K 









1 • t 






2. 


(no.5,ux,Fio.5' 






\ 


WHO' 
KNOWS 


' \ 


- 1 


























12 
t 














'l 
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Example 



Format Code 
On ^ Siirt^y : 

(2F^.5) 

On Monday: 

(F10.5,10X;F10.5) 
On Tuesday: 

(F10.5;20X;F10.5) 





« 


« 






f 




















SUNDAY 


MOKMY 


TUESDAY 















\ 



-Figure Fif-10. — Vaxiuub curd^ designs illtistrate possible 
us^s of the^ X-field 



weVere pmching our own data cards for the table of function values,- 
we niight naturally choose the ^irst design as the easiest. One never 
knows how many different uses such ^data might beT^JUt to. 

For the benefit of a human reader, however/ you may have, decided to 
separate the two fields by, say, 4 columns. Of course, we would wSift" to 
tell -^he computer ^hich alsO must read this car4 tl^t ye ♦have lef^ a 
space of U Ciolumns so we now must write the fbrmat code as eitft^ 

(F10.5; FlU.p) ^ ^ , 



or as 



J I ^ (F|10.5, UX, F1Q.5) . , . 

In thi first watr we merely treat coliamni^ 11 - l4 as the left (blanlt! 
portion of the second 



field*. In the second, approach we tell ^he computer 
to complet^^ ignore ijhese 'Columns. If, som|timp later, we decide to 

four intermediate collmns, like a sequence number, 
field becomes mandatoi 



punch sbmetjiing in the 
then the use of the X 
called columns 11 - lU on the second card a 
safer to use the X-field. , - 



* That is .why we have 
HO KNOWS. It is. certainly 



iou^ 



■Here is another card design^ this time with bwo WHqpOJQWS fields, 
s'ee how we can "sprinkle" the X-fields in arid around the F- fields for ♦ 
around I-:^ield^ If *we had them ITere) i^ any way we may chQose. i 
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k. The final example suggests that several different tabulated functions- 
for the same, values of 'X^ may be punched on one card. On Sunday we 
iljight yient tq do a table -look-up using values of Y in Cqlumns 11 - 20, 
on Monday jjerhaps vlth the values of Columns 21 - 30, and jso forth. 
The same FORTRAN program could be used each time,, provided that^ach 

. ^ time- we used the program, we change FORMAT number 51 in the program 
(Figure in the manner shown (Example k of Figure F^^-IO). 

Now that we^ve begun to apjjr^eciate ways in' which the "x-field^can be 
used for input format, let's see how it may be used fcj'ic output. ,^ 

. In the first place -ve can see that for FORMAT 5^ of' the program 
(Figure Fk-9) an alternative to • ' ' 

(1HD,5F10*.5) 

would be ^ 

* . - (1X,5F10.5) 

On output an'X-field inserts blanks. So IX ^the beginning inserts 
the blank needed as a cfarri^e control signal. We can also insert blanks to 
'Spread out the printed values, l^us, if we want the 5 values to be printed^ 
with arbitrary amounts if space in between eaah ntimber, we tould write 

. (IX, _F10.5, .2X, F10.5,, kX, F10.5, kX, FIO.% ^2X, F10*5) 
Tbis fomat would group X^^, Y^^ pair and the X^^^^, Y^^^^^^ pair each 
flanking the printed value of A Jja^Ser to help the eye do the grouping as 
suggested by the wiggles below ' • ^ ^' 

vwv— WW— vwv— Wv^^vwv , ^ 

^ . ' 2 • h h 2' 

.Of coTATse, the same "sp;'ead;' coulQ be achieved by simply lengthening each " 
F-field the^ right amount such'^s: 
/ 

. (IX, F10.5, F12.5, yih.^, Fill. 5, F12.5) 

Take your choice. * ' Cj* ' 

^ The X-field could also be used to heljPrn the spacing of headings across 
"'the top of tabulated data. " • ' 
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Exercises F4-3 Set A ' ' * . 

^— — — — — t 

Go back to Chapter 2 of ydur Flow Chart text. " See Figure, 2-7 there 
; which shovs tabulated results for our frien^y loop, Figure 2-8. Your job 
- is -to write a FORTRAIJ program; which would print 'a heading across the top Sf ■ 
_ the table to "name" the values below, as. siiggested in Figure ^-11. 



O 
O 

P' 

b 

'X) 



A 



• -4.3-'^ 



B 


C 


D.- 








10.0 


3.0 • 


11.6 


2.5 


. 6.1 


; 7.9 


5.7 


-3.2 


10.7 




o 
o 
o 
6 
o 
o 
o 
o 
o 

0 



EXTRA SPACE 



•r 4-' 



, Jigure. A table with printed heading 



Print tHe^ table 'heading at the top of a pa^ and fi;feure out a way, if ^ 
you can, to^.skip a line b.efore the 'first line of numbers is printed using 
single spacing thei^afte^'. Ychi will probably have to revise the flo^; chart 
of Fig-uTQ 2-8 first. ^ ^ . « 
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Fh-k Nested DO loops ' 



^ ^,4v;Just as we can have one loop, with its iteration, bo?c/ form "part of the 
commutation portion of another loop, we can have the one pO loop jDecome part 
of another DO ioo^» • v. ^ • . 

Examine the flow eha^t in Figure i|-29, begini^ing with Box k.^ We now 
shoy the "equivalent FORTRAN code in, Figure F^^ll. 







DO ao iE 1, M 

SUM(I)- = 0,0 * 










DO 10 J = 1, 






SUM(l) = A(I,J) + SUM(l) 


10 




CONTIMJE . . ' , 


. 7 




TOTAL = TOTAL + SUM(i) 


20 
^8 




CONTINUE 
PRINT 55, TOTAL ■ 



Figure Fl^-11, '"Nested DO-loopg 
We have deXiberatsely indented the statements after each JX)-loop to 



suggest the^ idea of hes^tlng 



Sujppose we call those Statements .which are 

repeated under cfontrol of a. DO- statement the "-range" or "scope of that DO 

-Statement or D9» l^op- Using this terminology, we can^ say tiiat the^scoj)e of 

the "outer" oj^ first DO loop includes ^n6ther DO st^envent^ and its Scope. 

If you again. look at Figufe F^l-ll you see that the loubly-^indented' stjkte- 

ments constitute the scope of th§ .inner DO stat^ment^r^hile all those which 

are at least -singly- indented belong to^ the scope of^jti^ outer DO statement. 

/ . ' * • fv * ' * * 

When executing these FORTRAN statements,^ -fhe computer ^an keep track, of 

what loop it is in at all times. x. When^^the ^ 

^ ]?0 20 I 1, M ;* ' 

is executed, the computer attempts ' to, repeat iihe scope of this DO M times. 
In doing eacH repetition, i'^ encounters, the inner DO sta^tement which means 
that the scope of this loqp must te repeated* JCT times 'before moving on to 
Statement 7, and com^leiting the ^cope-of the outer loop. 



i 



, Each time we emerge^^lrbm a loop by# suc^^dihg on a loop-limit test we 
say "a DO'loop has been satis^fied," * . . » Cr* . I 

Using "this terminology,^ we might say that, "in executing the sgop^ of 
the outer DO loop we must repeat an'inner DO loop until it is satisfiea," * 



Fk-k 



Now let's look at the whole program for summing the matrix entrfes. Here 
it is in Figure FK-12. .We have assumed the matrix can have up t& "50 rows 
and 50 columns. Data for the matirix entries are assumed to be punched up 
to if ixems peF card, in row by row order. . • ' ' /• 

Label , .Sjbatement or declaration 



50 

2 

51 



•10 

20 

8 
55 



SUMMING A MATRIX* 

DIMENSION A(50, 50) 

READ §p, 'M. N ' ' * 

F0RMAT(2I5) - 

.READ 51,: (Ca(I J),;j.= 1^ n) i i m) 
FORMAT (too. 5)^ ' ' ^' 

^AL =. 0.0 , - 
DO 20 I = 1, M 
SUM(I) = 0.0 
DO 10 J = 1^ N 

SUM(I) = A(I,J) + SUM(I) - 
CONTINUE 
TOTAL = TOTAL + SUM(l) 
CONTINUE 
PRINaj*^5,) TOTAL 

FORMAT (l^HDTHEaTOTALalSo, FIG, 5) 

STOP ' V 

END 



Figure Fi;-12*. 'The whole show 



Nested implied DO loops 

Letfs take a moment here to digress again on the' subject^f implied 
• DO loop.notation. In Chapter's we introduced you to the flow chart, notation 
like the ohe on Box 2 (Figure k-Z^) . 

'/'■■. ' . ■ 



The FORTRAN way/fs to write something veiy similar--nam^ly, 

2 READ 1, ((A(I,J), J = 1, N), I = 1, M) 
inner implied DC-loop 



outer ir/plled DO-loop 



. 'Youjcan now begin to see^ the method in this notational raadnesS'-if you^ ^' 
want to take 4;he time..^ . ^ 
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This is seen to be just a , slick -shorthand yithin a shorthand to represent 
.a simpler READ statement under control of 2 (nested) DO-^loops, Figure. P4- 13, 



70 
10 
2Q 



iX) go I = 1/M 
DO-10 J = 1^ N 
REAQ 70, A(I,J) 
FORMAT (F10.5') 
COmiNUE 

coiyrmiE - ^ 



^ . Figure Fi^-135r The "long, way", but not q\4t,e .equivalent 
^ ' * - • to Statement 2 ' ^ 

/ ^ In this case ve show here the equivalence with Stal^enient 2 is rtet quite 
complete because, only one matrix ent^ry can be put on one card whereas in ' 

- Figure 'F4-1§ v6 v^re^^able to specify \ items jDer card using FORMAT. We 

can correct this by having \ list elements in the" READ statement as shown 
^ in Figure FI4--IU. * • , ' . ' 



\ 



51 
>0 

20' 



DO 20 I = 1 M a ' ' * / 

DO 10 J = 1, N, ^ 
• KEAD 51, A(^ J), A(I, J+1), A(I, J+2),^A(3i,J+3) 
FORMAT (too. 5) 
CONTINUE • : 

COOTINUE ' T ' : ^ ^ . v ; 



Figure "The long way", but equivalent to 

St^ement 2 - ^ ' 




• i 



There is ?till ailotiier and even better "long way" as shown in Figure 
-F'^-l^. would recommend itoVTsy-^any;, ^ the other methods. It has both tl^e' 
' «,hort$ ,and the long look-. 

V . . . 

DO 20 I = 1, M 

READ 51, (A(I,J), J = 1, N) 
* 51K I * FORMAT (too. 5) 
20 1 CONTINUE 



^i^ureJPV-lJ.^ -Thfe "long-short" equivalent to 



"^co^ee'that we Jiave used a. single implie(^nDO-loop^in Vc^^ READ state- 
' *. ment and w^'rSpeei.,^e read statement mder control of a, DO statement. This 



method has the ^following virtue: 



The READ statement amounts to ao order to ready's whole'row of the 
matrix with up to four entries pefr card. If theilast entries of .a 



row don't quite fill' out a card— np matter— we can still start .the" 
ent-ries for the next rov-fr^esh 'on the next card in sequence,. This '\ 
way if- you have a fairly-big matrix to put otj cards, it's ^asy^to, f 
verify you've, punched entries .for each row correctly. You don't get ^ 
this, flexibility when you buy the, doubly, iinpiied DO, loop that is 
used in Figure Fk-12; there entries for each new row begin; if "neces-^ ' 
saiy, on tlfe same card with the last entry or, entries of the preceding 
card. * . f ^ 



Exercises Fk-k Set A 

1 - 8. For each of the flow charts^ycu constructed for^the exercises Xn 
Section h-k, Set A, write the equivalent FOI^TRAIf stateme?its. DonH 
bother writing declarations unless you fe^el/they ^dd.to your under- 
standing of the transliteration problem. 



-^Triply ^nest.ed DO loops ^re just as easy to* write as doub'ly nested 'ories. 
Figufe Pii-lg'^sliows hoyj the stickle'r in Figure 4-31 -would be coded in FORTRAN; 



Label Stateaient 



5 
50 
10 
20 
30 



DO 30 m = u 9. ^ 

D0'.20 IT = I, 10 ^. 
BO 10 Ili = 1, 10* * * 

ITMl = IT ^ 1, \. 

^ iim = lu - 1 

IF (100 + lO'^-ITNa + lUMiy 

7(IH**3 + ITM1^3 ^ lUMl*^) 
.10, 5, 10 - 
PRINT 50, ^H. ITMi; lUMl . 
FORMAT (3115; 
CONTINUE, ' 
CONTINUE 
'CONTINUE . 
STOP - ^ • 

END 



lFigiir6 The stickler in FORTRAN 



99 
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Since ^ DO-loop counter may not .begia witb ^^tero vi^lye we have had to be - 
a bit inventive to« accomplish the same objectives *of 'the -stickler algoritlim. 
As you.can see in Figure Fk-lSy we let. IT and . lU rang^^from 1 to 10 
insteg^of from 0 to 9. In the IF statement expression we use lOKL 
-(•def inedf as IT - l) and lUMl (defined as lU - l) *iij place of IT and " 
lU, respectively. ^ * 

Exercise F^-^ Set B • , * ' ^ ' 

I. .Write FORTRAN programs for the flow chart solutions ^ou olitaihed for 
Problem 7, Section i|-i4-/ Set B "of the Main Text. , ^ ' ^ 

Exercises i'k-^ Sfst c • ^ . ' *. 

^^7"- . . . ^ ^ ■ ■ 

1. Write a complete FORTRAN program for the Complete factorization Al^orithm^^ 
Figure Assume the input value foi» N* wili^ not be as large as^^ ^ ^ ' 
ipOOO. In convert;Lng to FORTRAN, you should" be on the lookout for^tvo ■ ;> 
interesting features.. First^ take aim on the test-'ijjortion. of the iter- 
ation box--the upper limitjnust be an integer. Second, notice that' th'e 
scope of the DO loop includes a "complete loop--but ^witliout ari iteratipp" 
box. ' ■ ' » , « ^ ■ . 

2. ^ , Write a complete FORTRAN program for the shu;|;tle-interchange sorting 

, algorithm^ Figure k-^h,. Assume you may wish to :^rt up^ 500 numbers, 
each of wKich.can be input' using kn FIO.5 f^i<ld— say up. to four values ^ 
■ .per card.\^^^ on the lookout for oc^^a^ional difficulties in converting 
. the flow cha^'- to the FORTRAN^^code--.especially in modeling Box 6 which _ 
\ shows d"ecrea«sing,.c6iinter s,et fhitially '^fto an expression. ' • 

3. Write a .aompletis FORTRAS^- program fgr.the sart algoiSithm sh6Vh 'in Figure^ 

' ' ' ^ » j ' \^ ' ' 

^4-35. Maker the same assumptions In this program ,:^M-t you are asked 

to make in .the preceding exercise. • 0 * 

k. Write :a compl^ifc^f^R^RAK program for finding the- largest decreasing 



. subsequence^ ^ase your program on the flow chart in Figure ^^--'39. 

Assume .th^ giveh sequence will not exceed ^100^ values' iti a/id that ,^ 

,th| 'values lih^elves -are to be punched on dat^^iAsy fb;ar values per 
card, governed by. Fio.5 field codes. ^ -.^ * 
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Chapter F5 
SUBPROGRAMS 



J 



^ P5*-r .FORTRAN subprograms * * • * 

FORTRM programs that correspond to reference flow charts are called 
FORTRAN subprograms. Ab ydu might expect, these subprograms loc^Vei^ much' 
like other IF'ORTRAN programs except for special statements at ;the beginning 
arid end corresponding to the funnel and return box of the 'flow ehart. In 
addition to these special statements there Ire^ ndmber. of conventions to be 
observed which will-^e explained in what "Allows. ^ . ' * 

A scabprogram is a self-contained unit which'^can be compiled separately/ " 
and then used by a nmber of programs or other'subprograms . One of the'^main J 
. advantages of separate comjnlation is that one can in. this way establish and 
develop a library of subprograms, which is av^xlable for.Mter use. . ^ 

Subprogram? that evaluate 'a function and report ^ single value to the 
main program are called functiog subprokrams . (A ^cond typ^ of subprogram 
will be encountered In Section ?5-i^. ) /For^ a- function subprogram, the i:iow 
, chart funne^ corresponds, to the first statement of the sdbprogram. This 
. . ^tatem^t begins. yith the^vord FUNCTION, ^ followed by the name . of ^the . JTuhdltJ^m 
and itSfar^^ent in parentheses, r'or example. Figure 5-V corresponds to * 

FUNCTION SQR06t(y) . - • ' ^ * 

Sinee special symbols iike -^TT are not available in FORTRAN, we' 
naturally replace such a symfeol witj;; an alpjiat^etic jiame for the function (in 
this case SQROOT). Names of function subprograms are to be chosen with some 
following cVventions. ,The firs,t character in M^he sub- 
* fwgram name must be I, J, K, L, M, or N if and only if the v^lue reported ^ 
is in the integer mode. The last character must not be ' F if the name has 
more than 3 cliai'acters-. 'otherwise the name can have ' 1 to 6 characters, 
the first being ^alphabetic. ' * ^ ' • 

You are alrea^^ f&mili^r 'with the use of predef ined^^iat hematic ai ^ ' ' 
. Ainctions-'like Sq^F, SMF, ABSF, etc (gee Table F2-2, ) FORTRAN sub- ^ 
".'♦Programs are different from these predefined mathematical functions even 
though our first example, , SQROOT, > served th^ same purpose ^s SQRTF. The 
^ difference is, hot just in the way the names are spelled '(the conventions are 
different fbr subprograms and .predefined functions), but. is mainly in thife ^ 
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F5-1-" 



fact that*the predefined mathematical functions are part of-/ the compiler 
system. The techriiques for^ adding, to the list of , predefined' mathenjatical 
functions are outeide of the scope of this book^but :&unctj.on subprograms " 
provide a way that, you can use to develop vhatevjgr set o^' reference programs 
you want. ^ * * I / ^ 

Just as reference flow charts must sooner or later' reach a return tipx, 

♦ / ^ • 

a subprogram^must eventually reach a RETURN statement.^ This statement, 
corresponding to Figure 5-5; ^consists just of the w^d VretURN". This state- 
ment need not appear in the^^Last line of the subprogram .(more than one RETURN 
statement may be used in a sing'le subprogr^), but it must beN;he la^t 'step 
in the <exec,ution of the ' subprogram. . * . ; ' 



Since the RETU^ ^tatement of FORtkain aoes not indicate what variable is 
to be reported to the n^ln program, another convention is needed -to identify 
the value to be reported. The convention is that the name of the functioit 
subprogram itself must appear ^at* least onc^ on the lef^ -side of an assignment 
statement and the value assigned is, upon return-^ the value reported to the 
main program. For this purpose only, th^ function name fs %b be^ thought of 
as a variable. . ' • » 

Tti all other ♦respects, a FORTRAN subprogr^ must conform to the require- 
ments of any ftiRTRAN' program. In particular, the l^^t line of 'a subprogram 
must "be folloveS by. 'an statement,- and ^JMENSION infomat^ion miist be given 
for each*^siib scripted Variable in .the subprogram. ' ' . * 

TI)e use of a 'function subprogram by a main program is shown in 
Figure F5-1. ' ' ^ ' . . > ' ' 







1 . ' 


0' 


1 


Z = A + SQROOT(X)-- £i3^ 


-j^ FUNCTION SQROOT(y) 
i • - 




■3 


IF(Z) 5, ^, 5 -/C " 


1 

1 - - . 




k 


■ z.= A - sqr66t(z)' /T) 


1 




5 


^^G) — - 


1 ^ 

"|-SQROOT = 
-pRETURH 






fi , 


1 END 

_P 


I — 




FORTRAN main* program 


1 function subpi^)gram 





^Figure F5-1. Use of a FORTRAIJ .function sjabprogram 



. Figure F5-1 ^is ititended to xiorrespond to Figure 5-6 of the flov chart 
text and to illustrate ohly those features related* to subprogram use. The 
first time the su)$prograin is required (in Statement Nuraljer 1 of the ,main 
program), ve go to the FUNCTION statement via Route 1. This statement direct 
that tije value of X be assigned to Y. It is essential that the name of. 
the argment^n the FUNCTION statement and the hame of the argument vhere'the 
subprogram is reqvieste^d be of'th^ same mode. That is, since Y represents 
a real variable, a request for .SQROOT(N) would^ be, in error si,nce N repre- 
sents an iijteger Variable. * ' ^ . 

When the subprogr^ has been executed a valuW-ha^^been assigned to 
SQEOOT and the return to Statement 1 of the mai^l'|)ro^|^am is by^oute 1. ' 
Upon return ,to the main program/ the value assigned to't SQROOT is added to 
A and the. result assjigned t® Z. 'Where the subprogram is again required - 
(in Statement Number K) ve are to go to'the FUNCTION stat^nt by Route ^, 
assign. to the Y in the subprogram (note they agree^^^oJe),,; execulSe 
the suUprogram, and return via Route* 2* to Statement i^. vitr the result of the 
subprogram assigned to SQROOT. / * ^ ' 

An actu'lal function subprogram for the square root can be prepared vith 
reference- to Figure 5-7, as shown in Figure F5-2. . . . ' 





C 


SQUARE ROOT 'FUNCTION SUBPROGRAM . 






FUNCTION SQR00T(Y) " ^ ^ 






G^^ 1.0 • 






H = 0.5 * (G + Y/G) * 






IF(ABSF(H - G) - .0001) 5, lO, 10 




' 10 


G = H 






GO TO 2 






SQROOT = H ' ^ . . 






RETURN 


r 




-END . . . • . ' 



Figure F5-2. Function subprogram for square roots 

• ■ ■ ' . ■ .* ' • • 

As you inspect this first function subprogram we take the opporturClty to warn 
^ ypu aw^^rom something you^re not likely to do ar^^how. ¥hdi;lis, don'.t try 
to assign anything to an ^gument (in'this case,* Y) of a function subprograi^. 
It is never necessary to do, this anyvay. With some FORTRAN compilers this 
will produce very strange results— so let the buyer bewarel (You will be 
given one more reminder of this danger area in Section F5-3--then you will be 
on yoTir own. ) * - * * f . * 
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Exeroises ■ . ' • * 

^ ' 1 _ ' / 

1-3- Write FORTRAN function- subprograms for the- flcTw bharts prepared in 
Exercises -5-1 mainHext. 
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F5-2 Functions, and FORTRAN 

^ The -f loH ^*h&rt' text' tells us that any flow chart which, when given a 
value will ]^oduce another value, can be viewed -as ^e.^ evaluation of some 
function. The statement is true for FORTRAN programs' as well as for flow 
charts. Mathematical functions exist^i^iich cannot be evaluated with a flow 
chart or by a FORTRAN prografa, but the common usage of the word funcjtion- in 
•computing is strictly limited* t'o those 'which can be evaluated with^ a flow 
chart. Moreover, we do ."not -^^Uy think of expressions which can be precisely 
evaluated solely by one 'of ^e &)ur basic arithmetic operations as functions 
althou|h, fact,; they arellh cbmputing, then, a function is commonly 
^th6ugh* of as a' relatio/isViip for which^^a reference, flow chart is used 

, • The,domaifK)f a function, i/ c'ompuling, is the'set of values that the 
arguaient,in*th-e fUnnel -of .the flow chart can take on. The range of a function 
is. the set of values that-can be reported to the main flow chart. In FORTRAN, 
the d6maln^can be either a set^of real numberV representable in a , computer or 
.a set of integers representable in a^computer. Which set of numbers is mean}; 
is indicated, in the usuarirvay, by the fi'rst letter of the name of the arg^- V 
ment. The range, too, can be' either from the set of integers or-.from the 
set of real numbers. Wilch is meant is indicated, in the usual way, by the 
spelling of the na^e of^the f^ijctfon. ■ ^ ' ' ._ 
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F5-3 FORTRAN^ fungt ions viith more than one argmen^^ ^ * ^. ^--w I ^ 

FORTRAN function subprograms chn have as many arguments as are nec.essary. 
The min function provides an exaqiple. The subprogram in Figure F5-3 corre- 
sponds to Figure 5-1^1 of the flov chart text^ 



MINIMUM OF TWO ARGUMENTS 

FUifCTION SUBPROGRAM 
FUNCTION FMIN(B, G) 
IF(B C) '2, 2; 3 
Z B , — . 
GO TO k 
Z = C 
FMIN = Z 
RKTUM 
END . 



/ 
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^ * i'igure F5-3. A function subprogram pf two' argumerrts 

Notice that we have changed za^ name of the function subprogram so^that it 
do^s^not begin with the letter M. That is^ the function subprogram, expects 
to ^ceiv^ two real values and to report^ real value. We do itot need to 
introduce the 'variable Z since< FMIN can serve tKe same role. For this 
reason^ Figure F5-3 can* be r^^ed by Figure F5-^.' Should Figure F5-2 b"e 
changed in a similar way? - F 
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MINMJM of two- ARGUMENTS 

FlAcTlON^'SUBPRC^RAf^ 
FUNCTION miU (B, C) 
IF(B - C) 2, 2, 3 
MIN = B 
GO TO i| 

FMIN = C ' i 

RETURN • . ' " ' y 

END 



* . Figure P5-^. Improved FMIN subprogram 

* The parameter list of a FORTRAN subprogram can contain, ijit eg er ;variables, 
and constants^ real variables and constants^ variables coataining alphanumerip' 
information, ^ the n&me^ of vectors, and^ the names of matrices. In every case 
there must be a one-to-one correspondence in the number and types of para- 
meters "between those of the parameter list where the subprogram is requested 
and' those of the parameter list in the J^CTION sta-fiement. - ^Confusion would 



/ . 

^ ^reign if we t;:ie*d to request the^subprogram^of Figure by writing some- , 

' " nhing^li'ke' ' ' " * ' " , ' "^^ " " ""^-'-^-^ ' 

' ' . T = PMN (A, -B/ C) ^ . ^ 

or , T = FMIN (M, P) 

A elaBsif ic ation of variables 

« * 

The distinction between^local and nonlocal variables Is- the same in 

• FORTRAN! as is d§s.cribed in the flow chart ,text. Thus in the two subprograms 

for' FMIN given earlier, the argiunents B and € are nonlocal variables. 

^ In the first subprogram Z is a local variable. 

Some compilers are written on the assumption you will nev^i* try to 
reassign values to nonlocal variable^ of a FUNCTION subprogram. Such compile 
^take advaxit^e of this fact in assembling, efficient t&rget code for your 
program: Hence/^an attempt to "out-fox" the compiler in this respect, just ' 
to save a -step in the program, or to save a storage cell or two, cotlld lead 
to trouble. Unless you are absolutely sure about what can ha:ppen you should 
avoid changing nonlocal variables in a subprogram. 

. An appreciation of.^he distinction between local and nonlbcal?tvariables 
can be had by recalling that subprograms can.be compiled- separately from 
programs that make use"-of *them. During cc^ilation a local variable is ' 
associated with a specific location 'in memory (except for an additive con- 
stant needed to account for possible shifting of the whole subprogram from 
^one place in memory to another.) On ^he other hand, . there'ls no way a sub- 
^ --progr^^can know prior to compilation of the main program, where nonlocal 
^ variables/may be located when the subprogram is executed.- 

Independence oP statement ^ numbers among programs » . . ' 

' ^ The realization tl^at compilation of a subprcJgram is an entirely separate 
^process from the compilation of a m^n program also brings out the fact that 
statement numbers in a subprogram are entirely different references from ^ 
statement niimbers ir^ another subprogram or in' a main pro-am even if the same*" 
numei^^l-has been used. ~ . , . _ . . . , . f . , ^ W , ' 

Composition of functioft referenc.es 

^ TJae situation wi^h respect to composition of function references is 
exactly as described in the flow ehar't te^t. Th±&^ just wh^t you havfe 
learned as composition of functions. That is, given function subprograms 




* >> » 

■ ' ■ 

defining Fl(X) and F2(x), ' one can write • 

f ^ Y = Fi(F2(X)> - • ' ; 

so long a's the Vange of F2 is a subse-^ of the domaia of ^ PI. "Correspondingly, 
the following function referenc^es are entirely |)roper: . t 







Fto:N(ABSF(A + B), 5.1*) , 


or 


xY = 


MEN(FMIN(F, AB^(T)")7 '~Qh— — . ^ 


. ^ 6r 


Y = 


FMIN(SQROOT(B * B - U.O * A * C), -B^ 


* or 
* 


Y = 


SQROOT(FMraCX,Y)) * • 


> F5-3 Set A . 







1-7. Write FORTRAN programs &nd function subprograms for the fJ.ov charts 
f prepared in Exerdises 5-^3, Set A, Main Text. ^ ^ ' , 



" Exercises F$-3 . ^et B 

1. Write a FORTI^N function- subprogram for the" GCD *functiuual ;*eference 
whose flow chart ygu prepared in. Problem 1 of Set B In the main i^t. 
can this subprogram KGClf; Why. must the name be changed for FORTRAN? 

2. Write a. FORTRAN function subprograiri called KGCF cofresj^pnding to the 
flow cljart you drew /or the GCF algorithm (Prob^lem 2, Sejb B, main text). 

•3. Write a FORTRAN \rogram. that corresponds to the flow chart for determining 

(a) the number of similar triangles 

^ ^- (b) the total perimeter of similar triangles *^ 

corresponding to the flow charts you prepared In Problem 3; Set B, 
main text. 

V • ■; ° " ■ , ' * - ■ 

* ^k. ' Write a FORTRAN program .that cojq^esponds to the algorithm for Problem 4, 
Set fe, main text. Try to e'stimate how much computation wlll*be ^involved. 
Computation can be measured in terms of the number of ^additions, ^ub- ' 
tractions and comparisons that must be made, counting each- as 1.^^^ 
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These problems are quite difficult. Ihe_ student vill be able to solve them^ 
only with considerable time arid effo^. . • > ' 
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F5-^ .>|ORTOAN procedures 

V. * ^ . 

^ ^ FORTRAN subprograms which correspond to reference flow charts for pron 

cedT^res are called SUBROUTINE subprograms. "Corresponding to the funnel of 

the^ reference floV chart is the .SUBROUTINE statement which takes the form 

suQ^h as: \ 

• • ^ • ^ * ■ SUBROUT^jPN^"sORT(N, V)^ 

The; subroutine statement consists of the word SUBROUTINE followed by the name 

o^^he procedure you are defiiung and a parameter list in parentheses.. - In. 

tliW SORT example we see that a^^.^ent ire vector is identified solely by its 
xime, V, in the parameter list. No attempt is ^lade to 'subscript V in 
the paramefter list but the subprogram will have to proyide DIMENSION inform-' 
.ati^ for this Vector. 

A^-v ^^^^^ ^ subroutine subprogram does not report a value in the same way 
^ffihat^ a 'function subprcJ^W does, the name of the procedure will 'not appear in' 
^''-tfie TOdy of the subprogram. Moreover, since this name will not be used as 
^^though it vere a local variable, its spelling (i.e., its first letter) has 
^no special significance.* The name of a procedure musf^not end in F if 
^^it is more than three characters long and all other ^conventions, are as they 
•^are for fmction subprograms". 

A subroutine subprogram corresponding to Figure 5-l6 Is given in 
' Figure F5-5. The -dimension of V has been given as 100 but, for each use, 
N is the actua^ number of confponents in the vector to be sorted. Thus, this, 
subprogram is usably for any vector having 100 or fever components.. 




SUBROUTINE SUBPROGRAM FOR SORT 
SUBROUTINE S0R*(N, V) - 
DIMENSION V(IOO) • 
K = N - 1 
.DO UO I = 1, K . 
M = I + 1 
DO 1^0 J =''M,'-N . 
IF(V(I) - V(^» kO, kO, h 
B = V(J) 
V(J) = V(I) 

yC?.) = B . 

'CONTINUE 

RETURN 

END 



Figure ^F5-5. Subprogram for sorting 
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oxes. Notice alsa that valued 



Look ^t Figures 5"l6 and F5-5 side by side, 
correspond almost exactly to the flow chart 1 
can- "be intentionally assigned to parameters cif a subroutine ^subprogram .- 'We 
warned you not to|do, this in function subprofirams h\xt this 'is 'Vow -a subr6utine 
subprogram produces its output. 



See how the FORTRAN ^atements^^ 



r 

at 



Use of a SUBROUTINE subprogram is accomi)lished with a CALL statement, 
% ^ '1 

analogous to the e:?%cute box of the flow chaH text. The CALL statement 

/ ^ 

consists of. the word CALL, followed by the name of the subroutine subprogram 

and a parameter list in parentheses, for example: ' 



CALL SORT(88, B), 



where B is dimensioned to have at least 



components in the main program. 



.Im FORCIHAN, the process' of referring to a sub'routine j.s called " calling the ^ 
subrout'ine " ' ^ , ^ ' , ' - 

A main program cabling SORT, corresponding to Figure 5-23> is shown 
in Fi|ur« 55"6, This progiram ^ , ^ ' * ^ 



101 
102 

103 



A PROGRAM TO ILLUSTRATE CALL STATEMENTS 
DIMENSION 3(100), C(lOO) * 
READ 101, K * 
FORMAT (13) 

READ 102, (B(l), I =;l,"K) 
FORMAT/5FI5.8) 
REAR 102, (C(I). I =" 1, K) 
CALL^SORT (K, B) 
CALL SORT (K, C) ' 
PRINT/103, (B(I), C(I), I = 1, K) 
FORMAT (IX, 2F I5.8) 
.STOP ^ 
END . ' \ 



Figure F5-6r A main J)rogram cal|Ling the sort procediire 

assximes that, is a three digit integer jin thp first three col^mns of the ^ 
first' data card; that the components ^ B and C are punchei^ five per -card 
in F15,8 format.' .It is also assumed that ^ '"'"•'^ ^^-^^ ^'^.r^oorj im\ 



K will not exceed iOO\ 
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Exgrc ises ¥3-k Set A • 

1_H 5. Fo* the flow charts 'prepared in Exercises ^-k, Set A main text, 
wite FORTRAN calLing programs and subroutine subprograms. N 



Exercises F^-l^ ( Set B 

1 - 3^JVj^5jjaow^-char-6^e^^ in Ex'^ire7?3;7set B main text for. 
'Problems 1, 2, and 3, write th4- FORTRAN ^function and subroutine sub- 
programs. | 



Exercises F^-h, 



For the 
Problems 1 



Set C 



flow charts jtrepared in Exerciser Set C maia text for 

through i, write the FORTRAN procedures and programs. In 
doing the ijrogram for the subroutine DEGREE, ke4p in mind that you must 
shift -the indexes for the coefficients, i.e., the coefficient a^ .must 
be associated with A(l), 'etc. The same problem must be handled in a ' 
similar fasiiion in coding the flow chart for SIMPLI^ and REDUCEMOD. 
(spelling O.K. for FORTRAN?) in Problems 2 and 3.' ' ' 



I i 4 
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F5-5 » Alternate exits and procedures for branching / 

Provision for alternate exits and branching from subroutine sulJprograms 
_4n^F0R^AK^mjxrors-J:h^ -di-scus^ioK In the" iTbv cttart'Text'; parameter 

is provided to indicate the result .of tests performed by the subprogram. 
Figures K5-f and F5-8 present programs corresponding to Ime *flov charts of* 
Figure 5^26 and Figure 5-27' ' / " 



SUBPROG'RAM to test EQUAlilTY OF 

TWO COMPLEX LUMBERS / 
SUBROUTINE C0MPEQ(A, C/ D, N) 
IF(A - C) 2, k I 
IF(B - D) i;, 3, ^ /' i 
N = Or- 
GO TO 5 ' 
N = 1 . 
RETURN 

END / 



Figure F5-T. ^^ua^ty of c^plex numbers in RORTRAN 



^PROGRAM SEGMENT- SHOWING USE OF COMPEQ 

C/^LL COMPEQ(X, Y, U^V, K) 
.tF(K) 3 
STATEMENT 3 IF UNEQUAL 
STATEMENT h IF EQUAL 



Figure A program segment to test equality of^ complex 

/ numbers 



Statement labels and function names as arguments 



FOJiTRAW II was not designed with such types of arguments in mind- 
fev if any processors have this capability, Mpre advanced languages c 
however, pennit full freedom in the use of such procedure arguments. 







• 
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Exercises F5-5 


• / 


Gene^ry^^^^f-hy^hsT^r-sT^ 


. Text; and Jhen^write-a FORTRAN ^subroutine s^l)progra2n to solve two . /' 


^ equations in 


two unknowns. " * , ^ . : 


\ ^ ^ . , ; * ' J . * " 

^2. Use the appropriate flow charts prepared for Problem 2, ExeVcisesl 5-5. to 


write a FORTRAN program and subroutine ^subprogram which will compute th6 


real roots of 


a quadratic equation* .Use the prqgr§m^to s.olve the follow-' i 


ing equations 






2x^ - 3x + y = 0 . ' ' 


• 


2 . ^ • 

■'^ 3ai^x - 6.2x - 11^.23 = 0 ■ 


i • » , 

3* Write ^a FORTRAN subrouting subprogram and calling program sta'bements 


corresponding to one of the' techniques you used for solving Problem 3, 


&cercises 5-5 




> 

f 


' ' " * \ 

f : ■ 


• 




■ . • \> 




' ' ' 1 ■ 


9 ' . 


^ 
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F5-6 Symbol manipulation in PORTRAN ■ 



J 



In Section F2-10 ve discussed .j^e' input and output of alphanumeric ' 



charaA:ter>s, No^ we want to find out how alphanumeric, data^ can De processed 
so that we will be ab^le to alter such input data^ as 

THE QUICK BROWN FOX JUMPED ^VER THE LAZY DOG, 
or 3.1^^159 ' •* 

or r + s^t + u(v + w)), ^ 



'Since we will wan1> to be able tfc^refe!' 4:0 each individ\^ element 
in^such character strings, we will associate a s^gparate variable with 

each ^element of a. string. This means ^hat a card, let lis, say, would 

I 

be read vith ^ FORMA^* code of 80A1. \ 

We are now ready to code a sub^putine subprogram, -.gor responding 

to Figure^5-33, for CHEKCH. !rhis_progfam_is -s^own in Figure ^"9* 

We notice that there is a ^DIMENSION statement giving the maximum length' 
of 5 as 200. Other t^ian for the DIMENSION statement the program 



A -SUBPROGRAM FOR'.GHJKCH " 
SUERDUTINE CHEKCH (N, S, M, C, IPof 
DIMENSION 3(200) 
DO 2 I = M, N 

IF(S(I) - 0*2, 3, 2 y ^ - 
CONTINUE ^ • 

IP = 0 

GO TO 5 ^ 
IP = I 
RETURN 

END , . ' 
' 1 



> Figure F5-9/ A FORTRAN, subprogram for CHEKCI^ 

•could accomO(fcte any length string. ^As presented in Figure F5-9; 200 is 

simply a maximum length and could be replaced by a^arger integer subject only 

t9 the amount of memory available. The' alphanumeric^ .data is associated with 

real variables, S and C,^ ;as- Irl Section F2-10. Finally we see that two 

characters like the values of S(lJ and C can be tested for. identity by 

» 

subtracting one from the other as if they were numbers, and testing for a 
zero difference* * *. 

^ . A FORTRAN subprogram corresponding to Figure is given in Figxire 

I5~10. Here we assume the string, S, has 200 or fewer elements and the 
substring," C, has" 20 or fewer elements. 

ERIC- - • li;/ 




A SUBPROGRAM FOR CHEKST 

SUBROUTINE CHEKST (N, S, % K/C, IP) 

DIMENSION S(200), G(20) 

~TW(r^jr+ir^y~3, i, n — 

CALL CHSKCH (N, S, L, C{l) , IP) . 

IF('IP) 5, 11, 5 - * 

IF(IP - N + K - l) 6, 6, 11 ^ ' 

IR IP + 1 

DO 9 J = 2, K 

IF(S(IH) - C(J)nJ;10A9; 10 

IR = IR + 1 ' : 

'CX) TO 12 - ^ 

l\= IP + 1 ' . ^ ' 

G(? ?0 2 * / - 

IP = O 

RETURN 

END 



"^^figure F5-10. A FORTRAN subprogram for CHEKST 



As the ifloii^/chart text has pointed out, the ability of* a subroutine sub^ 
'program to calX another subprogram has' far-reacRing significance. It is this 
ability that p1&,rmits the construction of subprograms of increasingly greater 
comple;?ity rronl^ simpler "building blocks". Chapter 8 and all larger problems 
* will make heavy use of. *this^uilding block 'ability? 

Exercises ' 

J. ii-. For the flow charts prepared in Exercises 5-6 , * maitjT texi:, write^ ^ 
FORTRA{f .|Srograins ^nd subprograms. 



' One fiarther point should be made here. It is often' very con-^enient ,to 
consider'th^Clength of a string, and the string itself as being a Single entity. 
One. way to-Ub this is to u£« an array, say STR, with the property ^h at its 
'first component STR(l) is equal to the length of the string, while its-re- 
main^.ng co.mj^nents STR(2), STR(3), and so on, are the characters themselves. 

r If we^denote the character* for which th^ search of CHEKCH is being made 
'by CHAK, '/"^nO. its positionr-(if found) by KP, the subprogram is shown in 
Figure F5.^11. . - ^ 
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F5-6 



1 



MODIFIED SUBPROGRAM F0R CHEKCH 
SUBROUTINE CHEKCH (STR/ M, CHAR, KP) 
DIMENSION STR(IOO) 
MM* = M + 1 ' 

jL^TR d) + 1 ^ :i \ 



DO 2 I = MM, N 
IF(STR(I) - CHAR) 2, 1, 2 
KP = I.- 1 ■ 
* GQ TO 3 . ' ♦ 

CONTINUE > 

KP/- 0 . .; i ; . _ ' 9" 

RETURN • ' , 

END 



Figvire F5-11, A "modified subprogram for CHEKDH 



The only complication is in adding 1 in two places^, and subtracting in 

another. These* steps are required because now the I . character in the 

th • ' 

string is the (l + l) - *cbmpon^nt of the array STR, ' ' • 
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Chapter;' P7 
SOME MATHEMATICAL APPLICATIONS ' 



F7-1 Root of .an eq^^ation by bisection 

*Now let us \a*ite a program find the root of%i'equ4tion ' y = f(x) in 

We s'hall write^the pilgrim, corresponding 



the interval = A 



and 



= 3. 



to the fl^v chart of Figure T-S*, 'in the, form of a subrou/ine.^ubprogram called 

^ZERO so tiat it can be used with many "different main programs*. Then ve ^hall . 

vrite a FORTRAN program whiCh calls ZERO. • ' ' 

J- 

. The given function f(x), ' ye shall call FmCTW} it will need to be 



acc\iracy of^ihe result 



defined as a FORTRAN function subprogram. Th^ desi 

is given as EPSI; so the bisection process will *be/ terminated when the 
length of the interval is ^fess than the value of^ E^SI. Inside the: subroutine 
the value of the root of the equation will b*e assigned to the variabJLe ' ROOT. 
The final yalu^ of ROOT .is p;'inted before ^turning to the calling ^ogram. 

PiguafT^-l shows one way to code thetZERO subroutine. Only three dummy 
arguments are used, ^ A,'B, and EPSlJ' ,mie first thing that's done by-the 
subprogram is to reassign the Wues of A and B to auxiliary varifWes 
n and X2 as a safeguard, to protect the values of arguments that m^ch A 
and B in the calling program. In cJase th4 bisection method is inap^^able 
we let the -subroutine print the messiage* ' * ^ 

The value of. ROOT could, of co^rse^^be carried back to tlie c^ing ' 
program by makihg ROOT' a dummy- varS^l^^and adding it to the parLneter list " 
in the 6penirig 'declaration. ' In this event the CALL statement would also show 
a parameter list of four actual argument^, .the last one being a variable that , 
:cepts the'computed value, of the root. ^ 
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SUBROOTINE , ZERO ( A , EPSI ) 
XI ^ A « V 

X2 = B 

n p FUNCT(Xl) 
IF(Y1 - JFWCT(5C2)) 6, 5^ ^ 
IF(Yl) 8, 7, S"" 
'JWT ='X1 
GO- TO 12 
ROOT = X2 

PRINT 6lX ROOT * 1^ 

F0RM4T(E20..8) . ^ . • 
RETURN 

(XI + 
- X2) 



X2) 

- EPSI) 11, To, 10 



XM = 0.5 * 
IF(ABSF(X1 
ROOT = XM 
GO TO 12 

IF(Y1 * -FUNCTCXM)) 13, lOj Ih 
X2 = XM- 
GO TO 6 
XI = XM 
GO TO 6 

PRINT 62 , ' . 

FORMAT( 23HnMETHODQISQlNAPPLICABLE) 



RETURN 
END 



Fi^re F7-1 \ 



E-fields : ^ . 

There are many things you probably want to. looi:^ at carefully in the Z§RO 
subprogram. We ijiight, however, digress temporarily to explain ihe strange 
^^^^JCield code that is used in- format ntmib^r 6l. 

The "E-field" ijy* in many respects s^imilar to the F-field for itiput and 
output of real data values. E-fields are used to describe the inj^t or out- 
put of real data values vritten in the* "E" notation like 

' ' . • ^ 0.5253E-5 / - 

, A number jjrinted under control of an E20.8 field' would occupy 20 columns. 



'If t] 



vaK( 



ue in memory is '-1.3^4-67 this number wo\jld appeeC): in print as: 

\ 

□ 000 □ *-0. 13i^■67000EDOl 
^ — * ' 

8 places 

■ 20 columns » 



The right-most four .columns would bemused for printing the exponent, the 
first of which is for the letter E,. the second for a minus sign .or blank> 
- ^nd th$ third and fourth for a two digit exponent (poVer of 10). There is 
a^zerp to the left of the decimal point of the" precision part preceded by a 
minus sigrf if the number is negative. The leading digit on the right of ^he 
decimal point is always mkde non-zero and the exponent is adjusted accord- 
ingly- ^ , ^ ... 

A number that is read into memory ui^der control of an Erfield like 
need not have the decimal point actually punched on the o^rd. If it 
is punched, the 8 ^in E20.8 is ignored just^^like in F- fields. Thus the 
number -1.3^167^ might be punched in several different ways on input using 
the E20.8 and still enter memory with the same value. Here are some 
examples. / a. 

Different ways o-to punch the same value on 20 j^olumns of a card. ' 



20 columns 



□ □□□□□□□□ - . 131^672^01 

□ □□□□□□ □□□□-i,3li.67E0 

□ □□ □ -l.3l+^6700000Baoo 

□ □□ □ -13^^6700000. E- 09 

□ Dij □□ -13li-670000EbOO 



Orfly in the last example where no decimal .point is punched will tlte ^ 8^ 
*in E20.8 ^be used to tell where the decimal point should be. In all other 
cases th^ computer t^kes the number just as it is punched and interprets it 
the 'same way it would any .FORTRAN constant which. is written the, same way. 

One word of caution— always be sure-^^if you use thl-s-^pe of input>-to- 
have the exponent part of the. number occupy the right-most four or fewer 
columns. Leave no 4:railing blanks afte;r the' exponent ^s th^y aVe generally 
inteirpreted as zeros* For example: 



' 20' columns 



i3lf67EDlD 




Here most FORTRAN, implementations would interpret this value as if it 



were 



,0- 
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.13i^67E10 , Instead of 
.13467Er- as intended. -* 
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Now suppose we wanted to i^se ZERO; to find the root of the equation 
3x^ - 7x - 2 = 0 which lies l^tween 1 and 2. Then. we. should need to , 
define a function subprogram/iPUNCT which calculates the vallue of 
3x^ - 7x - 2. If we choose EPSI •= 10 , we would need to^ call ^ERO by the 

it v.* " ' 

^> ^ r ^ CALL ZER0^(1.0, 2.0, 1.0 E - i*; ROOT) 

THeT FORTRAN program could be written as shown in Figure F7-2. , 



CALL ZERO (1., 2., l.OE-i;) 

STOP 
END 

C - PLACE THE DEFINITION OF SUBROUJTINE ZERO, HERE 
FUNCTION FUNCT(X) ^ r T 
/FUNCT = (3.OJ^X*X-7.0hX-2.0^ ^' 
RETURN 
END 

Eito 

Figure F7-2r 



Using ZERO on several different functions ^ 

Ideally we would prefer to have written ZERO with a list of dummy arguments 
like those shown in Figure 7-6 of the main text. Unfortunately FORTRAN II 
processors do not permit statement label arguments, and only a few permit 
function name arguments. This cuts down the flexibility we would' like ,^but 
we can still find good uses for procediires like ZERO after some slight modi- 
fication. C , " ' ' . * 

V 

Suppose we wish to use the ZERO procedure with a series of functions* 

■ • * ' %^ 

Though desirable, this would, be difficult with ZERO as written^ As things 

stand now, we would have to. reproduce the 2»ER0 and ^^package^' it with each 
function separately. What we need is sorn^ means by which to, identifj^ each of 
a series of functions so that we can communio/S-f e to ZERO which of the functions 
it is to use, when called. Although we must always call our function FUNCT, 
we will now consider FUNCT to be a functiqn of two variables,^ *x an^i k, 

where k ,is really just an index to identify the particular function. 

' ' » ' h , ' 

The revised FORTRAN program which we will. call ZEROK follows in Figure F7r 3 ♦ 
Notice the fourth argument is K, an index that is used to identify the parti- 
cular fCmction to be called on< 



10 

11 



SUBROUTINE ZEROK(A,B,EPSI,k) 
XI = A 
X2 = B 

Yl = FUNCT(X1,K) 

IF(Yl - FU1(CT(X2.K)) 6, 5, 1^ 

if(y:^) 8, 7; 8 , , . , 

ROOT = XI • ^ 



ROOT = XM • * 

GO TO 12 

IF(n*FUNCT(XM,K)) I3, 10, Ik 



RETURN 
END 



Figure F7-3. The subroutine ZEROK. - It^s identical with 
. ' that of Figure F7-1 except for the opening 

declaration and statements that contain 
references to FUNCT. 



Figure FJ-k shows. how we might use -^EROK to find the root of 3x^ - Tx - 2 = 0 
between 1 and 2, the root of x^ . l,x^ + . x + 3 = 0 'betwe^ ll apd 
0 and'the root of, x = cosCx^ between 0 and 1. We choose EPSI = lO"^. 



CALL ZEROK (l., 2., l.OEA, l) ' 
CALL ZEROK (-1., 0., l.OE-i^, 2) 
CALL ZEROK (O., 1., l.OEA, 3) 
' STOP 

END X 
PLicE THE DEFINITION OF SUBROUTINE ZEROK HERE 
FUNCTION FUNCT(X,k) 
IF(K < 2) 1, 2, 3 

1 FUNCT = (3.*X^X^.)*X-2. 
RETURN. * ^ ' 

2 FUNCT = (((X-lt.)*X+7.)*X*X-l.)*X+3. 
RETURN ^ ^ ' 

3 FUNC? = COSF(X) . ^ 
RETURN V 
EJID . 

END 



Figtire . TJ^k. 
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F7-1 



The final step in the development of ZERO as a procedure is to give it* 
an alternate exit if the n;ethod proves inapplicable. We shall use the 
technique illustrated in Figures 5-26 and 5-27 of the main text. That is, ve 
shall employ another c^utput argument^ vhich vill be. assigned an index 

vajue in the procedure. Upon return to the calling program the val^ue Of L 
wilj> be inspected to see if an alternate exit has been implied. Once ve add 
this feature it viUl no longer be necessary to have the procedure do any 
printing. We can add ROOT as an output argumeni> and carry its value back to 
'the calling .program. The calling program can^ in turn^ take all responsibility 
"for printing results or" diagnostic messages. 

We shall call the new bisect ^procediire ZEROKL. To complete the analogy 
I'^*^ith^he parameter list in the, funnel of the zero procedure in Figure 7-6^ 
\we let the parameters for ZEROKL be: 




,The subprogram is given in Figure F7-5* The argument "L/ is set to 0 
to indicate a normal exit and is set to 1 to indicate the ^iltefn^te exit. 
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SUBROUTINE ZER0KL(K,L,-A;B,EPSI,R00T) 

n = A 

X2 = B 

n = FUNCT(X1; K) 

IF(Y1 - FUNCT(X2;K))6,5,if 

L = 1 

RETURN 

IF(Yl) 8, 7, 8 
ROOT = XL - ' ^ 
L = 0 

RETURN ' r 

ROOT = X2 

GO TO 12 ' * 

XM = 0.5*(X1 + X2) 

IF(ABSF(X1 - X2) - EPSI) 10, 11, 11 

ROOT = XM 

GO TO 12 

IF(Y1J^FUNCT(XM;K)) 13; 10, ik 
X2 = XM 
GO TO 6 
XI = XM 

GO TO 6 ' . ' 

END 



Figure F7-5. A FORTRAN analog to Figure 7-6 



Exercises F7-1 

~" ~~~~~ ' i * , } 

1. Write a FORTRAN program to solve all of the equations^.given"in Exercises 
7-1, Set C, Problem 1 of the main text. Use tHi indicated intervals- 
and the indicated error tolerancejs. , Use the subprogram ZEROK given in 
Figure F7-3» Rtm the program and cqmpare your results with the han4 
calciilated ones. Also run the program for 3 with EPSI = lO"^. 

2. Write and run a FORTRAN program to carry out the function evaluations 
needed in drawing the graphs in Exercises 7-1, Set A of the main text. 
You can use sane of the same FUNCTION subprograms needed in Problem 1. 



3. Write and run a FOPTIWI program to solve thp alley problem (Number 6) 

in Exercises 7-1, Set t) of the main text. Then solve the problem to the 
. nearest hundredth of a foot if the ladders are 2^83 and 19.14 feet 
long, and the crossover point is 7; 17 feet above the ground. 

4 -.5. For each of the flow chart solutions you prepared for Problems 1 and 
k, Exercises 7-1/ Set b, employ the ZEROKL subroutine developed in Figure 
F7-5 and write the nec^ssary companion FORTRAN programs (main programs 
and FUNCTION. subprograms). . ; . 
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F7-2 The area under a curve : An example, y = l/x between x = 1 and x = 2 

Since the area under the curve 'y = l/x is of interest in defining - 
logarithms we begin by writing a simple F0BTRAN 'program for the calculation 
of the approximate area under this curve between x = 1 and x = 2. This 
calculation will provide an approximation to jgn2. We assume that an error 
tolerance EPSI is read in from a card and that calculation of the app»oxi- 
'mate area is to be carried out by doubling the numbei' of subdivisions each 
time and terminating the calculation when the absolute value of the dj.ff^rence 
of two successive approx mat ions 1^ less than EPSI. The 'program in Figure F7-6 
follows closely the flow chart of Figure 7-l6. Remember that f(x) = l/x. 



105 



30 

8 
9 

106 



DIMENSION a?(lOl) 
bREAH 105, EPSI 
F0RMAT(E20.8) 

T(i) = o.5*(i.o + 0.5) 

N = 1 
S = 0 

lUL ='2*^*N-1 
FNP = 2*^*N 
DO 30K = 1, NL,2 
FK = K 

S = s+i.o/(i.o+fk/fnp) 

CONTINUE 

T(N+1) = 0.5*T(N) + S/FNP 

IF(ABS(T(N+1) - T(N)) - EPSi) 9,, 

N = N+1 

GO TO 2 . 

AREA = T(N) 

PRINT 106, AREA, 

F0RMAT( 8HaAREAD= E20 . 8 ) 

STOP 

END 

Figure F7-6 



ft* 



Note that we had to introduce FK as .a floating point value of K and 
FNP as a floating point value of the number 2^, in order to be able to 
calculate ^ho- desired terms in the calculation of S since FORTRAN does not 



permi1:,us^re^ mixS.modes in arithmetic expression 
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Exercises F7-2 



1. In the above program it is implicitly assumed that. the calculation will 



terminate' "before N exceeds 100. 
(a) Is- it possibly for .N to exceed 



100? 



(b) What would feappen if it failed to tdVmi^iate before N exceeds -' 
100? . , \ , . 

(c) Add some stateittents to the above progranv to protect against this 
undesirable event, even if the error toleraiiceU^ not 'satisfied/ 

; , 1 ^ 

Pi^nt out a message in this case indicating -failure to satisfy the 
error tolerance. ' ' 



Criticize the above program for inefficiency. <Revisei it to make it more 
efficient by following the Aow chart of Figure 7-17 . ^ Also include a 
"safety" termination if N exceeds 100. Run yoirr revised program 
using first 0.01 andrthen 0.001 as^ values for EPSI. 



Instead of terminating the calculation of the approximate area when two 
successive approximations differ in absolute value by l^ss than EPSI, 
we could terminate the calculation after a fixed finite number of 
approximations have been Qalculated. Rehrlse the pro^ram^^of the previous 
problem to read in an upper limit for the number of iterations to be 
carried out and then to terminate when this is reached. Run your 
program for N = 15. . - 

Tell how to revise the program given in this section so thelt the 
calculation. could be repeated* for a series of values of EP§I each 
of which is read in from a card. 



Write a FORTRAN program for the calculation described in Exercise 7-2, 
Set C, Problem 6, main text. Use f(,x) = l/x. Run your program and 
compare results using n = 5, 25, 75, 125, 200. 



i 
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F7-3 Area urider curte : the general case • j 

We now consider the general casfe of finding an approximation to the area 
luider a curve y = F(x), at)ove .the x-axis and "between the vertical lines 
X = A and x = B. In order tp make the program as useful as possible we could 
write it in the form of a function or a subroutine subprogram. We choose the 
former here. The function F(x) is assumed to be defined as a FORTRAN 
function called FUNCT(X)» An error tolerance* EPSI is given and we termi- 
nate the calculation when the absolve value of the difference of two succes- 
sive approximations is less than EPSI. We follow the flow chart "of Figure 7-20. 



1^0 



FUNCTION AREA(A,B/ePSI) . 
M = 1 ' 
li = B-A * 

OLAREA = .5*H*(FUNCT(a)+FUNCT(B)) ' 
"M = ^^M 
H = H/2.P 
S = 0. 

DO 1^0 K = 1, M, 2 

JFK = K 

S = S + FUNCT(A + FK*H) 
CONTINUE 

AKEA = .5*9IiAREA + H*S 
IF(^SF(AREA - OLAREA) - EPSi) 9, 8, 8 
RETURN 

OLAREA = AREA 
GO TO 3 ' 
END 

.Figure F7-7 




If we want to use t his function ^subprogram to calculate and print the , 



approximate area under the 6urve y = 1/x, above the , x-axis, ^and between 
the lines x = 1 and x = 2 we might-AAse a tolerance of ^ EPSI = 10 ^ andi' 
then we could write the following program.: ^ " ^ * • 



2;'^ AREA(1.0, 2.0, 1.0E-1+) 
\PRINT k,Z 
k FORMAT (7H AREA = E20.8) 
STOP 
END , 

PLACE THE DEFINITION OF FUNCTION AREA HERE^ 
FUNCTION FUNCT(X) 
FUNCT =^ l./X 
RETURN 
END 



r 
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Exercises 



1; " (a) Write a ffORTRAN function subprogram AREA2(A,B,N) which calculates 
an approximation to the area under the curve, y =^F(x), above the 
X-axis, and between the lines/ x ='A and x t; B and which u3es 
^ a subdivision of the interval ' (A^,3) into N equal parts. Follow 
• the flow chart drawn in Exekises 7-3, Problem 1, of the main text. 
Test your^prggram for y / sinx between x = 0 and ^x = n yith 
N = 5000. CHow does your result compare with the area of a s^mi- 
circle of diameter ir?) 

J[b) Use function AREA2 to print out a table of natural logarithm 
values from 1 ^through 51 in In tea-Vals of 5. 

/ 

' ^2. Tell how the function subprogram AREA(A,B,EPSI) of^this section may. 
be adapted ^ protect against the possibility of an .endless loop by 
causing termination of the calculation if the number of subdivisipns 
exdl^eds N. If the calculation is terminated in "this manner without 
satisfying the accuracy criterion, a message should be printed in 
addition to giving the approximations to the area. , * O' 

• Write a FORTRAN program which includes the subprogramSj^» of this 

section artd the subprogram AREA2 of Problem 1. Then uset7T555'sub- 

\ * ^ programs to calculate approximations to the areas described below. 

First vse 1, 2, k equal -subdivisions ef the interval' and then use an 
^ error tolerance of EPSI = lO"^. Of course you will need to supply the* 
necessary FORTRAN function subprograms to define the functions: 

^^y^^^^^ y = .^3te9/x, above x-axis, between x = 1 and x = 3. 
(True area is log 3) ^ * 

(b) Below y = 3x ' + 2x + 1^ above x-axis, between x = -2 and x = 2. 

(c) -Below y = x^,^ above y = x^ between x = ,1 and x = k. 

t 

h. ^Tite FORTRAN programs that may be used in calling on AREA (A, B, EPSI ) 
to compute an approximate value of n to four decimal placesi (See 
Problem 6, in Exercises 7-3 of th^ main'text.) 
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fl-h^ Simultaneous linear equations ; Jtoielc^pinp^ a systematic nfethc;^ of solution 

.In this^ section of_the main text we explained carefully/how to solve 

systems ,of two and three simultaneous equations. 'Exercises/ F7-^ provided 

examples o'f the method. You should now be ready to write jL simple F^pTRAN 
^ program- for the ^lution of two simultaneous equations in -^wo unknowns. 

Exercises F7-^ * * / 

' 1. Follow the flow chart drawn in Exercises 7-^+, Set B of the main text, and 
write a corresponding FORTRAN pl-pgram for the solution of two simultaneous 
equations in two unknowns: 

^11^1 ^12^2 = \ . 



2. Uf/e the program of Problem 1 to solve the following systems of equations . 
Vn the computer. Make a hand- calculated check of your comp^uter results. 
•For systems (f) and (g) slide rule accuracy is sufficient. 



(a) 








= 5 


(e) 


5x + y = 2 








+ 


y = 


= 1+ 


* 


.3x - hy = 1 




(b) 


hx 


+ 


3y 


= 5- ^ 


(f)_ 


3.12i+x^+ 5.375x2 


= -I.23U 




.> 




hy 


= 7 




10.2i+6x^ - 5-21i+Xg 


= 3.7ii+ 


(c) 


3x 




hy 


= 12 


(g) 


5.128x^ - 3.87i+X2 = 


12.ij.2 






+ 


6y 


= 3 




(d) 












3'.8l7x^ + 15.157x2' 


= 3..W 


2x 


+ 


hy 


= -7 










3x 


+ 


y = 


: 2 






* 
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FT -5 Simultaneous linear 'equations ; Gauss algorithm ' . ' ' , 

In describing the solution of three' equations in three -ynknowns we 
d'escribed^^h qf the'essential operations in 'turn aifc drew a flow chart for 
each,' It will be instructive to ,b\iil(fup the FORTRAN program in the same ■ . 
gradual fashion/ We begin by divj.ding the first equation through by a^^ ^ 
as described in Figure ?-2U, The corresponding FORTRAN statements would be 

DO 100, J = 2; '3 ' \ 

A(1,J) = A(1,J)/A(1,1) 



100 



CONTINUE 
B(l) ^ B(1)/A{1,1) 



The elimination of x^' from the i^^ row, i = 2, 3 is described in 
Figul-e .7^25 and the corresponding FORTRAN statements would be 



^0 



DO 200, J = 2, 3 

A(I,J) = A(I,J) - A(I,1) A(1,J) 
CONTINUE 

B(I) = B(I) -^A(I,1) * B(l) 



Next ^e have to divide the new second equation by a^^ and then 
eliminate* from the third equation. Following' these simple examples, you 
should h^ little trouble writing the FORTRAN that's equivalent to 
figures 7-27 through 7-30, 



* £xerc^<i^s F7-3 Set A 

1. Write the FORTRAN statements corresponding to the flow chart of 
ligtrfe 7-27. • ' . . 



2v 



Note the similarity between the statements of ^Problem 1 and those 
corresponding to Figure l-2k. .Write a single set of FORTRAN statements 
toj cover bolih ^ies b;;^ollowing Figure 7-28. (Remember tha| k = 1 

^ite the FORTRAN statements for the flow chart of Figure 7-30, Be sure 
*to take account of the possibility that k may exceed 2. ' 

Now write the FORTRAN that's. equivalent to Figure 7-33. * 



F7-5 



Nex"t we wan1r to carry out the back solution in order to obtain x^, Xg, 
X, / in tuf*n. This is described in the flow charts 'of Figures 7-3^ and 7-35«* 
The FORTRAN statements corresponding' to the latter flow chart, Figure 7-35 
'migh-t Ije: - ' 



1100 
900 



DO 900 I = 1, 3 / 
II = If-I ^ ^ 

x(ia) = B(ii) ; ■ 

IF(I - 1) 11, 900, 11 
IlMl = II -.1 
DO 1100 'J = 1, IlMl 
Jl J 

X(I1) = X(I1) - A(I1, Jl) * xUl) 
CONTINUE * 



This set of ..-statements appears to be somewhat more complicated than the 
flow chart of Figure One reasbn is that in the flow chart the c^ontrolled 

, variables i and J a^ decreasing" whereas FORTRAN requires that controlled 
-variaBles always increase- We take fcare of this by making the substitutions 



II' l^. - I and Jl = l^. - J, Another reason is that in most versions of 

DO loop i^ 'alw^s executed at least once. To avoid th±s diffi- 



FORi;^ a 

Culty we bypass the inner J)& loop if I = 1. • ' 

> . Now just as the^c«mplete flow chart of Figure 7-33 was built up frem 
partial flow charts, so we can build up the C9mplete. FORTRAN progi^am corres- 
ponding to Figure 7-36 from the partial FORTRAN programs which we have Just 
discussfed and^hich you^ have written in Exercise F7-5, Set A. 



^^^^^ * ^ , ^ 

v^:^ ^ Exercises F7-5 Set B * • >^ 

"-^ 

^r;^, 1. Write a con5)lete FORTRAI^program for the Gauss Algorithm given in 

• Figure.. 7-36. ( . ' [ 

/ , ' • ♦ * 

2. JE^u^ the above program on your machine and use phe program, to solve the 
systems of simailtaneoua^^ linear equations represented by the following 
array s . 



ERLC 



(a) 3x + i^y + z = -7 
2x + l^y z = 3 ' 
3x -..5y + 3z = 7 

•(b) X + ^ - 2 = 1^ 
3x - ?y + i^z = 1 
K - 3y - 2z = 7 



13J 



(c) , l^x . ^ . 3z = 7 

' 3x - 5y + 2z = 1 
^ &c + y'+ 2z = 1 

(d) 2x* - y +1Sz = 3 
3x -'W + i^z = 1 
X + ^ - 5z ^ 7 

130 * * 



3. Now use the above program^ on your machine to solve these systems of 
'equations » ^ 

(a) ^.Iklx^ + 2.419:^ - 3M9x = 4.219 ^ ^ 
_ 6.21HX' - 5-.678xp +-!^.271x, = -52.17 • 

t&^lyi^ + 5.761x2 + S^.'Sl^x^ ^27.11^ V * ' 

(b) 27.mx^^ ^3 AlTx^ - 3.^:^3. = 5.617 
31.468x^ + 3.^8x2 -f h.7l9x^ = 31.421 
11.121xj^ - 3.171x2 + 5.3l4x = -17.121 



Solution^ of n equatior;s in n . unknowns . " ^ . 

; The gerieraliaation to n equations is quite easy if we follow exkctXy 
the pattern we just used for 3 equations.. , 



Exercise F7-5 Set C ' , , 

Revise your FORTOAN p^rogram for the *Gaus^s Algorith^o handle n equations 
and n unknowns according to the p2a6cedure* fl^sw chart you prepared for * ' 
Problem^^2y'7-5 ^et A of the main t^xt. ^st t^e procedure (subroutine sub- 
program) using the ^ i| by 4 .system giy^n in Problem 5 of 7-5 Set-B in- 
the main text. The pro(jedure should be capable of handling up to I5 
equati6ns in»'15 unknowns. • . « r 



*Efercises F7-5 ' Set D 



1.^ ' In Exercises 7-5^ Set C of the main text you were asked to insert 
"partial pivoting" as a capability of your flo1?..chart for the Gauss 
procedure. Show the corresponding ^hanges necessary to thd FORTRAN 
subroutine called GAUSS^ which you prepared in the previous exercise. 



/ 

7 



^"'^Otnflete the modification of GAUSS and use the nev subroutine t(? solve 
the following systems of simultaneous linear equations. Compare results 
with and without partial pivoting. » ' . * 

(a) 3x2, - ^x^ = * 

5x^* + 15x2 - 3x2 = ^ 

(b) '^2x^ - 3x2 + ^x^ = 7 ' ' • 

^^1 " ^^2 -^^^3 ^ -^-^ 
2x-j^ - 7X2 - 12X2 = 1 
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alphanumeric datfe, kk 

alternate exits from subroutines, 113 

area under a curve - * 

' from X = 1 to X = 2, 12V 
general case, 126 - . 

arithmetic expression, 27 

array- 
input and oujiput, 72 
storage, 70 

assignment meaning when there are^ 
type differences, 35 ^ 

assignment statement, 26 

bisection process, II7 
blani^ spacesf, kl 

branching, 51 ^ 

card layout^ 5 
carriage control, h2 
character set, 8 * 
compiler program, 2 

composition of function references, 107 
compound conditions, 65 ^ ' >^ 
conditional statement, 5I 
contiauation, ko - . • 

continue statement, 8(1 , 



DIMENSION declaration, -y7l,> 75 

DO statement, 79 

"domain of a function, I05 

double subscripts, 75 . ^ 

E-fields, 118 ^ 
exponentiation, 31 

_ field codes, 17 

^ F-field code, 19 

I-field code, 19 
floating point, 8 
FORMAT statement, I8 ^ 
Fortran II, 1 - ' ' 

f unct ion 

range of*, 105 
function names, 11 

as'' arguments, 112 
function reference, 27 
function subprogram, 101 
functions, I05 ^ 

of more ths0t one* argument, I06 

GSUss algorithm, 129 ^ 

GO TO statement, U 

greatest integer function, 30 

H-fiel<ib, 59 ^ . 1 



identifying remarks in FORTRAN output, 57 
IF statement, 51 

' rules for, 5© 
implied DO loops, 90 / ^ . , 

input -output statements, I5 
integer, 8 

division, 30 
variable, 11 
iteration, 79 

• 

labels ,\ 8, 11 

length of a statement, Uo ' 
local variables, 107- 
logical expressions, 68 
looking, 79 

mixed mode in arithmetic expressions, 29 
multiple branching, 65 

nested DCl loops, 96 
nested Jjnplied -DO loops, 97 
non-local variables, 107 \ 
number types, 11 

in arithmetic expressions, 27 
numerical constants, 8 ^ ' 

operator symbols, lU 
• order ^of computation, 3U 

predefined mathematical' functions, 12- 
PRINT ^statement, 22 
printer carriage control, ^2 
PROCEDURES, 109 

READ statement, 16' ' , ^ 

real variable, 11 

RETURN statement, 102 ^ ' 

simultaneous linear equations, 128 
source program, 2 

statement labels as argiuii^ts, 112 
storage of doubly-^sub^criptfed arrays, 75 
strings, 'llU 
subprograms, 101' 

parameter li^st of/ I06 
subroutines) 102 

branching from, 112 
' ^ call of, 110 % 
subscripted variables, 6$ 
symbol Manipulation ^In FORTRAN, 11^ 

table-look-up, 89 
target program, 2 
terminal statemer 
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unary minus, ^ > 



variables, 11 ^' \ 

X-field, 91 

ZERO subroutine, ll8 - 
ZEROK subroutine, 121 / 
ZEROKL subroutine, 123 
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